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The reaction Cr® (p, m) Mn® at 6.6-Mev proton energy produces isomers with period 21.3 
min. and 6.5 days. These are positron emitters with maximum positron energy of 2.2 Mev (21 
min.) and 0.77 Mev (6.5 day). The gamma-ray spectrum for each isomer has been investigated. 
There are 2 gamma-rays per positron in the short period activity. The gamma-radiation in the 
6.5-day activity is so intense that 95 percent of these disintegrations must be due to K- 
electron capture; this is followed by emission of probably two gamma-rays. In the remaining 5 
percent the positron emission is followed by three gamma-rays. The threshold for the 21-min. 
activity and the yield for both periods have been determined. 





AMPLES of metallic chromium have been 

bombarded with 6.6-Mev protons produced 
in the University of Rochester cyclotron. Posi- 
tron activities of two periods were observed and 
studied: 21.3 min. and 6.5 days. These periods 
have been observed in other reactions, and 
Livingood and Seaborg,' assign them to isomers 
of Mn®. This conclusion is consistent with 
results to be presented here. 


I. MEASUREMENT OF PERIODS 


The bombarded samples of metallic chromium 
weighed about 5 g. The surface of the metal was 
cleaned by grinding on an emery wheel. The 
initial activity for the 21-min. period was 67 
microcuries per microampere of protons for 
infinite duration of bombardment at 6.6-Mev 


*A preliminary report of this work was made at the 
Washington meeting of the American Physical Society, 
December 1938. See Phys. Rev. 55, 604 (1939). 

tA omg of work done at the University of Rochester 
during the summers of 1938 and 1939. 
assy Livingood and G. T. Seaborg, Phys. Rev. 54, 391 


proton energy. The initial activity for the 6.5-day 
period was 4.4 wC/ya under the same conditions. 

The two activities mentioned were identified 
chemically as Mn. A layer a few tenths of a 
millimeter thick was dissolved off of the Cr target 
by immersion in hot HCl. To the acid solution 
MnCl, was added for a carrier. NaOH was added” 
until the solution was basic, and it was then 
oxidized with Na:O2. The solution was boiled, 
diluted, and filtered. The precipitate showed the 
21-min. and 6.5-day activities. 

The periods were measured with an ionization 
chamber containing Freon-12 at 2 atmospheres 
pressure and an FP-54 amplifier arranged for 
automatic photographic recording of the data.* 
Figure 1 shows the decay curve of a Cr target 
which had been exposed to 6.8-Mev protons at 
0.07 wa for an hour. From this curve we determine 
the shorter period as 21.3 min. Subtracting the 
long period activity from the total, we see that 
the two periods mentioned account for practi- 
cally all of the observed activity, and there is no 


2S. W. Barnes, Rev. Sci. Inst. 10, 1 (1939). 
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Fic. 1. Decay curve of Mn®. Cr target bombarded at 
6.8 Mev for 60 min. at 0.074a. Measurements start 4 min. 
after end of bombardment. 


evidence for a 46-min. period, which activity 
occurs at lower proton energies.** We have also 
found the 310-day period assigned by Livingood 
and Seaborg! to Mn®™, but have not studied it. 
Figure 2, showing the decay curve of a Mn 
precipitate several days after bombardment of 
the parent Cr, gives the half-life as 7.4+1 days. 
The Cr target was bombarded for 45 min. at a 
current of 0.5 ya of 6.6 Mev protons. This 
period is slightly greater than Livingood and 
Seaborg’s value of 6.5+1 days, and since the 
310-day activity was not subtracted, it seems 
reasonable to assume 6.5 days to be the more 
precise value. 


II. PosiITRoN ENERGY 


The energies of both sets of positrons were 
measured by means of a cloud chamber in a 
magnetic field. Histograms are shown in Fig. 3 
and Fig. 4. The cloud-chamber tracks were 
mostly positrons, distributed in the usual beta- 
ray spectrum. The upper energy limits of the 
positron spectra are 2.2 Mev (21 min.) and 0.77 
Mev (6.5 day). The energy 0.77 Mev agrees 
fairly well with the value (0.7) given by Livin- 
good and Seaborg.! They do not report the energy 


3 Barnes, DuBridge, Wiig, Buck and Strain, Phys. Rev. 
$1, 775 agen: $38) ridge, Feeney Buck and Strain, Phys. 


Rev. 53, 447 
« Delsasso, een, Sherr and White, Phys. Rev. 55, 
113 (1939). 
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of the 21-min. positrons. We have measured the 
range in Al of these particles and find it to be 
0.973 g/cm?. Using the relation® 


energy = (range (g/cm?*)+0.165)/0.536 Mev, 


we get a value of 2.12 Mev. 


III. ExcitTaTION FUNCTION 


By bombarding a pile of thin films of chromium 
the threshold for the 21-min. activity was de. 
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Fic. 2. Decay curve of Mn precipitate. Cr target bom- 
barded at 6.6 Mev for 45 min. at 0.54a. Measurements 
start two days after end of bombardment. 


termined. The intensity of the 6.5-day activity 
was far too small for such a measurement. Cr 
films about 0.184 thick were evaporated onto 
2.5u Al foil. The Cr side of the foil was exposed 
to the proton beam so that the Mn recoil nuclei 
could not get away. The stopping power of the 
pile of eight foils was measured by placing it in 
the proton beam and observing the decrease in 
range of the protons. The stopping power of the 
pile was 6 cm air equivalent. Protons of this 
energy produce no activity in Al of such a 
period and intensity to be troublesome even in 
comparison with the activities of these thin 
targets. The incident proton energy at the first 
foil was 6.6 Mev. There was an 0.05 Mev energy 


5 E. E. Widdowson and F. C. Champion, Proc. Phys. Soc. 


50, 192 (1938). 
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loss in each Cr+Al film. The decay curves for 
the various films are shown in Fig. 5. The initial 
activities are plotted against the proton energy 
in Fig. 6. We find the threshold for the 21-min. 
activity to be 6.3 Mev. The actual reaction 
energy (threshold corrected for nuclear recoil) is 
6.3(52/53) =6.2 Mev. All we can say about the 
threshold of the 6.5-day activity is that it does 
not differ from 6.2 Mev by more than a few 
tenths of a Mev; at energies lower than 6.3 
Mev the 46-min. and 310-day activities are so 
great compared to the 6.5-day that a more 
precise determination of this threshold is hardly 


feasible. 


IV. GamMa-Rays 


Gamma-ray absorption curves are shown in 
Fig. 7. The first part of the curve for the 21-min. 
period has a slope which probably corresponds to 
annihilation radiation. In the 6.5-day activity 
the intensity of annihilation radiation is negligible 
compared to that of higher energy radiation. 
Using data given by Gentner® we estimate the 
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Fic. 3. Distribution of positron tracks from 21-min. Mn®. 


average gamma-ray energies as 1.1+0.2 Mev 
(21 min.) and 0.9+0.2 Mev (6.5 day). Since both 
gamma-ray spectra are probably complex these 
energy values are very uncertain. 

The source used in plotting the gamma-ray 
absorption curves (Fig. 7) was a thick Cr target 
covered with an Al sheet. Thus positrons 
emitted in any direction were stopped, and for 
each positron there must be two quanta of 


* W. Gentner, J. de phys. et rad. 6, 274 (1935). 
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annihilation radiation. The ionization current 
due to the hard component is twice that due to 
the soft (from Fig. 7), and the ionization per 
quantum of the hard radiation will possibly be 
greater than that of the 0.5-Mev radiation. This 
suggests that there may be about two 1.1-Mev 
quanta per positron. The assumption that there 
are two 1.1-Mev gammas per positron is also 
necessary to account for the observed threshold 
of 6.2 Mev, since the upper limit of the positron 
spectrum is only 2.2 Mev, leaving 2.2 Mev for 


.the gamma-ray energy. 


We can estimate the intensity of the gamma- 
radiation in the 6.5-day activity by calculating 
the ratio ['=y/(7+8) of ionization produced by 
gammas to total ionization, the former being 
measured with an Al absorber thick enough to 
stop all positrons. 

These ratios are 0.025 (21 min.) and 0.24 (6.5 
day). For the 21-min. activity y/8 =2.5 percent, 
of which two-thirds is for hard gamma-rays. 
For these only, y/8=1.6 percent. For the 6.5-day 
activity y/8=32 percent. Since the sensitivity of 
the ionization chamber does not vary greatly 
with energy in this region, and since the gamma- 
rays have nearly the same absorption coefficients, 
the number of gammas per positron in the 6.5- 
day activity is 32/1.6=20 times as many in the 
21-min. activity, which means 40 gammas per 
positron. 

We conclude that the 6.5-day isomer decays 
not only by positron emission, but also by 
K-electron capture followed by emission of 
gamma-radiation. Presumably 3 gammas accom- 
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Fic. 4. Distribution of positron tracks from 6.5-day Mn®. 














pany the positron, leaving the other 37 for K 
capture. Assuming 2 gammas per K capture, 
then there are about 18 K captures per positron, 
or positron emission occurs in only 5 percent of 
the disintegrations. 

We have looked for x-rays emitted by the 6.5- 
day isomer. There is evidently a soft radiation, 
which we were unable to identify, apparently 
because of the absorption in air. 


V. COINCIDENCE COUNTER MEASUREMENTS 


The data presented so far show that there 
might be some interesting features of the dis- 
integration of the 6.5-day isomer, and it seemed 
at least possible that these could be investigated 
by means of coincidence counters. Thick-walled 
alcohol-argon filled counters were used, and an 
amplifier with a resolving time (7) of (9.17 +0.25) 
X10-* min. The details of the counters, ampli- 
fiers, scaling, circuit, and the technique of 
measuring 7 will be described elsewhere. Counter 
data were taken for both short and long periods. 
In each case, measurements were made of single 
and coincidence counting rates with the thin 
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Fic. 5. Decay curves of thin films of Cr bombarded with 
protons. 
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(0.5 mm) source in both positions A and B of 
Fig. 8. In position A the source is placed between 
two Al slabs 0.25 in. thick which serve to stop 
all positrons. In position B ideally only haif of 
the positrons produce annihilation radiation 
which affects the counters. Actually, the number 
of annihilation quanta counted in position B js 
more than half of the total number, and this 
ratio must be ascertained by separate meas- 
urements. 
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Fic. 6. Excitation curve for the formation of the 21-min. 
Mn®. Threshold estimated at 6.3 Mev. 


If we let N,=single counting rate, position A; 
N2=single counting rate, position B; N3=coin- 
cidence counting rate, position A; N4=coinci- 
dence counting rate, position B; N=number of 
disintegrations per minute; »=average number 
of gamma-rays per positron (exclusive of anni- 
hilation radiation); S,=average gamma-ray 
sensitivity of the counter including the solid 
angle; S,=annihilation radiation sensitivity of 
counter; and K=a constant to be determined 
later (with the source in a vacuum, K = 1), we can 
write the following equations for the counting 
rates: 


N,=NnS,+2NSu,, (1a) 
N2=NnS,+KNS,, (1b) 
N3= NS,?n(n—1)+4NnS,Sa, (1c) 
Ni= NS,?n(n—1)+2KNnS,S.. (1d) 


To determine K, separate measurements were 
made using a source of Zn®. This element emits 
2.3-Mev positrons (nearly the same positron 
energy as the 21-min. Mn®) and no gamma-rays, 
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Fic. 7. Absorption curves in Pb of gamma-radiation from 
Mn®, both isomers. 


so that n»=0. By measuring single counting rates 
in positions A and B we can solve Egs. (1a) and 
(ib), obtaining K=2N2/N,. The result of such 
measurement is K = 1.06. 

We can now solve Eqs. (1) for m, N, Sa, Sy). 
We are at present interested in the value of m: 








[ Fatale el 
n= —_ ° ° 
2N2—KN,; N3—N, 


In measuring N,; and N,2 the background of 
the counters must be subtracted as well as the 
activity of the 6.5-day Mn® present. The single 
counting rates of the two counters were nearly 
the same, and for NV; and N2 averages of the rates 
for the two counters were used. Actually, N; and 
Nz could not be measured until the coincidence 
counting rates were too small to measure. In 
measuring NV; and N, the procedure was to start 
with a sample which gave about 1000 coincidence 
counts per min. in position A. Coincidence decay 
curves in positions A and B were then plotted 
(with cosmic-ray background subtracted). Since 
N,, N2, and r are known, the accidental counting 
rates are known, so one can calculate the true 
initial coincidence counting rates. The first part 
of the data (counting rate 1000 coincidences per 
min.) is useless, for these are mostly accidental 
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coincidences. At best, the method does not 
yield data of high precision with a source which 
has a 21-min. period, for the interval of time in 
which coincidence counts can be made is too 
short. The data for 21-min. Mn® are as follows: 


N, =395500+120 (counts/min.) 
N2=305500+90 

N;=360+17 (counts/min.) 
N,=231+11. 


The fractional probable errors in N; and N-~ 
were computed in the usual way: 0.67 divided by 
the square root of the total number of counts. 
Since N; and N, are differences of measured and 
accidental coincidence rates, the probable error 
was computed for each measured point, and the 
values of N; and N, are means weighted accord- 
ing to the inverse square of the probable error. 
Using Eq. (2) we find m=1.26+0.17. This 
suggests that positron emission is followed some- 
times by the emission of two 1.1-Mev gamma- 
rays, and at other times by one 2.2-Mev gamma. 
Presumably the high energy gamma-radiation is 
not very intense, for it did not show up at all 
in the absorption curve (Fig. 7). There is some 
further evidence for the higher energy gamma- 
rays in cloud-chamber pictures which show pair 
production by this radiation in a 2.5-mil Pb foil. 

The 6.5-day. isomer, which disintegrates 95 
percent of the time by K capture, can be investi- 
gated by comparison of the ratios N3;/N, and 
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Fic. 8. Geometry of counters and source. 


N,/ Nz for the two isomers. Taking averages of 
the data for each isomer, we have 


ISOMER N;3/Ni Ni/N; 
21-min. (10+0.5)x10~ (7.540.3) x 10~ 
6.5-day (3.5+0.08)k10-* (2.9+0.1)«x10—. 


Thus the coincidence rate for this activity is 
less than half that for the short period activity. 
Remembering that in the 21-min. activity there 
are 3.3 gammas per disintegration (including 
annihilation radiation) and never more than 4, 
it appears that for the majority of 6.5-day dis- 





integrations there are at most 2 gamma-rays. 
For disintegrations in which there is positron 
emission the energy available for gamma-radia- 
tion is 3.6 Mev. If some of this is the same as 
the 2.2-Mev radiation of the 21-min. isomer, 
the energy of the additional gamma-ray is 1.4 
Mev. The lower absorption curve of Fig. 7 gives 
the energy of the 6.5-day gamma-radiation, as 
about 1 Mev. If there were only one gamma-ray 
of this energy per K capture, an implausibly 
large amount of energy would have to be 
assigned to the neutrino. It is more probable 
that there are two gammas per K capture, one 
of them being possibly of about 3.2-Mev energy. 


VI. YIELD 


The initial activity of the 6-day isomer of 
4.4 wC/ya is given in Section I. The initial 
activity due to positrons only is 0.75X4.4 
=3.3 uC/ya. If all the disintegrations were by 
positron emission the initial activity would be 
3.3/0.05=66 wC/ya. Using the calibration data 
of the ionization chamber, we get the yields for 
each isomer at 6.6-Mev proton energy: 


0.40 x 10-* 
0.68 x 10-* 


radioactive atoms per proton 
radioactive atoms per proton. 


21 min. 
6.5 day 


VII. 


The data on the two isomers of Mn® are 
summarized in the following table: 











DISINTE- 

ISOMER GRATION 8+ KE THRESHOLD YIELD 
21.3-min. Br 2.2 Mev| 6.2 Mev | 0.40x10-¢ 
6.5-day Bt 5% | 0.77 0.68 x 10-* 

K 95% 




















The results of the gamma-ray measurements 
are somewhat uncertain, but all available data 
are consistent with the following possible scheme: 
the 21-min. isomer disintegrates with the emis- 





934 ARTHUR HEMMENDINGER 
































MEV Mn*? 
5.4 2! min. 
K 
bai wl / {> 6.5 days 
ry 
” 
» ‘ 
n 
% . 4 
¥ Q 
3 
— ao = —«¢ 
(3-4)-— 4 ’ 
% 
v 
“a 
g 
2.2 
1.1 
1@) d \ cr” 
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sion of a positron followed by two 1.1-Mey 
gammas or one 2.2-Mev gamma. The 6.5-day 
isomer, disintegrating by positron emission emits 
3.6 Mev of energy, probably in 3 quanta; dis- 
integrating by K capture, it emits 4.2 Mev of 
radiation in two quanta, one of which must be 
about 1 Mev. An energy level diagram corre- 
sponding to this scheme is shown in Fig. 9. 

In conclusion I wish to thank Professor L. A. 
DuBridge for making available laboratory facili- 
ties at the University of Rochester, and I am 
indebted to the whole staff at Rochester for 
their generous assistance, in particular Dr. S. N. 
Van Voorhis, Dr. S. W. Barnes and Dr. C. V. 
Strain. I wish also to thank Dr. L. M. Langer of 
Indiana University for many helpful discussions 
regarding the counter work. 
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Electrostatic Deflection Studies of Alpha-Particles; the-Alphas from Li‘(p, «)He* 


LeonarD C. MILLER 
Ryerson Physical Laboratory, University of Chicago, Chicago, Illinois 
(Received October 8, 1940) 


A technique is described for studying nuclear reactions in which short range alpha-particles 
are produced with a yield that is too small for the direct application of the electrostatic deflector. 
In preliminary experiments a strong polonium source is used to supply alpha-particles to the 
deflector, and a group of these alphas, homogeneous in energy, passes through the deflector and 
enters a variable pressure absorption cell and ionization chamber. Cut-off pressures determined 
for various alpha-particle energies result in a calibrated cell which can then be used to determine 
the alpha-particle energies in nuclear disintegrations. Application of this method to Li*(p,«)He? 
results in a Q value of 3.94+0.08 Mev. The calibration, extending in alpha-particle energy from 
1 to 1.8 Mev, is also used to determine the slope of the range energy curve. The average value of 
the slope found for this region is 200 kv per mm range. 





HE direct measurement of the energies of 

short range disintegration particles by de- 
flection in an electrostatic field, introduced by 
Allison! and his group, is suitable for high yield 
reactions. For the study of reactions in which 
short range alpha-particles are produced, but 
where the yield is too small to apply the electro- 
static deflector directly, the following technique 
has been devised. A variable pressure absorption 
cell, suitable for the particles, can be calibrated 
by using alpha-particles from a strong natural 
source. These particles are slowed down in air 
and then allowed to pass through the electro- 
static deflector which measures their energies. 
This cell can then be used to determine the 
energies of short range alpha-particles, since the 
solid angle through which the cell will accept 
particles is much greater than that of the 
deflector. 

This paper describes such a calibration and 
the application of the method to the alpha- 
particles from the disintegration of Li® by 
protons. The latter part of the paper shows how 
this calibration can be used to determine the 
slope of the range energy curve for alpha- 
particles in the energy region 1—-1.8 Mev. 


APPARATUS 


The electrostatic deflector and its accompany- 
ing high voltage source have been described in 
previous papers*? from this laboratory. The 

1S. K. Allison, L. S. Skaggs, and N. M. Smith, Jr., 


Phys. Rev. 54, 171 (1938). 
*L. S. Skaggs, Phys. Rev. 56, 24 (1939). 
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entrance slit to the deflector was 3 mm wide 
and 1.3 cm high, and was located 8.9 cm (the 
focal distance) from the edges of the cylindrical 
plates of the deflecting electrostatic field. These 
plates were spaced 0.635 cm apart and had an 
average radius of curvature of 25.400 cm. They 
subtended an angle of 90 degrees at their center 
of curvature and were 8.26 cm high. A collodion 
foil of 1.5 cm stopping power was placed over 
the entrance slit to seal up this end of the 
deflector. 

A strong polonium source’ deposited on silver, 
with cross section the same as that of the en- 
trance slit, was mounted on a screw arrangement 
in front of the entrance slit. Movement of the 
source with respect to the slit changed the 
portion of the range of the alpha-particles ex+ 
pended in air. 

Figure 1 shows the arrangement of the variable 
pressure absorption cell and ionization chamber 
which the particles entered after passing through 
the electrostatic deflector. The opening to the 
absorption cell, B, 3 mm by 2 cm, served as the 
exit slit of the electrostatic deflector. Over B was 
soldered a metal screen, 0.15 mm thick, with 
about 650 holes per square cm, each hole being 
about 0.2 mm in diameter. The screen was 
covered on the higher pressure side by a collodion 
window with an air equivalent of about 1.7 mm. 
The brass disk containing the opening B was 
seated accurately against shoulders on the cylin- 

*The author is indebted to Miss M. R. Jones of the 


Chemistry Department for preparing the strong RaD 
solution used. 
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drical body of the absorption cell. It was held 
in place by a collar and sealed by ‘a rubber 
gasket. The construction was such that the 
length of the cell was independent of the tight- 
ness with which the screws were drawn up to 
make the seal. The opening A, 6 mm by 2.3 cm, 
was covered with a similar metal screen but with 
no collodion film. Thus the pressures in the 
ionization chamber and absorption cell were the 
same. The absorption cell had a measured depth 
of about 0.59 cm and the ionization chamber of 
about 0.45 cm. The collection voltage of 530 
volts for the ionization chamber was applied to 
the outside brass case as indicated in the figure. 
The air admitted into the cell was freed of 
moisture and carbon dioxide by passage through 
tubes of calcium chloride and sodium hydroxide. 
Its temperature was measured by a thermometer 
whose bulb rested on the cell. The pressure was 
measured in an open tube mercury manometer. 
A correction was made for the temperature of 
the mercury. A volume of several liters was 
placed on the line to stabilize the pressure. 

The linear amplifier used in this work was of 
the resistance-coupled type with a 954 acorn 
pentode in the first stage. The output of the am- 
plifier was fed intoascale-of-sixteen circuit similar 
to that described by Stevenson and Getting.‘ 

For the study of the Li® reaction the high 
voltage set used was the Cockcroft-Walton circuit 














Fic. 1. Arrangement of the absorption cell and ionization 
chamber for taking number vs. aes curves of particles 
which von oy ssed through the electrostatic deflector. A pair 
of brass r ods (not shown 5 keeps the sylphon extended when 
the deflector i is evacuated. 


*E. C. Stevenson and I. A. Getting, Rev. Sci. Inst. 8, 
414 (1937). 
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described by Hatch.’ The arrangement of the 
target chamber and absorption cell-ionization 
chamber unit was similar to that used in a 
previous study‘ of this reaction in this laboratory, 


CALIBRATION OF THE ABSORPTION CELL wity 
THE ELECTROSTATIC DEFLECTOR 


The performance of the variable pressure 
ionization chamber has been described in a 
previous communication from this laboratory.‘ 
In our work the bias on the counting circuit for 
the main calibration curve was 13 volts, and 
in a manner to be described later it was deter- 
mined that a particle must produce about 4100 
ion pairs in the ionization chamber in order to be 
counted. There is a critical operating pressure P, 
above which the penetration into the chamber 
(in terms of air equivalent at standard condi- 
tions) required to produce a count should, 
theoretically, be independent of the pressure. 
P.. is the pressure at which the particle can pro- 
duce the number of ions necessary for a count 
before it reaches the end of the chamber. The 
reason for the independence of pressure can be 
briefly stated as follows. The actual length of 
path of a particle before it loses its ionizing power 
is inversely proportional to the pressure whereas 
the specific ionization is directly proportional to 
the pressure. ‘ 

The formula for the range of a particle as 
measured by the cell can be written 


R=(L/76)(288P/T)+C, (1) 


where L is the length of the absorption cell, 
T the absolute temperature and P the pressure. 
The first term represents that part of the range 
expended in the absorption cell, reduced to the 
standard conditions of 76 cm Hg pressure and 
15°C. C represents the air equivalent of the 
collodion window, assumed to be constant, plus 
that of the penetration into the ionization 
chamber. This latter part, then, is also constant 
so long as the pressure is above P,, which can 
be calculated from 


De 
No= f I(x)dx. 
0 


Here No is the number of ion pairs required to 


5G. T. Hatch, Phys. Rev. 54, 165 (1938). 
*G. J. Perlow, Phys. Rev. 58, 218 (1940). 
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Fic. 2. Curves for determining cut-off pressure as a function of the energy of the alpha-particles. 


produce a count, and is equal to 4100 for a 13- 
volt bias; I(x) is the specific ionization function ; 
D.=DP./76 where D is the linear depth of the 
ionization chamber. If we use the curve of 
Holloway and Livingston’ to represent the func- 
tion I(x) we find that P.=27.0 cm Hg. For this 
calculation the maximum ordinate of the specific 
ionization curve has been set at 60,000 ion pairs 
per cm. Only those number-pressure curves which 
did not start ‘to drop off until the cell pressure 
was above 27.0 cm Hg were used in this work. 
The quantity Ny» was determined in the 
following way. A thin polonium source was 
placed so that its alpha-particles passed through 
the absorption compartment and entered the 
ionization chamber near the end of their range. 
The pressure was then lowered until rough 
observations with the oscilloscope and counting 
circuit indicated that the particles traversing the 
ionization compartment were producing approxi- 
mately the critical number of ion pairs corre- 
sponding to the selected bias. With the pressure 
held at this value, the distance of the polonium 
source was readjusted until the pulses showed a 
maximum height. Then an accurate measure- 


™M. G. Holloway and M. S. Livingston, Phys. Rev. 54, 


18 (1938). 


ment of the cut-off bias was made by counting 
at this pressure and varying the bias to get a 
number bias curve. Such curves were taken at 
various pressures. Let the pressure associated 
with a given cut-off bias be Po. If the linear depth 
of the ionization chamber is D, then its equiva- 
lent depth at Py) is Do=P .D/76. The critical 
number of ion pairs No corresponding to a given 
bias may then be found from 


ve-| f az (=) 


where @ is to be chosen to give the greatest 
value for the integral. 

Figure 2 shows a set of curves taken with the 
electrostatic deflector. Each curve was taken at 
a fixed voltage on the deflector plates, with the 
ordinates giving counts per second and abscissae 
the pressure in the cell as read off the manometer. 
The energy of the alpha-particles passing through 
the deflector is given! by E=40.00V where V is 
the deflecting potential in volts. An end effect 
correction® equal to —0.0093E must be added on 
to E. The dispersion is seen to be quite high, 


max 


8S. K. Allison, L. S. Skaggs and N. M. Smith, Jr., 
Phys. Rev. 57, 550 (1940). 
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a 100-kv energy change resulting in a pressure 
change of about 5 cm Hg. Most of the straggling 
is range straggling, with small contributions due 
to noise straggling and ionization straggling. 
The electrostatic deflector itself contributes a 
negligible amount, since calculations of its re- 
solving power show! that with 3-mm slits the 
energy spread of the particles passing through it 
is of the order of 1 percent. 

Since the particles leaving the deflector are 
essentially homogeneous in energy, these curves 
can be treated as ordinary number vs. range 
curves for a group of alpha-particles all of which 
have the same initial energy. In our case the 
pressure at half-maximum can be used as cut-off 
pressure for the particle of mean range in each 
group; the values, uncorrected for temperature, 
are indicated in Fig. 2. Curves for the same 
alpha-particle energy taken on different days 
showed variations in cut-off pressures, amounting 
in the extreme case to about 1.8 cm Hg. These 
variations were attributed to such factors as 
changes in the gain of the amplifier or in the 
constants of the scaling circuit. 

Sets of curves were taken on eight different 
days. On all but two of these days one of the 
curves was for an alpha-particle energy equal to 
1.437 Mev. The average value of 288 P/T, where 
P is the cut-off pressure and T the absolute tem- 
perature, for these curves was 44.72 cm Hg. On 
the assumption that the gain of the amplifier, 
constants of the scaling circuit, etc., might be 
different on different days but essentially con- 
stant during the course of one day’s runs, all 
cut-off pressures for the 1.437-Mev curves were 
adjusted to 44.72 cm Hg, and the pressures for 


TaBLe I. Values of (288/T)P for various alpha-particle 














energies. 
a-PARTICLE UNADJUSTED VALUES 
ENERGY OF (288/7T)P 
(MEv) (cm Ho) ADJUSTED VALUES 
1.730 60.62 60.68 
1.671 57.13 56.79 
1.553 51.10, 51.05 50.42, 50.71 
1.496 47.11* 47.50* 
1.437 44.24,* 44.27, 45.40 44.72 
45.06, 44.79, 44.06 
1.379 40.74* 41.22* 
1.333 37.51*, 39.31, 39.37 37.99*, 39.76, 38.69 
1.205 33.62, 34.04, 33.62 34.07, 34.07, 32.94 
1.086 28.65 28.68 








* Old 954 tube. 
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Fic. 3. Energy vs. cut-off pressure. The open circles 
represent points determined with the old 954 tube in the 
first stage of the amplifier. 


the other curves for each day’s runs were shifted 
by an amount equal to the shift required in the 
1.437-Mev curve of that day. Table I shows the 
unadjusted values when corrected for tempera- 
ture and the adjusted values. 

Figure 3 shows a plot of energy vs. cut-off 
pressures. After the three points indicated by 
open circles were taken it was necessary to re- 
place the 954 tube in the first stage. At this 
point the gain of the amplifier was adjusted so 
that the cut-off pressure for 1.437 Mev occurred 
at about the same point as with the old tube. 
No changes in the amplifier or in any other 
parts of the circuits were necessary after this, 
An examination of the cut-off pressures shows 
that the curve of Fig. 3 would not be altered 
essentially if the unadjusted points were used. 


. Use or CELL For Stupy oF Li®(p,a)He® 


The absorption cell-ionization chamber device 
was placed opposite a thin target of lithium 
fluoride to detect the alpha-particles at right 
angles to the incident proton beam. The auxiliary 
apparatus, such as the amplifier and scaling 
circuit, was the same as that used in the calibra- 
tion experiments. Moreover, care was taken to 
keep unchanged all the wiring connecting various 
circuits together, such as the leads from the out- 
put of the amplifier to the scaling circuit, etc., 
since changes in these might affect the operation 
of the circuits. 

A metal disk with a 3-mm hole in the center 
was placed over the absorption cell slit; this 
opening and the target spot (5.4-mm diameter) 
served as the limiting aperture widths with a 
separation between them of 4.0 cm. The proton 
beam of about 5 microamperes was fed into a 
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current integrator. The bombarding voltage was 
342 kv. This is rather close to the observed 330- 
kv resonance in F!*, but comparison of the curves 
obtained with those obtained in a previous study® 
of this reaction in this laboratory, made at a 
bombarding voltage of 370 kv, shows that the 
alpha-particles from fluorine were avoided. The 
curves obtained at 370 kv, which is farther 
above the 330-kv resonance, gave the same ratio 
of counts on the plateau to counts on the base 
line as that shown in our curves. 

In order that the cut-off pressure should fall 
within the pressure range for which the cell was 
calibrated it was found necessary to raise the 
bias to 21.0 volts. Figure 4 shows a plot of the 
number of counts vs. pressure, the number of 
counts per point being about 220 on the plateau. 
The base line represents the long range alphas 
from Li’ in addition to the He*® particles and 
background noise. Though the points are some- 
what scattered, the pressure at half-maximum is 
fairly well defined at 58.5 cm Hg. When corrected 
for temperature (28.9°C) this becomes 56.0 
cm Hg. 

In order to determine the energy corresponding 
to 56.0 cm Hg at 21.0-volt bias it was necessary 
to take more data with the electrostatic deflector. 
The energy vs. cut-off pressure curve at this 
higher bias should be parallel to the one at the 
lower bias, since it is only the value of C in 
Eq. (1) that is changed. To check this, two 
points were determined at the 21.0-volt bias, 
one with an alpha-particle energy equal to 
1.437 Mev, the other 1.612 Mev. These are 
plotted in Fig. 3. It is seen that the line they 
determine is parallel to the 13.0-volt bias line. 
A check was also made at this time at the 13.0- 
volt bias, and a value of 288P/T of 44.03 was 
found for an energy of 1.437 Mev, which agrees 
within experimental error with the main curve. 
From Fig. 3 it is seen that at the higher bias a 
cut-off pressure of 56.0 cm Hg corresponds to an 
alpha-particle energy of 1.79 Mev. This gives a 
Q value for the reaction of 3.94+0.08 Mev, the 
greatest part of the error resulting from the un- 
certainty in the cut-off pressure. 

The result is in excellent agreement with the 
other® determination of this Q value reported 
from this laboratory. As pointed out in the report 
of the previous work this value gives greater 


consistency with various series of reactions than 
the previously accepted value, and can be used 
with other well-established Q values and masses 
to give a He*® mass of 3.01688+0.00011 mass 
units. 

The present determination is independent of 
the slope of the range-energy curve and of the 
length of the absorption cell. Moreover, it does 
not depend on the Q value of the reaction 
Be*(p,a)Li® as does the previous determination. 
In a sense, then, the present work can be con- 
sidered as a check on the beryllium reaction. 


APPLICATION TO ALPHA-PARTICLE 
RANGE-ENERGY CURVE 


The data of Table I can be used with Eq. (1) 
to give information about the alpha-particle 
range energy curve in this region. If L and C 
were known the curve could be determined 
absolutely. If only LZ is known the data can be 
used to give the slope of the range-energy curve. 
In those experiments® where energies are deter- 
mined by comparison with particles of known 
energy, only the slope of the range-energy curve 
need be used. 

In determining the proper value of L to use 
in the case where the pressure in the ionization 
chamber is variable, one is confronted with such 
unknown factors as the point at which the 
collecting field of the ionization chamber be- 
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Fic. 4. Counts vs. pressure for alpha-particles from 
Li® (p, a) He*. The base line represents the He’ particles 
and the long range alphas from Li’, plus the background. 
The horizontal arrow indicates the half-maximum. The 
lines and numbers at the top indicate where the half- 
maximum would have to occur to give the various Q values. 
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Fic. 5. Determination of the effective length of the absorption cell. The figures between the curves are the computed 
values of the cell length using adjacent curves. The other figures are cut-off pressures (pressures at half-maximum). 
Locations of plateaus for the four curves at the lower pressures were determined by solid angle calculations. An average 
value of effective cell length of 0.740 cm is indicated, with some evidence for a variation in cell length of about 3 percent 


over region of pressures used. 


comes effective and the constancy of that point. 
In addition the variation in the recombination of 
ions at different pressures may be significant for 
such a device. There are also the usual uncer- 
tainties in determining the length of a short 
absorption cell such as the irregularities of the 
screens and the protruding of the collodion 
window through the screen apertures. Because of 
all these factors, the measured length (0.59 cm) 
of the absorption cell is only an approximation 
to the effective length. 

The following method® was used to get the 
effective length of the cell. An auxiliary chamber 
was attached in front of the absorption cell and 
this chamber was kept evacuated. It was thus 
possible to reduce the pressure in the absorption 
cell below atmospheric without rupturing the 
absorption cell window. A thin polonium source, 
mounted on an accurate metric micrometer, was 
placed in front of this chamber (which had a 
front circular opening of one-half inch diameter 
sealed by a collodion window of 13-mm stopping 
power) along the axis of the absorption cell- 
ionization chamber device. A number-pressure 
curve was taken for a given position of the source. 
With a cut-off pressure of P at temperature T 





the range of the particles is given by 
R=x(p/76)(288/T)+L(P/76)(288/T)+C, (2) 


where x is the length of path of the particles in 
air and p is the atmospheric pressure. For 
another position of the source, that is, for another 
value of x, we get another equation for R, since 
the range of the particles is the same. If C is 
assumed to be constant, the two equations can 
be solved for L. If » and T do not change 
between any two curves the value of L is inde- 
pendent of the absolute values of x but depends 
only on the difference of the x values, which 
can be accurately read on the micrometer. If p 
and T do change, the absolute values of x enter 
only in a small correction term. 

Preliminary sets of such curves gave an 
average value of L of 0.74 cm. Since the dis- 
crepancy with the measured value was large the 
cell length determinations were repeated with a 
collodion film on the screen between the absorp- 
tion cell and ionization chamber. This made it 
possible to keep the ionization-chamber pressure 
constant at atmospheric so that the device was 
being used in the conventional way. Two deter- 
minations were made under these conditions. 
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The values obtained were 0.61 cm and 0.60 cm, 
in good agreement with the measured length of 
0.59 cm. These results were taken as evidence 
that with no window between the two compart- 
ments, part of the ionization chamber was 
effectively part of the absorption cell. Though 
the entire explanation for this is not known, it 
seems Clear that the correct length to use for the 
electrostatic deflection experiments, in which the 
window was not present, is 0.74 cm. 

A final set of curves to determine the cell 
length without a film on the ionization-chamber 
opening was taken. In these experiments the 
alpha-particle source on the micrometer screw 
was moved only one-half mm between curves, in 
an attempt to find out if the effective cell length 
varied with pressure. A study of this variation 
requires the use of curves whose cut-off pressures 
are close together. A difference in position of the 
alpha-particle source of one-half mm gives a 
difference in cut-off pressures of about 5 cm Hg. 
However, if one uses two such curves to compute 
the effective cell length, an error in the difference 
of cut-off pressures of only 3 mm would result in 
a 6 percent error in the cell length, that is, an 
error of about 0.04 cm. Since the difference in 
cut-off pressures for any two positions of the 
micrometer source is reproducible to within 
about 3 mm, this method can determine the 
variation of cell length to an accuracy of about 
6 percent. 

The set of curves is plotted in Fig. 5. The 
locations of the plateaus for the four curves at 
the lower pressures could not be determined 
experimentally, as they would occur at pressures 
much below the critical pressure P. mentioned 
above, and the ionization chamber would not 
record all of the particles. These plateaus can, 
however, be calculated from solid angle con- 
siderations, since the number of counts on the 
plateaus for these experiments changes only 
because of the change in the solid angle which 
the opening to the absorption cell subtends at 
the source. The three plateaus that were deter- 
mined experimentally checked these calcula- 
tions. Though some of the points on these 
plateaus were taken at pressures slightly below 
P., the general appearance of the curves makes 
it certain that practically all of the counts were 
being recorded. The gradual decrease from the 
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plateaus which the curves of Fig. 5 exhibit is 
characteristic of number-distance curves for ordi- 
nary thin sources. The effect appears exaggerated 
in Fig. 5 because of the scale used for abscissae ; 
10 cm Hg correspond to only about 1 mm 
distance in air at standard conditions. 

The figure between any two curves indicates 
the cell length computed by using those curves. 
The average value of 0.74 cm for the cell length 
over the range of pressures used is substantiated. 
If the cell length is calculated by using the 
difference in cut-off pressures between the curve 
at the highest pressure and that at the lowest 
pressure, a 3-mm error in this pressure difference 
results in only a 1 percent error in cell length, 
that is, an error of less than 0.01 cm. This, then, 
is a measure of the accuracy with which the 
average effective cell length is known. 

Though the accuracy with which this method 
can determine the variation in cell length is not 
high, there does seem to be an indication in 
Fig. 5 of a downward trend in the effective cell 
length as the pressure is lowered. A variation in 
this direction of about 3 percent for the range 
of pressures used in the electrostatic analyzer 
data was assumed. The data in Table I was 
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Fic. 6. Points on the range-energy curve in the region 
studied. The data were matched to the Cornell curve of 
Holloway and Livingston at 1.437-Mev energy. 


used as follows. For pressures greater than 55 cm 
Hg L was set equal to 0.750 cm, for pressures 
between 35 and 55 cm a value of 0.740 cm was 
used for LZ, and for the two lowest pressures a 
value of 0.730 cm. An absolute range of 0.758 cm 
was assigned to the 1.437-Mev alphas, as deter- 
mined from the latest Cornell’ curve. This gives 
C a value of 0.323 cm. The results are plotted in 
Fig. 6. The points are seen to agree well with the 
curve of Holloway and Livingston. The spread 
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of the points is not serious, since the dispersion 
of the method is so high that the greatest spread 
corresponds to less than +0.01 cm in range. 
The chief uncertainty is in the value of L used. 
The results can be said to give an average value 
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by integrating their specific ionization curve) 
that the ratio between energy loss by ionization 
and by excitation is constant. Our result, then, 
is an indication of the validity of this assumption 
for this region. 
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The Velocity Spectrum of «-Particles 


Roy RINGo 
Ryerson Physical Laboratory, University of Chicago, Chicago Illinois 
(Received September 18, 1940) 


A 60° magnetic a-particle spectrograph has been constructed and used on a number of 
natural a-particle emitters, with photographic plates to detect the particles. The energy of the 
main group of protoactinium a-particles was found to be 5.053+0.007 Mev. With the thinnest 
sources of Po, ThC, and ThC’ a-particles, as many as 10 percent of the particles were found to 
have energies 5 percent or more below the energy of the maximum of any one group. 




































1. INTRODUCTION 


T is desirable to determine the energy set free 
in nuclear transformations with the greatest 
possible precision. This can frequently be done by 
finding the energies of the charged particles 
produced in the transformations. The most satis- 
factory method of measuring the energy of these 
particles is by their deflection in electric or 
magnetic fields. The difficulties involved in pro- 
ducing large electric fields make their use less 
suitable than magnetic fields for the deflection 
of the heavy particles of mass 1 or more and of 
energy greater than 2 Mev. Focusing magnetic 
spectrographs, for measuring the energies of 
the natural a-particles, have been built by 
Rosenblum! and Rutherford, Wynn-Williams, 
Bowden and Lewis,’ using the focusing properties 
of 180° deflection to increase intensities. Rosen- 
blum, working with the large electromagnet of 
the Académie des Sciences at Bellvue, detected 
the a-particles by the blackening they produced 
on photographic plates. His spectrograph had 
~ 1, Rosenblum, J. de hys. et rad. [7] 1, 438 (1930). 


2 Rutherford, Wynn-Williams, Bowden and Lewis, Proc. 
Roy. Soc. A139, 617 (1933). 


some disadvantages. The magnet could not be 
moved, hence it could not be used on reactions 
produced by artificially accelerated high energy 
particles. The method also required very strong 
sources. 

Rutherford and his group used an annular 
electromagnet and detected the particles in an 
ionization chamber connected to an amplifier and 
counter. Rays of different velocities were brought 
into the chamber by varying the magnetic field. 
The type of magnet used was limited to 180° 
deflections, this required the source and detector 
to be in the magnetic field. In addition the 
sensitivity of the detector was limited by the 
background counting rate produced by contami- 
nation and cosmic rays. These two limitations 
restricted the usefulness of the instrument mainly 
to a-particles from strong natural radioactive 
sources. The energies of all the a-particles from 
the natural radioactive elements with disintegra- 
tion rates equal to or greater than that of radium, 
have been measured by these two instruments. 

In this investigation a spectrograph was used 
which deflected the particles through 60° and 
which has the source and detector about 40 cm 
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from the magnetic field. This follows a suggestion 
of Herzog* that the energies of a-particles could 
be measured conveniently by the use of such 
deflections. The a-particles are detected in this 
apparatus by means of a photographic plate in 
which the tracks made by the individual particles 
are found by microscopic examination. This 
method records every particle and has practically 
no uncertainty due to background. There is also 
the advantage that exposures of indefinite length 
can be taken, permitting the use of very weak 
sources. 


‘2. DESCRIPTION OF THE APPARATUS 


The two core pieces of the magnet, on which 
the coils are wound, are 37 sq. in. in cross- 
sectional area and 9 in. high. The yoke is a hollow 
iron rectangle made of four pieces, 24 in. X10 in. 
in cross section, providing a return path for the 
flux on both sides of the core and windings. The 
core pieces are mounted inside this rectangle one 
attached to the top piece and the other to the 
bottom. The pole pieces and the assembly which 
holds them in alignment, fit between the core 
pieces and are easily inserted and removed. The 
iron parts are all made from Armco Ingot 
forgings, chosen for their high saturation flux 
density. 

The two pole pieces are trapezoids at the back 
surfaces (8$ in. and 124$ in.X3 in.) which fit 
against the core, and taper down to curving pole 
faces 1;% in. wide, as indicated in Fig. 1. They 
are 2} in. high and tapered at an inclination of 
45° to the pole face. Four columns set in holes in 
the pole pieces as shown in Fig. 1 keep them in 
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Fic. 1. A horizontal cross section through the spectro- 
graph, taken at the level of the middle of the gap between 
the i. The vacuum line and the magnet windings 
(on the core) are omitted. 
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+R. Herzog, Zeits. f. Physik"89, 447 (1934). 
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alignment. The columns are steel below the 
surface of the pole pieces and brass above. They 
were pressed into their holes on the sloping sides 
of these pieces and the whole face side of each 
piece, including the tops of the columns, was 
ground off to the same level. The two pole 
pieces are clamped together by means of brass 
rods, threaded at both ends, running through the 
columns and held by steel nuts in holes counter- 
sunk in the back surfaces of these pole pieces. 





Fic. 2. The path of a charged particle through a mag- 
netic field with the boundaries shown. The direction of the 
field is perpendicular to the plane of the paper. 


The pole pieces are separated from each other by 
brass washers, } in. thick, between their columns. 
This insures the pole faces being parallel and 
separated by a constant distance. 

The windings consist of 1000 turns of No. 12 
wire on each core piece, cooled by copper sheets 
between every other layer of the windings. This 
arrangement carries a current of 2 amp. with 
little heating after an hour’s operation. This 
current was sufficient for the a-particles studied. 
A current as large as 8 amp. can be safely carried 
if a desk fan is turned on the coils. 

The whole instrument weighs about 1200 
pounds and is mounted on a truck so that it may 
be moved easily. 

The magnet current was drawn from a 48-volt 
battery of lead cells. It was measured by a 
potentiometer and held constant, to within 1 
part in 3000, by varying a rheostat in series with 
the magnet. 

The source and plateholder were designed 
according to the equations for the focusing of ion 
beams given by Herzog.’ The relation between 
the position of the source and the position of the 
focus for an ion beam with the central ray 
incident normally on, and emerging normally 
from the boundaries of a homogeneous magnetic, 
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Fic. 3. The number of a-particles from the thicker Po 
source reaching 0.64 mm? of the plate, plotted against the 
distance from the low energy end of the plate. In the region 
of the maximum 1 cm change in position represents 1.5 
percent change in energy. The exposure time was 1 hour. 


field is 
(Y—g’)("—g") =f", (1) 
where I’ is the distance of the source, or object, 
from the edge of the field, /’’ is the distance of the 
focus, or image, from. the other edge, g’ is the 
distance from the edge of the field to the focal 
point on the source side (i.e., the point where a 
parallel bundle of rays, incident normally to the 
edge of the field on the image side, should 
focus), and g” is the corresponding distance on 
the image side. In this case g’ = g’’ = 14.03 cm and 
f is the focal length, the distance from the principal 
plane to the focal point, which for 60° deflection 
is 2g’. For the case l’=/’’, (1) gives /’=42.09 cm. 
For a particle going through the field in a path 
whose radius of curvature p differs from that of 
the particles following the middle of the pole 
faces, the angular deflection and the focal position 
are different. For this more general case we have 
from Herzog’s paper 
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Equation (1) holds also. The symbols have the 
meanings defined above or shown in Fig. 2. 

In the instrument under discussion, with e’ =0, 
the maximum radius of curvature of any path 
which will pass through the diaphragm at the 
exit from the field, when it is open 1 in., is 25.8 em 
and the minimum is 22.9 cm. From this, and 
Eqs. (1)—(4), we have /’’=36.0 cm and 51.9 em 
in the two cases. By finding the direction of the 
ray by geometrical construction on a drawing of 
the pole face, and laying off /” along it, one 
obtains the position of the focus for these 
extreme rays. The plateholder is designed so that 
the plate passes through the positions of the foci 
of the rays of maximum, middle and minimum 
radii. 

The source mounting is 42.1 cm from the 
entrance to the magnetic field so that I’ will 
equal /” for the rays entering and leaving the 
field normally. The plateholder is in the position 
indicated in Fig. 1. It is set at an angle of 21° 
from the normal to the edge of the field. At this 
angle of incidence the a-particles leave tracks, 
each one of which is observed as a row of de- 
veloped grains on the developed plate but which 
has sufficient depth to be distinguished from a 
surface scratch. Because a vacuum of 10-* mm of 
mercury is sufficient, it is convenient to break the 
vacuum every time the plate is changed and no 
special equipment is used for inserting the plates. 
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Fic. 4. The number of a-particles from a ThC source 
on 1 mm? of the plate, plotted against the distance from 
the low energy end of the plate. 
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A vacuum good enough to remove the danger 
of change of charge of the particles during passage 
through the magnetic field is all that is needed, 
for scattering is unimportant with such energetic 
particles. At a pressure of 10-* mm the mean free 
path for the gain of an electron is 10° cm, and for 
loss of an electron by He? particles, which are 
present in the radiation from even thin natural 
radioactive sources to the extent of one percent, 
the mean free path is 10‘ cm.‘ Either process thus 
is quite unlikely in the 26 cm of path in the field. 

The photographic detection of the particles 
follows a technique recommended by Wilkins.® 
The plates used were Eastman a-particle plates 
1 in. X8 in., which were developed 4 min. 30 sec. 
at 18°C in Eastman x-ray developer. The plates 
were then cut into three pieces after the separa- 
tion of fiducial marks on each piece was measured, 
and examined in a microscope using 450 X mag- 
nification and bright field illumination. The 
distribution of a-particles over the plate is found 
by counting the number of a-particle tracks in a 
definite vertical range as the plate is moved along 
horizontally. 


3. PRELIMINARY OBSERVATIONS 


The behavior of the spectrograph was first 
studied with a-particles from polonium and 
thorium C. The distribution of a-particle tracks 
after an exposure of one hour to a polonium 
source, is shown in Fig. 3. The source was a silver 
wire, 0.3 mm in diameter, which had been dipped 
in a radium D solution in which polonium was 
growing. This source emitted 6.7 X 10° a-particles 
per sec. in all directions. The graph shows the 
horizontal position on the plate as abscissa and 
the number of tracks observed in a strip 2 mm 
high and 0.32 mm wide, at that position as 
ordinate. Figure 4 shows the number on a strip 3 
mm high from a thorium C source emitting 
approximately 7X10‘ particles per sec. in all 
directions. The exposure time was two hours. The 
ThC source was prepared by collecting the active 
deposit of thoron produced by some radio- 
thorium, on a charged platinum wire 0.2 mm in 
diameter. 

In Fig. 4, the two lines photographed by 
Rosenblum! and referred to as ap and aj, are seen 


*E. Rutherford, Phil. Mag. 47, 277 (1924). 


‘T. R. Wilkins, J. App. Phys. 11, 35 (1940). 
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to be separated. The two very faint lines recorded 
in his photographs, a2 and a3, do not show up, 
probably because of insufficient exposure. Many 
tracks of particles with velocities less than those 
of the two main groups appear. Similar tracks of 
low velocity particles from Po are seen in Fig. 3. 
These would not have been expected from the 
clear backgrounds shown in Rosenblum’s photo- 
graphs and their origin will be discussed in 
Section 6. 
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Fic. 5. The ratio of the radius of curvature in the mag- 
netic field, corresponding to a point s on the plate, to 
the radius of curvature of the particles arriving at A, 
plotted against the distance from the low energy end of 
the plate, s. This ratio is equal to the ratio of the momenta 
of the particles arriving at s to those arriving at A. A is 
arbitrarily chosen as the point at which the Po @-particles 
happened to fall on a plate used in drawing the curve. 


4. CALIBRATION OF THE SPECTROGRAPH 


The strength of the magnetic field was meas- 
ured by a method of Ellis and Skinner,* which 
was modified slightly. The charge produced by 
removing a small search coil from the field was 
balanced against that produced in a search coil 
in a long solenoid when the current in the 
solenoid is interrupted. A ballistic galvanometer 
was used for the balance. The results obtained by 
this method were checked by measurements of 
the Hp of a-particles of known energy. The 
results from these two methods of the measure- 


*C. D. Ellis and H. W. B. Skinner, Proc. Roy. Soc. 


A105, 60 (1924). 











ments of H agree within the probable error, which 
was about two percent, because of stray fields, 
lack of precision in measuring the angle of 
deflection and variation in the magnetic field 
corresponding to one setting of the magnet 
current. It was found that a current of 1.5 amp. 
produced a field of 11,900 oersteds between the 
pole faces, 3 amp. a field of 21,600 oersteds, 
and 6 amp. one of 25,000 oersteds, at which 
current the iron was nearly saturated. With a 
field of 25,000 oersteds, a-particles of enérgy as 
large as 19.5 Mev could be brought on the plate. 

The dispersion, defined as the separation of 
particles differing one percent in momentum, was 
obtained by measuring the distribution of par- 
ticles from several natural radioactive elements, 
Po, ThC, ThC’, and RaC’. A wire coated with Po 
and one coated with the thorium-active deposit 
were mounted together in the source position 
and the separation of the a-particle groups 
measured. The ThC a-particles (E=6.086 Mev, 
Hp=3.532 10°) were 103.5 mm from the Po 
particles (E=5.303 Mev,? Hp=3.302X105). 
Hence the average dispersion is 15.3 mm for a 
one-percent change in momentum. 

To each point on the plate, distant S from one 

end, corresponds a definite radius of curvature 
(p) in the magnetic field. Thus if we plot a curve 
of the position, S, of tracks of several groups of 
a-particles of known energy, for some one value, 
H, of the field, against Kp for these particles, 
where K is an arbitrary constant, we may read 
off the curve, for any two points on the plate, the 
ratio of the p’s corresponding to the a-particles 
that caused the tracks appearing at those two 
points. Since the momentum of a particle is 
proportional to p, this ratio is also the ratio of the 
momentum of the two a-particle groups. Such a 
curve is given in Fig. 5. The ordinate is the ratio 
p/pa=Kp and the abscissa, the distance S along 
the plate, A is arbitrarily chosen as the point at 
which the Po a-particles happened to fall on one 
exposure. 

Two points on the curve in Fig. 5 are given by 
the results with ThC and Po mentioned above, 
and the slopes at various points, which are 
inversely proportional to the dispersions at these 
points, are given by the separation of the two 
main groups of the ThC a-particles. The slope is 
found to be 1.5 times as large at the position of 


946 ROY RINGO 



























































180 
@- PARTICLES hes Po 
150 
120 
4 
r 
wn 90 
x 
Vy 
< 
e« 
~ 
s 60 
cs 
z 
30 
° 
ry oO 
0 ~~ : a 
65 60 95 110 125 140 155 


MM. FROM END OF PLATE 


Fic. 6. The number of a-particle tracks from the thinner 
Po source on 0.64 mm? of the plate, plotted against the 
distance from the low energy end of the plate. Exposure 
time was 4 hours. | 


the Po a-particle tracks as at the position of the 
most intense ThC a-particle group on the plate 
on which both are recorded. From these data, the 
calibration curve which is not far from a straight 
line, can be drawn as shown in Fig. 5. 

In Herzog’s theory,’ the dispersion, D, is given 


by 
f 


p= 1+ 

100 l'—g’ 
(= 2p/100 in our case). However this is dispersion 
in a direction perpendicular to the beam; to get 
the dispersion along the plate, it is necessary to 
divide by the sine of the angle between the plate 
and the beam, 21°. This gives D = 13.6 mm for one 
percent change in momentum at the position on 
the plate where the beam leaving the field 
normally strikes, that is 8 cm from the left-hand, 
or low energy end of the plate. This is very near 
the measured value of 13.7 mm. 

In using the calibration curve it is assumed that 
the dispersion is independent of the absolute 
value of the field. This is true if the stray field 
does not change too much in proportion to the 
field proper as the latter is changed. In this 
spectrograph with a gap of } in. between the 
pole faces and operating well below the point of 
saturation of the iron, the change in dispersion 
due to stray field is negligible. This was checked 
by the results obtained with ThC’ and RaC’ 
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Fic. 7. The number of a-particle tracks from a ThC’ 
source on 0.64 mm? of the plate plotted against the dis- 
tance from the low energy end of the plate. Exposure time 
was 4 hours. 


a-particle sources together in the spectrograph. 
The ThC’ a-particles (E=8.778 Mev?) reached 
the plate 15.05 cm from the low energy end and 
the RaC’ a-particles (E=7.683 Mev’) 4.96 cm 
from that end. These values fit the calibration 
curve. : 

The resolution obtained so far, is at best 
barely equal to that theoretically expected with 
the width of sources used. The width of the lines 
is in theory’ equal to pf? sin }@ where 8 is the 
angle between the normal from the source to the 
boundary of the field and the path deviating the 
most from that, which still passes through the 
diaphragms in the field. The width comes out to 
be 0.2 mm. This width could not be attained 
because of the width of the sources, which, being 
over 0.3 mm, gives images over 1 mm wide. 
Some lines were obtained whose width at } of 
their maximum was 1 mm. These give a resolu- 
tion (equal to the dispersion X 100/line width) of 
1300, but most lines were much broader, probably 
because of absorption of energy in the source for 
some particles, and unsteadiness of the magnetic 
field, which, on the long exposures taken with the 
thinnest sources, varied by more than 1 part in 
1000. For the immediate purpose of calibration 
and for many measurements this broadening of 
the lines is not a serious disadvantage however. 


™W. E. Stephens, Phys. Rev. 45, 513 (1934). 
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The fraction of the total of emitted a-particles 
which gets through this spectrograph at the 
maximum usable field width of 2 cm, is 1 out of 
70,000. This compares with 1 out of 20,000 
detected in Rutherford’s magnet? with a dia- 
phragm 3 cm X1 cm. The theoretical and attained 
resolutions of Rutherford’s magnet were slightly 
larger than those of this magnet. 


5. ENERGY OF THE PROTOACTINIUM a-PARTICLES 


Because of the method of detection employed, 
this instrument can be usedl to measure the 
a-particles from quite weak sources. One such 
source whose a-particles were measured was a 
preparation of protoactinium, kindly lent by Dr. 
A. V. Grosse. It contained seven percent PaO; 
mixed with zirconium and hafnium oxides. A 
source which emitted 1200 particles per sec. was 
prepared by depositing this material on a 
platinum strip 1 mm X5 mm in area. The value 
for the high energy limit of the a-particles from 
Pa obtained by comparison with the a-particles 
from Po, put in the spectrograph at the same 
time, was 5.053+0.007 Mev. This is to be com- 
pared with a value 5.063 Mev, found by using the 
Pa a-particle range, 3.673 cm, and the Po a- 
particle range, 3.925 cm, in air at 15 C and 760 
mm pressure, given by Geiger,® the slope of the 
Cornell University range-energy curve® in this 
region, and the value 5.303 Mev for the energy of 
the Po a-particles? used with the data from the 
spectrograph. 

The thickness of material necessary to obtain a 
source of the strength used, caused such a 
broadening of the line by retardation of the 
particles in the source that nothing could be told 
of any possible fine structure. This illustrates an 
inherent limitation of the apparatus when used 
on weak or impure radioactive materials. 


6. Low ENERGY a-PARTICLES 


In addition to the widening of the peak of the 
a-particle lines observed, and ascribed in part to 
loss of energy in escaping from the source, many 
a-particles of much lower energies are observed 
(cf. Fig. 4). These have been observed before by 


8H. Geiger, Zeits. f. Physik 8, 45 (1921). 
*M. G. Holloway and M. S. Livingston, Phys. Rev. 54, 
18 (1938). 








I. Curie!® and J. Schintlmeister,"' and are 
noticeable on curves published by Rutherford.’ 
The results in earlier work were all obtained with 
sources 100 or more times as thick as the com- 
parable sources used in this work and were 
generally ascribed to absorption in the sources. 
The results obtained with this spectrograph are 
shown in Figs. 3, 6 and 7. The density of deposit 
on the source producing the spectrum shown in 
Fig. 3, was 9.93X10-* g per cm?, equivalent in 
stopping power to 5.310-‘ mm of air. Figure 6 
shows the distribution of a-particles from a Po 
a-particle source of density 5.29X10-* g per 
cm? or 2.8X10- mm air equivalent, deposited on 
a flat strip of silver 1 mmX5 mm X0.1 mm. 
Figure 7 shows the distribution of particles from 
a ThC’ a-particle source of 0.5X10-" g per cm?, 
deposited on a similar platinum strip. 

The Po sources were prepared by precipitation 
on silver from a RaD solution whose only 
contamination was believed to be inactive lead. 
The ThC’ was prepared as was the ThC source 
described previously (Section 3). The sources all 
appeared clean when first prepared and used but 
darkened later, presumably because of oxidation 
resulting from ionization by radiations from the 
radioactive materials. 

In the case of the thinner Po a-particle source 
2.6 percent of the particles were retarded by an 
amount equivalent to 200 to 300 kev, for the 
other Po a-particle source, 4.3 percent, and for 
the ThC’ a-particles, 4.0 percent. To retard the 
particles by an amount equivalent to 250 kev a 
path in air approximately 1.8 mm long is required. 
This means a layer of air 510° atoms deep or 
ys aS many atoms of the radio-elements if the 
ordinary free paths for the a-particles are as- 
sumed. Assuming a very favorable case for the 
absorption in the radioactive material, that it is 
all concentrated in a small area and forms a cone 
with base diameter equal to the altitude, we find 
that such a cone, high enough to account for the 
retardation observed, would contain 3 times as 
much material as there is Po on the stronger Po 
source, 50 times the amount on the other Po 
source and 5000 times the amount of thorium- 

active deposit present on the ThC’ a-particle 
source. 


101, Curie, Ann. de physique 3, 299 (1925). 
"J. Schintlmeister, Sitz. d. Akad. Wien, Abt. 2a 146, 
389 (1937), 
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The possibility that the active material might 
have deposited in clumps, was tested as follows. 
For a short time, the sources were placed in 
contact with a photographic plate which was 
protected from infection by the radioactive ma- 
terial by a thin sheet of collodion. The plate was 
then developed and examined under a micro- 
scope. The a-particle tracks were found to be 
quite uniformly distributed in the case of the Po; 
no centers were observed. In the case of the ThC 
many distinct clusters were seen while more than 
? of the area touched showed very few tracks, 
However the total amount present here is go 
small that these clusters, by the argument of the 
preceding paragraph, cannot account for the 
retardation observed. 

Because the retardation does not depend much 
on the average density of deposit, an increase of 
20 times in the density of Po increased the 
fraction retarded only 60 percent, and because 
there are hardly enough atoms present in the 
most favorable case to account for the retardation 
observed, it is believed that absorption in the 
source cannot be held responsible for these low 
energy particles. 

Scattering, on the pole faces, or in air, is 
excluded as a source of these particles, as it 
should also cause an appreciable number to 
appear on the high energy side of the peak, by the 
alteration in direction it causes, and this region is 
remarkably free of tracks. However, as a pre- 
caution against scattering the pole faces had been 
coated with Aquadag, and diaphragms inserted, 
which collimated the beam vertically as well as 
horizontally, so that the number of particles 
reaching the pole faces was reduced. 

If absorption in the source is also excluded as 
an explanation of these low energy particles, 
these experiments then show that there is present 
among the a-particles of homogeneous velocity, a 
considerable number, leaving the radioactive 
atoms, or at most small groups of such atoms, 
with a continuous distribution of energies going 
down as low as ten percent below the energy of 
the high energy limit, or lower. This number may 
include ten percent of the total number of par- 
ticles in energy ranges five percent below the 
maximum. 

In conclusion the author wishes to express his 
gratitude for the advice and guidance of Dr. A. J. 
Dempster, who suggested this problem. 
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Families of odd nuclei differing from each other only by a-units show some evidences of a 
certain type of regularity in their stability against 8-processes. The assumption that a similar 
regularity holds for all such families permits the unique assignment to the correct isotope of 
many known radioactive periods and the prediction of many vet unknown lifetimes. 


Y assuming the existence of certain regu- 
larities among the half-lives of the radio- 
active nuclei as functions of their atomic num- 
bers, the assignment of discovered half-lives to 
the proper isotopes may be helped and the 
periods of yet unknown isotopes may be esti- 
mated. The existence of such regularities might 
be expected because of the well-known regu- 
larities among the energies of nuclei, if a unique 
correspondence between life and energy release 
could be assumed. Such a uniqueness is, of course, 
actually destroyed by the possibility of ‘‘allowed”’ 
and “forbidden”’ transitions. However, one can 
hope that the nuclei chosen for comparison with 
each other (henceforth referred to as a ‘‘family’’) 
will behave so similarly toward such disturbing 
factors as to result in a dependence on atomic 
number smooth enough to be useful. The ap- 
proach is almost entirely empirical; the known 
half-lives of a family are compared and whatever 
regularity is found is extrapolated to predictions 
about unknown and uncertain members. 

The criterion we have adopted for the forma- 
tion of a family is that the members differ only in 
the number of a-units (of two protons and two 
neutrons each) they contain. We have chosen this 
in preference to other possible criteria! because 
of the fact that the occurrence and even the 
abundance of the stable isotopes shows great 
similarity for alternating values of the atomic 
number, Z, throughout the periodic table and 
makes promising the extensibility of the scheme 
even to the heaviest nuclei. It also has as a result 
that the stable isotopes of each family follow each 


1Our scheme is consistent with Professor E. Wigner’s 
treatment (Phys. Rev. 51, 106 (1937)) of nuclei differing 
by four particles as having the same multiplet structure 
in a first approximation. This may help in maintaining 
a similarity of behavior with regard to “forbidden” 
transitions. 
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other consecutively as Z increases.? The finite 
lifetimes of the family should then be expected to 
fall off smoothly as Z increases or decreases from 
its stable values. 

Our specific procedure was to plot the loga- 
rithms of the known half-lives against their 
atomic numbers. Results are shown in the figure 
for several families for which the data are most 
abundant. The families are distinguished by such 
symbols as (3e) and (50), in which the number 
represents the neutron excess (isotopic number) 
while the letters ‘‘e’’ and ‘‘o’’ denote families 
with even and odd Z, respectively. The plots 
show some evidence of the expected behavior, 
especially among the positron emitters (nuclei 
with a larger than stable Z) of the families 
(—1e), (—10), (7e) and among the negatron 
emitters of (3e). All these are families of odd 
nuclei. All the families of even nuclei without 
exception display extreme irregularity, even at 
the beginning of the periodic table as is demon- 
strated by the family (00) in the diagram. We 
must, therefore, at the outset confine our.method 
to the odd nuclei. The completely random be- 
havior of the even nuclei serves only to make the 
few definite evidences of regularity among the 
odd isotopes appear less accidental. 

The diagram shows not only the families of odd 
nuclei giving the best agreement with our 
expectations but also those containing the most 
numerous exceptions to the expected behavior. In 
the case of every such exception, some question 
concerning the experimental evidence can be 
raised, as is shown in the following paragraphs. 


? This also is true to a large extent of nuclei differing only 
by one neutron and one proton, but exceptions are fre- 
quent (e.g. J=1, 3, etc.,). Moreover, whatever advantage 
the argument of footnote 1 gives would be lost under this 
criterion. 
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Fic. 1. Logarithms of half-lives, 7, as functions of atomic number Z for several families, as denoted by symbols such 






as (—1le) (ie. J=—1, Z even. To get mass number, double Z and add J). Most of the data are from Livingood 
and Seaborg’s summary (see footnote to Table I). The upward arrows (f) denote stable isotopes. The full circles 


(@) represent the best established radioactive isotopes, of classes 


A and B according to Livingood and Seaborg, 


The crosses (X) denote isotopes not previously definitely yer The assignments best agreeing with our scheme 


were chosen for these cases (see Table I). The small dots (- 
55, 691 (1940)) with the assumption that the transition is ‘‘allow 
the datum has been put in doubt from the standpoint of our scheme as discussed in the text. T 


represent estimates from Barkas’ masses (Phys. Rev. 
” A diagonal dash through a pee signifies that 
e squares are new 


assignments of some of the newly doubted isotopes. The dashed curves were used in making the pane of 


Table II. Cr®, of the family (7e), was rang predicted to have a 2- to 3-hour period on the basis o 
aboratory by McDaniel and Dickson (Phys. Rev. 57, 351 (1940)) and 


the ye my verification was done at our 
a'so by others (Amaki, Iimori and Sugimoto, Phys. Rev. 





The final decisions must, of course, be left to 
experimental test. 


(le).—The activity assigned to Ni*’ fits better in the 
family (3e) as Ni®. Only the failure to observe it in neutron 
capture by Ni®* seems to have prevented the latter assign- 
ment; a low capture cross section may be responsible 


instead. 
(lo) and (30).—Co*’ is one of only three exceptions to 


the rule that stable members of a family follow each other 
consecutively when ordered according to Z.* This makes 


3 A second exception is the rare isotope, Rh™, isobaric 
with Ru’, Only the finding that Ma*’ is stable would 
help this case. The third exception among the odd nuclei 
is due to the absence of Ce" from the stable nuclei of 
family (25e). This rare earth isotope may yet turn up when 
the mass-spectrographic investigation of Ce is pushed to 
the accuracy usual in other parts of the periodic table. 
A similar exception exists in family (130) due to the near- 
stable isotope Rb*’. Only the rather farfetched assumption 
that this long period is isomeric to a much faster decay 


this scheme; 


57, 752 (1940)). 





it desirable to verify further its mass-spectrographic 
observation,‘ especially since its isobar Fe’ is stable. 
Moreover, the supporting evidence offered by Sampson, 
Ridenour and Bleakney,’ that the two negatron periods 
which arise in Co through neutron capture must be due to 
two stable isotopes, seems now to be nullified. Co** produced 
in several other ways emits positrons with a totally dif- 
ferent period. The existence of Co*’ as a target of trans- 
mutations has, therefore, not yet been proved. 

Only the supposition that Co*’ is stable has prevented 
the assignment to it of the activity presently given to 
Co; the transfer would help the agreement in both the 


from the ground state seems able to make it consistent 
with our scheme. There may be some support for this in 
the fact that Sr*’, which differs from Rb* only in that a 
proton replaces a neutron, has isomeric states. 

a 936) B. Sampson and W. Bleakney, Phys. Rev. 50, 732 
a . _ Ridenour, and Bleakney, Phys. Rev. 50, 382 
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TABLE I. New or confirmed assignments of known periods. 

















ASSIGNMENT PERIOD FORMERLY 
K* 18 min. K‘43 44 
Ge*® 6-10 days Ge*®®. 71 
Ge™ 195 days Ge®?, 6% 7 
Ge™* 37s. Ge™ 
As” 88 min. As". 73 
Se*! 19 min. Se7® 8 
Kr”? ~1 day Kr? 8 
Rb# 42 min. Rb? 8. 83, (87) 
Rb* 200 hr. Rb?® 81, 83, (87) 
Sr™* 70° =min. Sr® 
Sr”! 6 hr. Class D 
Zr”? 6 min. Class F 
Mo" 17 min. Mo. 93 
Ru” 90 min. Class E 
Agi” 40 sec. Agi07, 109 
Sn! 40 min. Class D 
Xel33 9.4 hr. Class D 
Eris 12 hr. Erisain 
Er!” 7 min. Eri69,171 
Os! 40 hr. Os!91.193 
Hg? 25 hr. Hg? 206 








1See J. J. Livingood and G. T. Seaborg, Rev. Mod. Phys. 12, 30 
(1940). 


(30) and (10) families. The agreement would become com- 
plete if, further, the activity which has been tentatively 
assigned to V*’ were to be attributed to the ground state 
decay of V** instead. The possibility of this isomerism was 
recognized by its observers.® 

(3e) and (7e).—On the basis of our scheme, the ground- 
state decay of Cr®*! is expected to be too slow to be easily 
observed (a period of some decades). That the far shorter 
period actually observed may be from an isomeric state is 
supported by the appearance of a highly converted half- 
million-volt gamma-ray with it. A similar explanation 
seems also necessary for the too short Ni® period in family 
(7e). 


6 Walke, Williams and Evans, Proc. Roy. Soc. A171, 360 
(1939). Note added in proof—In agreement with our con- 
clusion is L. A 


A. Turner, Phys. Rev. 58, 679 (1940). 
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For most of the families not shown or specifi- 
cally discussed, the data are not sufficient to make 
a decision about the validity of our scheme. In 
four of these families (19e, 250, 25e, 37e) there 
occur pairs of well established nuclei which in no 
case conflict with our scheme. 

In Table I are listed assignments of known, but 
previously unassigned, activities made to agree 


TABLE II. Approximate half-lives predicted for 
yet unknown isotopes. 











IsoTOoPE PERIOD IsoTOoPE PERIOD IsOTOPE PERIOD 
BY 0.01 sec. Ca*® 1 sec. Kr77 1 hr. 
cs 20 sec. Ca‘? 10 sec. Kr® 10 yr. 
N" 0.1 sec. Ti*® 5 hr. 43% © 

F2 0.5 sec. Cr* 1 hr. Ru! 1 min. 


Na*® 20 sec. Cr 20-200 yr.| Ru’? 10 sec. 
Mg** 0.5 sec. Mn*® 1 min. Pdi 1 hr. 


Al* 10 sec. Co” 20 hr. Pd"3 1-60 min. 
Al* 1 sec. Ni? =. 2 min. Cd" 3 min. 
Si® 5 sec. Ni®® 30 hr. In13 30 min. 
p29 5 sec. Ni® 20 yr. Sn? 10-50 min. 
pss 10 hr. Cu®® 20 min. Te'*® <20 sec. 
p% 5 sec. Zn@ = 1 min. Csi 3 yr. 
S$ 2 sec. Ge? 2 min. Ba"? < 2 min. 
S37 10 sec.-| Ge7® 20 min. Celt ce) 

cl 1 min.| As® 2 min. Sm“ > 1 yr. 
A® 1 sec. Se® 1 min. Gd! 40 min. 
K37 1 sec. Br® ~=5 min. Os! 2 min. 
K*® <20 min. 














with our assumptions. Table II contains esti- 
mates of yet unknown half-lives as obtained from 
extrapolations and interpolations along curves 
such as those shown in Fig. 1. Both tables contain 
values having varying degrees of certainty, de- 
pending on the amount of data available con- 
cerning the corresponding families. More assign- 
ments and predictions will be possible as more 
information becomes available. 
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HE theory of charged mesons of spin 1 leads 
in its present form to divergent or otherwise 
unreasonable results not only in the majority 
of nuclear problems, but also in the treatment of 
the interaction of mesons with the electromag- 
netic field. The simplest instance is the scattering 
of mesons by a point charge. Massey and Corben! 
and Oppenheimer, Snyder and Serber” have shown 
that in the Born approximation the cross section 
for this scattering does not decrease with the 
energy E of the meson, if E>wc*, but tends to a 
finite limit. Still one could hope that the exact 
solution of the problem would lead to a reason- 
able result. However, this is not the case. 

The meson wave function has 6 independent 
components which correspond to the components 
of the vector potential and of the electric in- 
tensity of the electromagnetic field. Let us in- 
vestigate the stationary solutions (V~e“*) of the 
Proca equations in a central field of force for 
h|w| >c* (unbound states). The complete set of 
these solutions consists of three independent 
eigenfunctions for any given set of values of w, 
j and m, such that |m|=j and j=1 (for j=0 
there is only one eigenfunction). Here jh denotes 
the total angular momentum of the meson and 
mh its projection on the Z axis. The three eigen- 
functions correspond to three possible orienta- 
tions of the spin s of the meson with respect to 
its orbital momentum /. However s and / are not 
constants of the motion: One of the eigenfunctions 
is of the form 


Wy, =f,(r)- V;7a(8, g)-e*', A=1, 2---6, (1) 


where m,=m or m.=m+1, while two others are 
of the form 


vra= tft (r)- Vi41"(8, ¢) 
+fx-(r) Vi-™™(0, g) fe, (2) 


where f,+ and f,~ are not independent of one 


another. 


1H. F. W. Massey and H. C. Corben, Proc. Camb. 


Phil. Soc. 35, 463 (1939). 
2 J. R. Oppenheimer, H. Snyder and R. Serber, Phys. 


Rev. 57, 75 (1940). 


If the potential energy V(r) of the meson has 
no singularity at r=0 and if 721, then for fixed 
w, j and m there is one regular eigenfunction of 
the type (1) and two of the type (2). However, in 
the field of a point charge e’V=ee’/r there js 
again one regular eigenfunction of the type (1), 
but the solutions of the equations for f,+ and f,- 
are of the form: 


fypt~re 4s. ers, 


CGH ee! 


pc? r 


where 


s and 7=1', 

i.e., n= +1 or n»=-+7. Not two, but only one of 
these four solutions is finite at r=0; the three 
others lead to an infinite value of the energy 
E= fWHy-dr of the meson.* 

Thus the regular solutions of the Proca equa- 
tions in a field of a point charge do not form a 
complete set of functions and the problem of 
the Coulomb scattering of mesons has no solution. 

The reason for this breakdown of the theory 
may be sought in the neglect of the finite size 
of elementary particles. If that is the case, the 
Coulomb scattering of mesons with energies 
E2hc/ro must substantially depend on the size 
ro of elementary particles. This will also be true 
for the bremstrahlung and probably also for the 
production by mesons of fast 8-electrons. Thus 
the estimates of the cross sections for these 
processes based on the results of the Born 
approximation cannot be valid for E 2 he/ro. 

We hope to discuss in a further note some 
further consequences of the peculiar behavior of 
mesons in a Coulomb field (in particular the 
bound states of the mesons). 

The present note is a summary of a part of a 
paper which will appear in the Journal of Physics, 
edited by the Academy of Sciences of the USSR. 

3 If one calculates the charge density p of mesons with 
the help of the least objectionable nonregular solutions 
(n= +7) one finds that in the neighborhood of the point 
charge e’p~r~*'? and that the sign of p is uniquely deter- 
mined by the sign of e’ and is opposite to it. This means, 
that if e.g. e’>0 an infalling negative meson may 
caught by e’ in a spiral orbit, whereas an infalling positive 


meson produces pairs, the negative components of which 
are caught by e’. 
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The Electromagnetic Properties of Mesotrons 
H. C. CorRBEN* AND JULIAN SCHWINGERT 
Department of Physics, University of California, Berkeley, California 
(Received June 27, 1940) 
AS A general theory describing particles of unit spin and arbitrary magnetic moment is developed 
d and applied to the motion of such particles in a Coulomb field. In the particular case of mag- 
of netic moment unity (Proca theory), the exact equations for the radial components of the wave 
' functions possess regular solutions only for those states (j=/#0 and j=0, /=1) in which the 
in ; ; ; 
: orbital angular momentum / is a constant of the motion. For particles possessing a magnetic 
1s moment of two mesotron magnetons, the radial equations are free of singularities for all states 
), but two: j=1, /=0 and j =0, /=1. The cross section for a fractional energy transfer to electrons 
- by energetic mesotrons is calculated for the various simple possibilities of mesotron spin ¢ 
(0, 3, 1) and magnetic moment yu (arbitrary except for zero spin), and it is shown that only for 
o=}, »¥1 (in particular, .~=0) and ¢=1, u=1 is the cross section of the correct magnitude 
and form (i.e., essentially independent of the mesotron energy) to account for observed burst 
phenomena at energies greater than 210" ev. Criteria for the validity of these formulae 
indicate that the region of applicability for the theory ¢=1, «#1 is more limited than that for 
o=1, w=1 (Proca). 
of 
e I. INTRODUCTION ciated with bursts of energy greater than 
: . : : 2X10" ev may be adequately described by 
y HE electromagnetic properties of nuclear x ; y q "yg thes 
; : assuming that the cross section for a given 
systems depend on the interaction of meso- ; 
1- ; ; ; fractional transfer of energy to the soft com- 
trons both with heavy particles and with the ce 
a ; ; ponent is independent of the mesotron energy. 
electromagnetic field. Indeed, so simple a quan- ; , , 
of ; ‘ . A cross section of this nature possessing the 
tity as the magnetic moment of the neutron is , ‘ , Le 
1. : correct magnitude would indeed be obtained"? 
completely determined only by both of these . ; ; 
y : : : ‘ if the electromagnetic properties of the mesotron 
interactions. Data obtained experimentally from eb 
e : . were those deduced from the Proca equations 
electromagnetic nuclear properties do not there- “iin pone" 
e ‘ ; by the application of the Born approximation. 
fore determine independently the mesotron- ; “as : 
S ; ; To examine the validity of the Born approxi- 
heavy particle coupling and the mesotron- : ; eee" 
e ne : ; ; mation, one may consider the limiting case of 
radiation interaction, and it would be of : . . 
e scattering by an electron of infinite mass, the 
advantage to study each of these forms of : : ; 
e ; : : process involved then being the scattering of 
interaction directly. Information concerning the ; ; 
s ; mesotrons by a static Coulomb field. An exact 
influence of the electromagnetic field on meso- ; : 
e : : se solution of this problem would enable’ one to 
trons may be obtained by the investigation of . a aa . . 
n ; x cs ; judge the reliability of the Born approximation 
recoil electrons resulting from collisions with the ; : 
; : as applied to this type of process. Unfortunately, 
mesotrons forming the penetrating component of ; : : 
€ ; “tee ; as will be shown in Section III, when one 
cosmic radiation. This process presumably ; 
yf , attempts to follow this program one encounters 
accounts for the fraction of the soft component : " : 
e , cai “- singular equations which admit of no complete 
(apart from secondaries arising from the disin- , "aaa 
: : ; set of regular solutions. This implies either that 
tegration of mesotrons in the atmosphere) which ; 
a ry: sae. ; the theory is wrong or that the model of the 
is in equilibrium with the hard component. The : : 
3 , scattering of a mesotron by a Coulomb field is 
| observed dependence of cosmic-ray bursts on : : 
; : : : too great an abstraction of the physical process 
their magnitude and the surrounding material a clei th emnamiinaiin tik setaienaiie 
? indicates that mesotrons may, with appreciable — — ‘ 
t probability, transfer a large fraction of their — 1 Oppenheimer, Serber and Snyder, Phys. Rev. 57, 75 
‘ energy to the soft component. Phenomena asso- (1940). 
y 2 H.S. W. Massey and H. C. Corben, Proc. Camb. Phil. 
e Soc. 35, 463 (1939). 
e * Commonwealth Fund Fellow. 3 A. Proca, J. de phys. et rad. 7, 347 (1936); N. Kemmer, 
h t National Research Fellow. Proc. Roy. Soc. A166, 127 (1938). 
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The electromagnetic properties of a particle 
with given charge and mass are essentially deter- 
mined by its spin and its magnetic moment. One 
may therefore use the experimental evidence 
available from cosmic-ray measurements to 
determine these two characteristics of the 
mesotron. Of the three simple possibilities for the 
spin (0, 4, 1) the case of zero spin may be 
excluded, for the impossibility of associating an 
intrinsic magnetic moment witha scalar mesotron 
implies an energy transfer cross section varying 
inversely with the mesotron energy, thus con- 
tradicting the experimental evidence. Were the 
spin of the mesotron 3, one could exclude the 
possibility of a magnetic moment of one mesotron 
magneton, as predicted by the Dirac theory, for 
this would lead to a cross section of similar 
character.‘ However, an alteration of the mag- 
netic moment changes the energy dependence of 
the cross section to that demanded by experi- 
ment. On the basis of cosmic-ray evidence, 
therefore, one cannot exclude the possibility of 
a mesotron of spin } and magnetic moment 
different from unity (in particular, zero) although 
such evidence as is available from nuclear phe- 
nomena indicates that this is not likely. 

The current theory of mesotrons postulates a 
spin of unity and a magnetic moment of one 
mesotron magneton. The physical adequacy and 
the mathematical difficulties of this theory have 
already been mentioned. The sole remaining 
simple possibility is that of a mesotron of spin 
unity and magnetic moment different from that 
assumed in the Proca theory. The following 
pages are devoted to the development of a 
general Lagrangian theory of such particles 
(Section II) and the application of the theory 
to the motion of mesotrons in a Coulomb field 
(Section III) and the problem of electron recoils 
(Section IV). 


II. GENERAL THEORY 


The most general bilinear Lagrangian density 
which involves only a four-vector ¢,, and the 
electromagnetic potentials A, in the form of the 
first-order gauge covariant derivative 


Dio,= [d,+ (ie/hc)A » lus 


*H. J. Bhabha, Proc. Roy. Soc. A164, 257 (1938). 
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together with the conjugate complex quantities 
$x, D4, is 
L=A.,(D,8")(D,o") + °9%e, (1) 


a quantity which automatically satisfies the 


. . . Py ae 
condition of gauge invariance. Here the A,, are 
numerical tensors which can depend only on the 
values of the metric 


uv=0 (uv), 
=1 (u=rv=1, 2, 3), 
=-1 (u=v=4). 


A general form of such a tensor is 
» ue ” ’ 
Zs =2 Ber + Berge + Reker (2) 


where f and ¥ are constants. The choice of unity 
for the first coefficient on the right-hand side of 
(2) is not an essential limitation. To guarantee 
the reality of ZL, one must impose the condition 


Ay. =A, (3) 
which implies that 8 and y are real. 


The equations of motion derived by variation 
of (1) are 


A‘,D,D,0" = Kb. (4) 
and the conjugate complex equations. Inserting 
the form (2) for the A. one obtains 


‘ D+Dyde 4 BD,D,¢” + yD.D,¢" = Ko, 
or 


D*Dyb.+ (B+7)D.(D,9") 
+ (te/hc)B$’Hy._=«°bs, (5) 


where we have used the commutation relation 
D,D,—D,D, = (te/he) Hy 


and the H,, are the strengths of the electro- 
magnetic field. 
Specializing Eqs. (5) to the field-free case one 
obtains 
OO ube t (6+ 1)9«(9,0”) = Kb. (6) 


Solutions of these equations may be written in 


the form 
de = We + 06 Y; 


where y satisfies the condition 0,¥°=0 and ¢ 
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is a scalar. Hence 
(#0,—«)¥e=0; (B+7+1)0d,9=K. (7) 


Writing «= Mc/h we may identify the first of 
Eqs. (7) as that satisfied by a particle of spin 
unity and mass M, and the second as that of a 
scalar particle of mass M(8+~7+1)-!. One there- 
fore has a theory of a particle which may exist 
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B+y+1=0, (8) 


which implies ¢=0. Were this restriction not 
imposed, there would be a finite probability for 
a particle to change its mass and spin in an 
external electromagnetic field, a phenomenon for 
which there is no experimental evidence. 

With this restriction, the coefficients assume 
the form 









in two states, with which are associated different 
spins and different masses. The theory may be 
restricted to particles of spin 1 by the choice 


uy uy a? ue, ar 
A,,-=(g Lor — Bre) + (Loks — 2rhe), (9) 
to which corresponds the Lagrangian density 





L=$P Fe + 6,¢+7 | (D,d*)(D,0’) — (D4) (D,0")}, (10) 


where F,,=D,¢,—D,@,. This Lagrangian may be simplified by noting that the expression in brackets 
may be rewritten 
— 0,(¢D,¢’) +9, (6*D,¢’) — (ie/2hc) Hy ( Go” — $4"). (11) 


The first two terms of (11) may be disregarded, for, possessing the form of divergences, they con- 
tribute nothing to the variation of the Lagrangian. We have therefore arrived at a convenient 
expression for the Lagrangian density, 


L=}P,,F" +6," — (ie/2hc)yHw(oo" — $4"), (12) 


which differs from that of current theory by the addition of a term involving the explicit appearance 
of the electromagnetic field strengths. This term may be interpreted as corresponding to an additional 
magnetic moment for the mesotron. 

The linear equations of motion which follow from (12) are 


D+F, = x, + (ie/hc)y¢"*Hy ; D,¢,—D,4, = F,,. (13) 


The electromagnetic current four-vector of the mesotrons which follows by variation of the 
Lagrangian with respect to —A, then assumes the form 


jt = (te/he)[6,F" — Fg, +-7d,(oo" — oe") ], (14) 


which differs from that of the usual theory by the term in y which is independently conserved. To 
show the interpretation of this additional term, we rewrite the expression for the current in the form 


j* = (ie/he)[(6,Do" — 6, Dd’) + (y +1) 0,(d*0" — 8°") + o*(D, 0") — &*(D.4’)]. (15) 


The second term represents a polarization current arising from the magnetic moment. Indeed, its 
coefficient is to be identified as the magnetic moment of the mesotron, which thus has on this theory 
the value 1+ mesotron magnetons. The first term has the appearance of a convection current yet 
may not receive this interpretation for it is not by itself conserved. It is, in fact, conserved only 
upon the addition of the terms involving D,¢@*—terms which have no simple physical interpretation. 
The importance of the proposed generalization now becomes apparent, for, by proper choice of + 
these terms may be made to vanish under certain conditions. 
One may evaluate D,¢’ by applying the operator D, to the first of Eqs. (13), with the result 


«Do = (te/2hc) (y — 1) Hy FY — (4aie/he*)yJ,9", (16) 
where J,, the current-vector of the external charges, has been introduced by means of the equation 
Hy» = —(42/c)J,. 





































A stills TILED salen 





956 H. C. CORBEN AND J. SCHWINGER 


Thus in regions outside those occupied by the external charges (i.e., J, =0) the choice y=1 implies 
D,¢’=0. It would then appear that the choice of y=1, which in the absence of external charges 
guarantees the subdivision of the current into independently conserved conduction and convection 
currents and imparts to the particle a gyromagnetic ratio of 2, places the theory in close analogy 
with the Dirac theory of the electron. 

This simplification finds its counterpart in the equations of motion 


D+D,¢,—D,(D,¢") - xo, + (te/hc) (y+1)¢H,, (17) 
which are obtained from (5) by writing .8=—(y+1). In comparison with the rather complicated 
equations of the current theory (y=0) viz.: 

(D+D, — x*)¢, = (te/he)o*H,, — (ie/2hc) (1/ x?) D, (Hye F*?) (18) 
one obtains, for y = 1, 
(D*D, — x*)o, = (2ie/hc)o"H,, — (4rte/he?x?) D,(J,¢*), (19) 
which in regions external to charge and current distributions reduce to 


(D*D, — x2), = (2ie/he)o"Hy». (20) 


To complete the theory, it is necessary to construct the stress-energy-momentum tensor and 
show that it satisfies all conditions that may be reasonably imposed on it. The real tensor con- 
structed by the recipe 


T*, =2{ (AL /0F ,.)Dbe+(AL/8 Fue) Dibe+(AL/OH yo) Hye} — 64L (21) 

satisfies the conservation equation 
0,7*,= —j*Hy (22) 
and thus may be interpreted as a stress tensor. For the Lagrangian (12), it assumes the explicit form: 
Tp = FueD ib? + F yeDr! — (ie /he)yHva( Gud" — Ou) — Burk. (23) 


In order to guarantee the conservation of angular momentum it is necessary that a theory yield a 
symmetric stress tensor. Although (23) is not symmetric, it may be replaced by the symmetric tensor* 


Or = Pye F iy? +P a Fu? +4°( Gib, + bby) + (te/he)y | Hou Grd" — br) + Hor( Gud" —$°bu)} — Burl, (24) 
which also satisfies the conservation equation 


0,.0%,= —j*Hy, 
since 
@+, = T+, — 9, { Fg, + F9g,}. (25) 


Although the energy density derived from this stress tensor is not positive definite, it does become 
so in the absence of an external electromagnetic field, which is all that may reasonably be required 
of the energy density derived from a theory of particles with integral spin. 


III. STATIONARY STATES IN A COULOMB FIELD 


The importance of a rigorous treatment of the equations representing the mesotron in a Coulomb 
field has been stressed in the Introduction. An attempt to construct exact solutions of these equations 


forms the basis of this section. ; 
The spherical symmetry of the Coulomb problem introduces the possibility of defining conserved 


5 This symmetric tensor may be also derived directly from the Lagrangian by the prescription 
OM, =2{ (OL/dF yo) Pye + (AL /OGy)br+ (OL/8 Fya) Fra t+ (OL/8by)br+ (OL /AH ye) Hye} — 5,*L. 
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angular momentum integrals. Appropriate definitions, in terms of the stress tensor 9,, are: 
Mum f Ou-x.8,0dr, (26) 

which, by employing the definition (24) of 0,, may be written 

My=(/e) f (Pal —i)(ee—x18))6'— Fu i(eiO.—n18)8\dr 


+(i/e+) f \(—) Pubi— Pus) -i(Fubs— Fad \dr. (27) 
The first term, with its associated operator 
L jp = —1(x 04 — X40 j) 
is obviously to be interpreted as the orbital angular momentum, thus endowing the additional term 


with the properties of spin angular momentum. The operator representing the total angular mo- 
mentum may therefore be defined by the equation 


J jx! = Lj! —1(8; Sim — 5¢'5jim)o™, (28) 
from which one infers the operational definition of J?, the square of the total angular momentum, viz. : 
J*o.= (L?+2)o.—21Li9/. (29) 


To express these formulae in a more amenable form, we may introduce space-vector notation, 
representing the space-components of ¢* by the vector ® and the time-component ¢* by the scalar 
function ig. The angular momentum definitions thus assume the guise of vector formulae: 


J ®=L,O+ie.X®, etc, J*H=(L?+2)®+2iLX®, (30) 


wherein e, denotes a unit vector in the z direction. The analogous definitions involving operations 
on the scalar ¢ are simply 
Jyv=L.ig; Je=L'e. (31) 


Stationary states of the system may be characterized by m and j(j+1), the proper values of J, 
and J?, respectively, and by the parity, the eigenvalue of the reflection operator R which is defined by 


R®(r)=@(—1r); Re(r)=—¢(-72). (32) 


This latter equation is a consequence of the opposite behavior of vector and scalar functions under 
spatial reflection. The construction of stationary state wave functions is facilitated by the con- 
sideration of auxiliary ®-functions which are also eigenfunctions of L* associated with the proper 
value /(/+-1). By virtue of Eq. (30), such auxiliary wave functions satisfy the condition: 


((G+1) —ll+1) —2)®@=21(LX®). (33) 


The detailed consideration of this equation is deferred to the Appendix, wherein it is shown that 
there exist three types of auxiliary ®-functions which may be symbolically expressed in the form: 


l=j: © =LF(r)P;”, 
l=j+1: @={(j+1)(t/N)+ilt/r) XL} FOP, (34) 
l=j—-1: @={—j(r/r)+i(r/r) XL} FO(n)P,”, 


involving three arbitrary radial functions. The general solution associated with the eigenvalue 
is obtained by linear combination of these auxiliary ®-functions. The existence of the parity quantum 
number permits a further classification of these solutions, for the parity associated with a state of 
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orbital angular momentum / is (—1)', thus assigning to the state /=j a parity opposite to that of 
the state /=7+1. The state /=j and a linear combination of the states = j+1 therefore constitute 
independent stationary states. No further sub-division will be possible in general, for the spin-orbit 
coupling introduced by relativistic requirements destroys the constancy of orbital angular mo. 
mentum. Turning to the scalar function ¢, we observe from (31) that 


= G(r) P;", 


introducing an additional arbitrary radial function G(r). For /=j7 the requirement that ¢ and © 
have opposite reflection characteristics can be made compatible with the equal parities of these two 
functions only by demanding that g=0. These remarks find their complete expression in the follow- 
ing formulae for the two distinct types of admissible solutions: 


 (@=LF)P @ = { (t/r) Fil) +i (4/1) XL] Fal) } Pi”, 
=]: l=jx+1: (35) 
¢=0, g=G(r)P;, 
which are written in terms of the linear combinations 
Fy(r)=(j +1) F(r)-—jFO(r),  Fe(r) = FO (nr) + FO(r). (36) 


These formulae express all the information available from the general symmetry properties of the 
system ; to proceed further we must resort to the specific dynamical equations. The equations (18) 
of current theory, specialized to represent the motion of a mesotron with energy W in a static field 


described by a scalar potential V(r), read 


{v= +(—*) -(—) |e e aad =| — a [te “I 
as Moa are ar} 1 


W+eV Mc edVr eh W+eV dV WteVr dg 
ea ie ° 











(37) 


—_ ——— — om 2 


he dr r M*c* = te ao fic or 











when translated into vector notation. These equations must be supplemented by the subsidiary 
condition contained in (16), viz.: 


W+eV eh dV(W+eVr deg 
--@-——}. (38) 


r= 
hc M*c* ae ic fr or 








V-@®+ 


This set of equations then yields the radial equations appropriate to the two types of solutions (35). 
The treatment of the state /=j affords no difficulty, for ¢=0, (r/r)-®=0, V-®=0 (L= —ir XV), 
thus effectively reducing the set of Eqs. (37), (38) to 


 foe(%EY-C)o-* 


d? 2d Kil+i) sW+eV Mey? 
[4 +(— ) -(=) | rn)=0 (38) 
dr? rdr h 


This is simply the canonical form of the Klein-Gordon equation, which requires no further attention. 
The radial equations for the second type of solution (/=j+1) are readily derived with the aid of a 
few elementary lemmata. The gradient operator finds its most convenient expression in the form: 


rd rxL 





or, in radial form 
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with the derivative relation 
oe 20 L? 
ail eee Caer (40a) 
or? ror fr 


Further 
r Tr 
L*®= “GHD t2Fit24G+1) Fe) +i(- XL) GG+1)+2F) Pi, 
r r 


which is most easily obtained by reverting to the original expression of ® in terms of eigenfunctions 
of L?. It is now possible to write the terms of the first of Eqs. (37) as linear combinations of the two 
angular functions (r/r)Pj" and i[(r/r) XL]P,", the coefficients of which yield the desired radial 
equations for F; and F;. The radial equation corresponding to the latter part of (37) may be written 
down by inspection, while to perform a similar operation for (38) requires only the observation that 


r r a 2 jG+1) 
V-@= VF Hiv--XL Fs} Pn=| (—+-) Ft Pal Pm 
r or r r 


r 





The complete set of radial equations thus obtained is: 


o=—Pt JG) 6 —* “(0 “), 





dr hc dr 
- le i dV on € ail e W+eV dV 
wap <am, GG anf nem g-, 41 
- * her a cdr tc he ar -_ 


dF; 2. j(j+1)_ W+eV e dV 
+> +—_—G 





—-+-F, 2 =— 


dr fr r he he dr’ 
where 2 symbolizes the Klein-Gordon operator 
dq 2d j(j+1) (= a) -(=) 





Q=—+4-—- 


dr? rdr r? oa 


and Q is defined by 


~) on hc 





Mc dG W+eV 
( Fi. (43) 


It should be noted that if 7=0 the second equation of the set (41) disappears, for this equation was 
inferred from the coefficients of i[(r/r)XLJP,;", which identically vanishes for zero angular 
momentum. 

In order that regular solutions of these equations exist for the Coulomb field, V = Ze/r, it is neces- 
sary that the functions F;, F; and G vanish at the origin at least as rapidly as some power of r, e.g., 
r, so that, from Eqs. (41), Q@ must exhibit a radial dependence proportional to r‘+". To test the con- 
sistency of these requirements with the definition of Q, we may employ the equations obtained by 
substituting the fourth member of (41) into the first and third members, namely: 


W+eV h \* j+1) 
-ro( i) 
hc 2 


dQ h \* j(j+1) W+eV 
7 t,0-6 +(x) r {c+ hc Po}. 
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In states of nonvanishing angular momentum (j#0), the second terms of these expressions for Q 
are dominant near the origin, thus imparting to Q a radial dependence two powers of r less than that 
of the functions F;, F2, and G, contrary to the mathematical regularity requirements. Therefore, no 
regular solutions exist for states possessing finite angular momentum. Whether any solutions com- 
patible with the physical requirements exist remains an open question. For j7=0, these Eqs. (44) 


become simply aay 40 2 
e 
Q; G=—-+-Q, (45) 


hc dr r 


= 





which, coupled with the definition (43) of Q, imply that 


d*Q Phan +(—)- (*) l= . a 


dr? rr ae 








the Klein-Gordon equation for a state of unit angular momentum. This result is quite satisfactory, 
for a state of zero total angular momentum is rigorously a state of unit orbital angular momentum, 
The general conclusion to be drawn from these considerations based on the Proca theory may be 
summarized by the statement that regular solutions exist for the Coulomb field only for those states 
in which the orbital angular momentum is a rigorous constant of the motion, namely /= 70 and 
I=1, j=0. 

The technique developed for the expression of the Proca equations in radial form is immediately 
applicable to the equations of the modified theory (y=1). In vectorial notation these equations 
comprise the following set: 


W+eV Mc 2erdV 4neh 
{v4 +(— )- -(= —) |e =—-—-——y— V(T(r)¢), 




















her dr M?*c3 
W+eV Mc 2edV .. 4re 
|v +(— )- -(= —) Je= macy (W+eV)Tine. (47) 
he drr Me! 
W+eV 4reh 
V-@+ v -™- 
he M?c3 


wherein I'(r)=J*/c denotes the density of the static, spherically symmetric charge distribution 
under whose influence the mesotron moves. The treatment of the state /=j7 need not detain us, 
the Klein-Gordon equation being obtained as before. The radial equations describing the coupled 


states /=j+1, viz.: 











2 zt 2edV 4neh d 
QF, =—(Fit+j(j+1) F2) -— —G- —(T(r)G), 
r? he dr M?*c3 dr 
2 4reh 1 2edV 1e 
QF,= Ait -Ir(r)G, 2G=——F,- (W+eV)I(r)G, (48) 
r? M*c8 r hc dr M*c! 








dF, 2 i(j+1) W+eV 4reh 
sti G r 


el ak 2 is 


dr r hc M?c? 


may be obtained, mut. mutand., from the previous set (41). A complete solution of these equations 
can be effected in regions of space unoccupied by external charges (I'(r) =0). Within such domains, 
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the radial equations assume a particularly simple aspect : 


QF, =(2/r)CFitj(j+1)F2t+ZaG], QF,=(2/r?)F:, G=(2/r?)(—ZaF), 





dF, 2 ji+1) W+eV 
eens men | F:.+———-G =0 (49) 
dr rr r hic 


when the Coulomb field (eV =Ze?/r=Zahc/r) is introduced explicitly. A process of diagonalization 
reduces these equations to a canonical form. 
The linear combination | 
AFi+jG+1)F2+ZaG (SO) 


provides a solution of the differential equation 


OF = (20/P)F 


if \ is chosen to represent one of the roots of the eigenvalue equation 


AMA—1) =jG+1) —Z?a? (51) 
viz.: 
Ar=$+([(j+3)?—Za?}!; Ae=$—L(G+3)?-—Z?a*}}. (52) 
Associated with these roots, therefore, are two functions: 
(Ar —A2) Hr =A Fitj(j+1)F2+ZaG, (A2o—A1) 2 =AeFit+j(j +1) F2+ZaG, (53) 


which satisfy the differential equation 


d? 2d W?—M%! 2WZa dA+1) 


dr? rd? h?c? he r r? 


Git, t:=t,) 


the form assumed by the Klein-Gordon equation in a Coulomb field. The two Eqs. (53), supplemented 
by the fourth member of (49), viz.: 


F,=Git+52, j(j+1) F2+ZaG= — (A2Fi1+A1F2) 


hc({ d 1+ A, d 1+): 
G=—-—}—Sit Fit—Fet+ 
Widr r dr r 


,= 





(54) 


(S5) 








Fer, 


constitute a complete solution of the problem in terms of §; and §e. 

In states of nonvanishing angular momentum, there exist two linearly independent solutions 
associated with §§; and with §2, which we label by a quantum number “‘/’’ = 7+1, j—1, respectively. 
In the nonrelativistic limit, ‘‘/’’ becomes the orbital quantum number. §; and 2 may not be con- 
sidered as independent in states of zero angular momentum, for a relation between them is implied 
by the Eqs. (55) specialized to the state 7=0. The last two equations—equations for G—now serve 
to determine the arbitrary constants in §2 in terms of those in §: thus providing only one independent 
solution. To show this in detail, we observe that, if {§, is a solution of (54), 


d 1+’ W Za 











F.u=—Ft+ F,—-— —F,, (56) 
r r fic X 
whence 
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in which the relation \;A2=Z?a? (j7=0) has been utilized. Since \1(A1— 1) = A2(A2—1), Fa,-1 and Sas 
satisfy the same differential equation, thus enabling one to impose consistently the condition 


Oa,-1+9a,-1=0, (58) 


which exhibits the relation between j§; and §2 demanded by the coexistence of Eqs. (55) for j=0, 

Although the equations of the modified theory (y=1) are completely soluble in regions devoid 
of external charges, a comparable situation does not exist within charge distributions (I(r) #0), 
While the equations presumably possess regular solutions in regions of smoothly varying charge 
density, no serious attempt has been made to obtain them. In the limiting case of a point charge, 
the presence of the singular terms in I'(r) imposes boundary conditions, which may be inferred from 
the requirement of the existence of physically admissible charge distributions. According to Eqs, 
(15) and (16), the mesotron charge density p is given by 


2e | W+eV _— - 4reh 
p= ———1 _ —_(®- @+ $e) —(®-V)$—(®-V) ¢+——I(r) Ge}, 
hel he M?*c3 


which would predict an amount of charge 
—2Za(eh/M*c*)G?(0) (59) 


to be contained within an arbitrarily small neighborhood of a point charge Ze. For solutions of the 
type j= ‘‘l’’+1, G(r) ~r* in the vicinity of the origin, while for those of type j= ‘‘/"’—1, G(r) ~r' 
Since \; is not an integer, G(r) tends to either zero or infinity as 7 is diminished. In the latter case, 
which occurs for j7=0, “/’’=1, and j=1, “7’’=0, the amount of charge (59) in the immediate neigh- 
borhood of the origin, and hence also the total charge of the system, becomes infinite. This meaning- 
less situation implies that these two states are completely forbidden. 

It is of some interest to note that the energy levels of the allowed discrete states are represented 
by the Sommerfeld formula, all states of given total angular momentum and principal quantum 
number coalescing as in the Dirac theory. This result is inexplicable on the basis of the formula 
derived from the Larmor and Thomas precessions, which would predict in addition a dependence 
on the orbital angular momentum quantum number /. The origin of this discrepancy is to be sought 
in the inadequacy of considering the spin as a simple vector, for the coupling energy between spin 
and orbit consists not only of the usual Larmor-Thomas terms but in addition terms which are 
nonlinear in the spin. The equivalence between the Larmor-Thomas formula and the corresponding 
prediction of the Dirac theory is to be attributed to the necessary linearity of the spin terms for a 
particle of spin 3. 

The general program of this section has been an attempt to obtain rigorous solutions of the equa- 
tions for a particle of unit spin in a Coulomb field. While for the equations describing Proca particles 
this attempt has been abortive, the modified theory discussed in this paper provides regular solutions 
for all states save two, these being completely forbidden. The net result of the modification has been, 
therefore, the concentration rather. than the removal of the singularities. The lack of a complete 
set of regular solutions thus compels us to have resort to the Born approximation in order to give 
meaning to these two theories. 


IV. ELECTRON-MESOTRON COLLISIONS 


Information concerning the electromagnetic properties of mesotrons may be obtained by the 
study of various scattering processes. However, it would be difficult to infer such information from 
processes such as mesotron-proton scattering and radiative emission induced by mesotrons, which 
involve in addition the mesotron-heavy particle coupling. For collisions between mesotrons and 
electrons, however, this latter type of coupling is presumably absent, and the energy transfer is 
dependent only on the electromagnetic interaction between the particles. Using various mesotron 
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models, we shall calculate the cross section for this energy transfer, developing a general analysis, 
based on the Born approximation, without use of the explicit form of the mesotron charge and 
current densities. The results are applied to mesotrons of spin ¢=0, magnetic moment p=0, o=}, 
yu arbitrary (in particular 4=0) and o=1, uw arbitrary (in particular .=1 and w=2). 

Let us consider a transition of the mesotron from a state of energy W® and momentum P*® to a 
state of energy W and momentum P. With this transition are associated charge and current densities 


p= pis exp [i (P°—P)-r—(W°—W)t]/h}, j=jisexp {i[(P°—P)-r—(W°—W)t]/h} (60) 
which produce a transition electromagnetic field 
4rh® exp {i[(P°—P)-r—(W°— W)t]/h} , _ exp {iL (P°—P)-r—(W°— W)t)/h} 


= ji V =4rh*p; . (61) 
cc) [PO P22 W)? ‘| P°—pi2—--2(W— W)? 








This field induces transitions of the electron from a state of energy and momentum E°, p® to the state 
E, p. The spin direction of the electron in these positive energy states may be characterized by the 
indices \°, \=1, 2. If the electron wave functions are taken to be 


initially : up exp [i(p°-r—E%)/h], finally: up* exp [i(p-r—Et)/h], (62) 


the condition (up»*, up”) = (up, up’) =1 corresponds to one electron per unit volume. The transi- 
tion probability per unit time from the state (A°, p®) to the state (A, p) of the electron is then given by 


1 
(w (o1;—~2-ju)ue™ ) 
Cc 


To find the total electron transition probability, we sum over the final spin and momentum states 
of the electron. The summation over the final momentum states has merely the effect of imposing 
the condition of momentum conservation. Summing over the final spin directions, one has 
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-3(W+E—W°—E°)ép+p, Po+ps. (63) 
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c 


which is independent of the initial spin of the electron. To calculate the cross section for a given 
incident mesotron, we sum the transition probability (63) over all relevant final states of the mesotron 
and divide by the product of the mesotron particle density and the relative velocity of the electron 
and mesotron. The cross section thus defined is relativistically invariant and is most easily computed 
in the center of mass system (P°+p°=0), since in this system the energy of either particle does not 
change during the collision. The resulting cross section for the scattering of mesotrons through an 
angle 3 into a solid angle dw is 


9 


d f - 70 0.4 2 2 po? .:.,2 9 JF e 
=a wee) EG pist+P *jis)?—c P®™ sin? 19(61° erie) 


dw 
sin‘ $3 








(65) 


if the mesotron particle density is also normalized to represent one particle per unit volume. Here we 
have used the facts that the relative velocity is (c?P°/W°)+(c?P°/E®) in this coordinate system and 
that the number of mesotron states per unit volume and per unit total energy range is 
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For a particle of spin one-half and magnetic moment unity the charge and current densities are given, 
according to the usual Dirac theory, by the four-vector 


je = — (eh/2M) { Gd" — (0"b) ¥) +39,(b(y"7’— 7") ¥)}. (66) 


Here ¢ is defined by ¢=yt*y‘, and y* is given in terms of the more familiar Dirac matrices a, 8 by 
y= —iB, y‘= —iBa‘. The two terms on the right-hand side of (66) are independently conserved and 
correspond, respectively, to a convection current and a polarization current arising from the mag- 
netic moment. In close analogy with the theory for particles of spin unity, outlined in Section II, 
one may develop a theory for which the coefficient of the polarization current is arbitrary, cor- 
responding to a particle of spin } and of arbitrary magnetic moment. We consider here the particular 
case for which this magnetic moment is zero. The charge-current four-vector then assumes the 


simple form 


j* = (ieh/2M)[y* Bay — (a*y*) By], (67) 


an expression which differs from the corresponding result for particles of spin zero only by the presence 
of the matrix 8. Using plane waves for the mesotron wave functions, viz.: 


¥i= up exp [i(P*-r—W%)/h], ws=ur* exp [i(P-r—Wt)/h], 


one finds that the transition charge and current densities corresponding to a mesotron of spin 3 and 


magnetic moment zero are 
(W°/Me!)(we>, Bue), jeg — (ur, cr) 

pis=—e c*)(up*, Burp’), jis = —e— —( uP’, Burp’). 68 

’ : 2W) Me? 


The wave functions are normalized by the condition 


(up, up”) = (up, up) =1, 


which corresponds to one mesotron per unit volume. 
Using the relation 








Mc! Pp 
x | (ur, Bur’) |?= 1+ sin? 38 (69) 
A=1,2 w” M?c? 


/ 


and substituting (68) in (65), one finds that the cross section for the scattering of a mesotron through 
an angle # into a solid angle dw is 


e! We \* do [ ap): 
oles) seller) 
4ctP*\ yr+k/ sint tab we 


2 cPY, Mc p® 
—c?P® sin? 10| (24 )+ t].[ + sin? 3 | (70) 
Ww Ww w” M*c? 











The corresponding cross section for a particle of zero spin—and hence zero magnetic moment— 
differs from the above expression only by the absence of the last factor, which, for the case of spin }, 
arises from the spin summation. It is to be noted that at high energies this factor entirely alters the 
dependence of the cross section on the mesotron energy. 

To obtain from (70) the cross section for Coulomb scattering it is only necessary to consider the 


limiting case of infinite Z°, which yields 
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in which we have set W°=eMc’. At high mesotron 24 (71) becomes 


do=- -(— ~)— (72) 
Mc*/ sin? 1g 


which differs from the corresponding formula deduced from the Proca theory® by the absence of 
a factor } cos? $0. 

The above formulae represent the angular scattering cross section in the barycentric system ; we 
are primarily interested in the cross section for a given fractional energy transfer to an initially 
stationary electron. It is therefore expedient to introduce the ratio, f, of the kinetic energy gained 
by the electron to W)>=eMc’, the total energy of the incident mesotron, measured in the laboratory 
system. The relation between f and the angle of scattering in the center of mass system is simply 


2w(e?— 1) 


f=fmax sin? 48, where fue Perry (w=m/M) (73) 


represents the maximum fraction of the mesotron energy which may be transferred to an electron 
by impact. Employing the relation between W, and the electron and mesotron energies in the bary- 
centric system, viz.: 








w(e+w) 1+ we 
E%= Mc?, W= Mc*, (74) 
(1+ w?+ 2we)! (1+w?+ 2we)! 


we may rewrite the cross section as follows: 








deo=r 








eae 7 
— e[s0- p+—-s-49)-— J] (75) 


At high energies this is independent of the energy of the incident mesotron. For a mesotron of arbi- 
trary magnetic moment »#1 and spin 3 the leading term of this expression is multiplied by a factor 
(u—1)*. The corresponding result for mesotrons of spin zero is obtained by omitting the last factor 


in (70): : 
do=nr —(1—f)-——f I, 76 
(—) = e—1 f? —<|- a-/) we? 76) 


which for high mesotron energies e decreases inversely with e. 
For mesotrons of spin unity, one again uses plane waves for the mesotron wave functions, viz.: 


Cc y= 
2Ww? w® 











4 
o,=ne(—) e° exp [7(P°-r—W%)/h],  oi= — ihe exp [7(P°®- r— W%) /h] 


with similar expressions for the final state. The normalizing condition 
e®- eo — (cP°- e°/W)?= 1 


then corresponds to one mesotron per unit volume. 
The transition current and charge densities are given by 


cP°-e° cP-e cP-ec(P®°—P)-e° cP®-e® c(P°—P)-e 
pis= ~e| (e?-e- )+i0+n/ - )} 
we w w w w w 
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TABLE I. Values of fx(f) for particles with various spins 
and magnetic moments which collide with the electron. 











TYPE SPIN MAGNETIC MOMENT K(f) 

I 0 0 2n-(1—f) 

I] 1 : 2n-*(1—f+4f%) 
Ill 3 ul (u—1)*(1—f) 
IV 1 1 $(1—f+4f?) 

Vv 1 ul $(u—1)*n(1—f) 








and on substitution in (65) one may obtain the 
cross section. For general y and e the expression 
obtained is complicated but for large « and y #0 
the major contribution to the cross section arises 
from the terms proportional to y? which cor- 
respond to longitudinal-longitudinal transitions. 
One finds, for fwe>1, f< fmax=1—(1/2we), and 
70, 

do = $x(e?/ Mc*)*y*we(1 — f)df, (78) 


which is to be compared with the corresponding 
cross section for y=0 


do=jn(e/Mc)*(1—f+3f)df. (79) 


It is to be noted that for y#0 the cross section 

increases linearly with the energy of the incident 
mesotron. 
Writing 

da=aox(f)df, oo=x(e?/Mc*)? (80) 


for the cross section for a particle of energy 
«Mc? and mass M to transfer energy feMc* to an 
electron of mass m=wM, we may summarize 
the above results for 7=few>1 in Table I. 
Only cases III and IV lead to cross sections 
essentially independent of the mesotron energy, 
as required by experiment. In addition, it has 
been shown! that in order to account for the 
observed size’ of bursts of energy >2X10"° ev 
it is necessary for the cross section to be such that 


.= J (f)fetdf~4. 


For the case of spin } and magnetic moment u.¥1 
«=0.2(u—1)?, 
which is of the correct order for |w—1| <2. 


7M. Schein and P. S. Gill, Phys. Rev. 55, 1111 (1939); 
M. Schein and V. C. Wilson, Rev. Mod. Phys. 11, 292 
(1939); H. Carmichael and C. N. Chou, Nature 144, 325 
(1939); Bhabha, Carmichael and Chou, Proc. Ind. Acad. 
Sci. Al0, 221 (1939). 
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Because of the uncertainty in the experimental] 
evidence the value of » cannot be assigned more 
definitely, but it is to be noted that, for the 
particular case 4=0 the value of «x is 0.2, which 
is very close to the corresponding value (0.22) 
obtained from the Proca theory. It is therefore 
impossible to distinguish between cases III and 
IV from cosmic-ray evidence alone. 

All these results having been obtained by use 
of the Born approximation, it is important to 
examine the conditions under which they may be 
expected to have validity. A necessary condition 
for the validity of the Born approximation is the 
smallness of the coupling between the mesotron 
and the fields arising from both the electron and 
the zero-point fluctuations of the electromag. 
netic field. A previous investigation of these 
questions by Oppenheimer, Snyder and Serber! 
has disclosed that the more incisive condition 
arises from the requirement that the coupling 
between the mesotron and the zero-point fluc. 
tuations be small. To express this requirement in 
quantitative form, we consider the mesotron- 
electron collision in the center of mass system, 
employing wave packets of approximate linear 
dimensions h/P. The kinetic energy density 
~Pc(P/h)* must then be large compared with 
the coupling energy density —(1/c)j-A. The 
potential A describes the zero-point fluctuations 
of wave-length ~h/P and possessses the order of 
magnitude P(c/h)!. The form of the current 
density j depends on the assumed mesotron 
model and is of the order of magnitude ec(P/h) 
<(P/Mc)", where n assumes the value zero for 
particles of type I, II of Table I, n=1 for types 
III, IV, and »=2 for type V. This varied mo- 
mentum dependence of the current density is of 
course intimately related to the different energy 
dependence of the cross sections for these several 
types of particles. The condition of smallness of 
the coupling energy density compared with the 


' kinetic energy density is then expressed quan- 


titatively by 
ai(P/Mc)"<1. 


For the various types of particles this condition 
assumes the form 


(I, II): al<, 
(III, 1V): Wo<M?%c?/am=2X10" ev, (81) 
(V): Wi<M?2c?/ atm =2 X10" ev 
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in its limitation on the mesotron energy Wo. The 
domain of validity of the modified theory (case 
\) is thus more limited than that of the Proca 
theory, but it is still of sufficient extent to include 
most of the experimental region. 

The authors wish to express their gratitude to 
Professor J. R. Oppenheimer for his continued 
advice and encouragement. 


APPENDIX 


The basic equation determining the simultaneous eigen- 
functions of L? and J?, Eq. (33) of the text: 
((i(J+1)-U+1)—2)@ =2i(Lx ®), (Al) 
may be readily solved by elementary quantum algebra. 
The operations of scalar and vector multiplication with L 
yield, respectively, 
(jG+1) -10+1))L- ®=0, 
Ki+1)®—L(L- ®) =}(7(j+1) -1l+1)) 
X(7(G+1)-10+1) —2)®, 
from which one obtains the characteristic equation 
(jG+1)-—10+1))(GG+1) —H+1)? 
—2j(j+1)—210+1))=0 © (A3) 
by elimination of L-@®. The admissible values of the 
orbital angular momentum are thereby restricted to 


(A2) 


b=j—-1, j, j+1, (A4) 


a not altogether unexpected result. 

Sufficient information is available from the relations 
(A2) to determine the properties of the solution /= 7. The 
equation thus obtained: 


jJG+1) @O=L(L-®), (AS) 


when reformulated in terms of the scalar function u, 
defined bv 


L-@®=j(j+1)u (A6) 
yields the general form of ®, viz.: 
@M=Lu, L*u=j(j+1)u. (A7) 


The second equation of the set (A7) implies that u is a 
spherical harmonic of order j, Pj", multiplied by an arbi- 
trary radial function F(r). The spherical harmonics are 
assumed normalized to provide the conventional matrix 
representation of angular momenta. One therefore obtains 


l=j: ®=LFir)P,". (A8) 

The choice of the spherical harmonic P,” guarantees that 
J, shall have the eigenvalue m, for 

J &,=L&,=m?,,. (A9) 


The relations thus far developed (A2) provide only the 
information that 
L-®=0 (A10) 
when restricted to the states /=7+1. To proceed further 
one may utilize the equations 
({(j+1)—10+1)+2)(r/r)-@ = —2L-((r/r)xX®), 


GG+1)—10+1))(2/r) X ® =2iL((r/r)- ®), (Alt) 
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obtained by applying the operations of scalar and vector 
multiplication with r/r to the basic equation (Al). The 
introduction of the scalar function v, defined by 


(t/r) -@= — 3(j(f +1) -10+1))o, (Al2) 


permits a more convenient expression of Eq. (A11), namely 
(r/r)X ®= —il», 


Ly =}3(j(j+1)-10+1)) 
X(GG+1)—-Hl+1)+2)e (A13) 
= j(j+1)v. 
The relations (A12) and (A13) suffice to determine the 
general form of ®, for 


® =(r/r)((r/r)- ®)—(r/r)X[(t/r)X ®] 
= {4(1(2+1)—j(j+1))(t/7)+2((r/r) XL)}v. (A14) 


In virtue of the second of Eqs. (A13), the scalar function 
v may be represented as the product of the spherical har- 
monic P," with an arbitrary radial function which is 
independently assignable for the two states /=j+1. As 
before, the choice of the harmonic P;” guarantees that J, 
shall have the eigenvalue m. The ®-functions associated 
with these two states, expressed in terms of the two 
arbitrary functions F“(r), F®@(r) then assume the sym- 
bolic form 


b=jti: B= {(j+1)(t/r)+4(t/r) XL} POW) P;", 
l=j—1: ®={—j(r/r)+i(r/r) XL} F@(r) Pi. (A15) 
To exhibit clearly the angular dependence of the ®-func- 


tions for the several states, we shall write in detail the 
behavior of the components 


©, +10, 
, = id, ’ 
®, 





ta (CG—m)(G+m+1) PP?) 
! C(j-4+-m)(j—m+1) PP LF), 
m Ps 
b= j+1: (A16) 
2j+1\! ~[(j-+1-4+m)(j+2-+-m) PPT | 
note [(j+1—m)(j+2+m) PPpn FPO: 
L Cij+1}*—m*} = Pm 
l=j-1: ; 
2j+1\! [(j—m)(j—m—1)} Pi ' 
et ) —[(j+m)(j+m—1)PPpa fh): 
L [j?—m?*}} Pj 





These results are not completely novel, being intimately 
related to the angular representation of electromagnetic 
multipole fields.* The electromagnetic equations are 
identical with the field-free Proca equations in the limit 
of zero rest mass, the quantities ® and ig being identified 
with the vector and scalar potentials in this limit. The 
vector ® satisfies the differential equation 


(V2+k)@=0 (k=W/he), (A17) 


8 W. Heitler, Proc. Camb. Phil. Soc. 32, 112 (1936); 
W. W. Hansen, Phys. Rev. 47, 139 (1935); S. M. Dancoff 
and P. Morrison, Phys. Rev. 55, 122 (1939). 
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multipole of order 27. Using the recursion relations 
isa t- Sia = (25+1/kr)fi, (A19) 
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which, for our purposes, is more conveniently expressed as 
@ 2 d Poy ©=0, 
& rdr 


employing Eq. (40a) of the text. All three possible types 
of radial functions associated with a given orbital angular 
momentum / thus satisfy the differential equation 


@ 2d l(i+1) 
tg thio. (ANB) 
dr? 
The solution ®; for 1=7 (or more exactly the real part of 
Det"); 

@;=Lf,(r)P" 


represents a transverse wave field (V-@;=0) which is to 
be identified as the vector potential describing the radi- 
ation emitted by a magnetic multipole of order 27. Ap- 
propriate linear combinations of the solutions for /=j+1 


Din. = (G1) (8/1) +4((2/7) XL) } fig P, 
@j_1= | —j(t/r) +4((t/r) XL) } fj-aPs, 


may be found to represent the longitudinal and transverse 
parts of the vector potential associated with an electric 


(+1) fi+1— Ffi-1= — (2j+1/k)(d/dr) f;, 


one may verify that 








2jt+icrd rxL 
@1= Bir @1= -— f-fim (A20) 
k \Wwdr 
= —(2j+1/k)vfP;", 
by Eq. (40) of the text. Similarly 
j j+1 
®, — 2 ®;,, woke ®;_; 
2j+1 2j+1 
_fSIGt) r rxLi 
= =fj—i — <r) Pp (A21) 
ae r 
= —(i/k)y X(Lf,P;). 


These functions therefore satisfy 
VxX@,=0, V-®,=0, 


and are thus to be identified as the longitudinal and trans. 
verse fields of an electric 2j-pole. 
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In Section I a general formula for a, the coefficient of recombination of ions in gases, is de- 
veloped. It covers all ranges of pressures and temperatures and can be made to include all types 
of recombination (preferential, initial and volume ionization). The evaluation of the general 
formula, (1.1), depends on the relative values of three linear quantities: the mean free path of 
the ions \, the mean distance between ions of different signs # and the well-known parameter 
ao=e@/(RKT). In Section II the case Xap is treated. The mechanism of recombination then 
depends on the ratio ao/7. If ao/f is large the migration of the ions under their mutual attraction 
prevails and (1.1) leads to Langevin’s formula; in the opposite case, ao/r<1, diffusion is the 
decisive feature and 1.1 leads to a formula which is practically identical with that of Harper. In 
Section III the case \>>ap is treated by a method previously developed by the author. Under 
certain restrictions (1.1) reduces to Thomson’s formula, but in general a depends on the concen- 
tration of the ions. In Section IV it is shown that (1.1) is in fair agreement with such experi- 
mental data as are available for the region of transition between the cases treated in Sections I] 


and III. 


HE problem of recombination of ions in 
gases under varying conditions of ioniza- 
tion, pressure and temperature has proved to 
be much more complex than was originally 
anticipated. Various types, such as preferential, 
initial and volume recombination are involved.! 


1 For the definition of these types and a survey of the 
Fe subject see the excellent treatment in: L. B. Loeb, 

Fundamental Processes of Electrical Discharge in Gases 
(John Wiley and Sons, New York, 1939). 


It has been recognized in recent years that the 
two most important theoretical formulae for the 
coefficient of recombination, that of Langevin 
and that of Thomson, have separate domains of 
applicability. Harper? and Loeb* have given 
formulae which bridge the expressions of 
Langevin and Thomson in a formal way, but it 


? W. R. Harper, Phil. Mag. 18, 97 (1934) ; 20, 740 (1935). 
*L. B. Loeb, Phys. Rev. 51, 1110 (1937). 
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seems desirable to establish a general theoretical 
expression from which either of them, and 
eventually others, can be derived. It is the object 
of the following paper to give a theory that 
covers the whole range of physical conditions 
and that can be made to include the various 
types of recombination. 


I. THe GENERAL FORMULA 


We start by defining three probabilities. Let 
n, and n_ be the instantaneous average volume 
densities of positive and negative ions, respec- 
tively. The ionization is considered to be homo- 
geneous in a macroscopic sense and no preferen- 
tial directions are assumed to exist. Still, in order 
to include preferential and initial recombination, 
we shall admit microscopic inhomogeneities. 
Then let p:(r)n_dr be the probable number of 
negative ions which are located in an element of 
volume dr at a distance r (irrespective of direc- 
tion) from a positive ion. 

By p2(r)dr we define the probability that the 
nearest neighbor to a positive ion is to be found 
at distances between r and r+dr. Evidently we 


have to have 
f p2(r)dr=1, 


if 7m is the nearest distance of approach for two 
ions of different signs. 

Further, let »3(r) signify the probability of 
recombination. That is, of dN, pairs of nearest 
neighbors of different signs which are at distances 
between r and r+dr, p3(r)dN, are supposed to 
recombine on the average. This probability 
concept replaces the strict criterion for recom- 
bination which other authors‘ have tried to 
establish. 

Finally, we have to define the velocity of 
approach of two nearest neighbors of different 
signs. If they are at a distance r apart, u(r) is to 
signify the average component, in the direction 
of r (from the negative to the positive ion), of 
the relative velocity. This is supposed to be the 
resultant of all acting influences, i.e., diffusion 
and electrostatic attraction. It might also be 

* (a) J. J. Thomson, Phil. Mag. 47, 337 (1924); (6) L. B. 
Loeb and L. C. Marshall, J. Frank. Inst. 208, 371 ti932; 


W. R. Harper, Proc. Camb. Phil. Soc. 28, 219 (1932 
31, 430 (1935). 
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made to include the action of an external field, 
but we limit ourselves to the case where there is 
no such external field. 

The expressions of m,, m_, p: and p, will, in 
general, be functions of the time, as it is our 
object only to calculate the instantaneous coef- 
ficient of recombination, a. 

With these definitions it is easy to write down 
a general expression for a by generalizing 
Langevin’s argument. In doing so we first fix 
our attention on such positive ions as have their 
nearest negative neighbors at distances between 
r and r+dr. There are p2(r)n,dr positive ions of 
this class per cc. We consider each of them sur- 
rounded by a sphere of radius r. The number 
of negative ions crossing one of these spheres 
during a time interval Af will be given bv 
4nr*pi(r)n_u(r)At. But not all these ions will 
ultimately move up to the positive center and 
recombine, but only the fraction p3(r) of them. 
Hence 


4nr’n,n_pi(r)po(r) ps(r)u(r)drdt 


is the number of pairs (per cc) of the specified 
class which are to be considered as lost by 
recombination during At. Adding over all dis- 
tances r and dividing by n,n_At we obtain 


a=4r f Pilr)pol(r)ps(r)u(r)r’dr = (1.1) 


as our general expression for the coefficient of 
recombination. 

P Before we try to evaluate (1.1) for special 
cases, a few words must be said regarding the 
physical significance of the various functions of r. 
First of all, the probabilities p:(r) and 2(r) 
determine the particular type of recombination. 
They are not independent of each other. 

If the distribution is homogeneous also in the 
microscopic sense, i.e. if ~:(r)=1, and, further- 
more, if the diameter of the ions can be neglected, 
the probability 2(r) is known. P. Hertz® has 
shown that, in a random distribution of points 
without extension (m per cc), p2(r) is given by 


b2(r)dr =4xnr* exp (—4anr*/3)dr. (1.2) 


This law is modified, though only slightly, if, 
instead of points without extension, ions of 


* P. Hertz, Math. Ann. 67, 387 (1909). 
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finite diameter have to be considered. It is, 
however, evident that this modification will be 
insignificant if the average separation of the ions 
is very large compared with their diameter, as 
in all cases of gaseous ionization. Then, also, rm 
can be replaced by 0. 

If p1(r) is a given variable function of r, it can 
be shown by an argument which is exactly 
analogous to that of Hertz that 2(r) will be 
determined by 


p2(r)dr=4xr?n_pi(r) exp (—n_P(r))dr, (1.3) 
where 


P(r)= | 4xr*pi(r)dr. (1.4) 


'm 


Thus when /;(7) is determined the probability 
p2(r) may be calculated. Bradbury*® has recently 
determined ;(r) for the cases of preferential and 
initial recombination and, therefore, it would 
not be difficult to evaluate a for these types of 
recombination. However, we shall limit ourselves 
to the fundamental case of volume ionization 
with n,=n_=n. Not even then is it permissible, 
in the strict sense, to replace p:(7) by 1. Owing 
to the electrostatic attraction there is a greater 
probability of finding a negative ion in an 
element dr near a positive ion than in an equal 
element at a greater distance. This is conclusively 
shown by the theory of ionic atmospheres as 
developed by Debye and Hiickel.’ This theory 
is applicable to our case, but for all cases of 
gaseous ionization the characteristic distance of 
the Debye theory is so large in comparison with 
the mean distance between two ions that this 
inhomogeneity may be disregarded. 

For these reasons we shall, in all that follows, 
take p:(r) =1 and shall also replace r, by 0 which 
amounts to determining p2(r) by (1.2). 

The functions p3(r) and u(r) depend on the 
mechanism of recombination which is effective 
under given circumstances. Their determination 
constitutes the specific difficulty of the problem. 
Before starting to effectuate this determination 
for particular cases we wish to emphasize that 
these evaluations are only approximate. Whereas 


our formula (1.1) claims to be exact and to cover 


*N. E. Bradbury, J. App. Phys. 11, 267 (1940). 
7 P. Debye and E. Hiickel, Physik. Zeits. 24, 185 (1923). 
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the whole ground, the following calculations are 
only intended to give a survey of possible cases, 


II. SpectAL CASE: \ SMALL IN 
COMPARISON WITH do 


The mechanism of recombination depends 
upon the relative importance of Brownian 
movement and migration under the mutual 
attraction of the ions. Whether the one or the 
other process is prevalent can be inferred from 
the ratio of two linear quantities: 7, the mean 
distance between an ion and its nearest neighbor 
of opposite sign, and a characteristic length 
defined by 

do=e?/(KkT), (2.1) 


where ¢ is the elementary charge, K the dielectric 
coefficient and k the Boltzmann constant. If 
ado/? is large, migration is decisive, while, if 
ao/? is small, Brownian movement predominates, 
This fact has been pointed out by Loeb and 
Harper and will be shown by our treatment in 
a somewhat different way. 

A further distinction has to be introduced 
according to whether Brownian movement can 
be adequately described by the differential 
equation of diffusion or not, and this distinction 
depends upon whether the mean free path of 
the ions, \, is small or large in comparison with 
ado. 

Boltzmann’s fundamental equation must be 
considered the basis for every rigorous kinetic 
deduction. The author* has shown that the solu- 
tion of Boltzmann’s equation has to proceed 
along entirely different lines if the mean free 
path \ is small and if it is large compared with 
the linear dimensions concerned in the problem. 
The decisive length, in our case, is ad. If, then, 
\ is small compared with a» the usual treatment 
by the differential equations of transport is 
adequate, though it represents only a first 
approximation. 

From the study of the differential equation 
of relative Brownian movement Onsager® has 
deduced the following expression for the prob- 
ability p;(r): 


ps(r) =1—exp (—ao/r). (2.2) 


8G. Jaffé, Ann. d. Physik 6, 195 (1930). 
*L. Onsager, Phys. Rev. 54, 554 (1938). 
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It follows that do is, in a certain way, the 
length which decides the fate of a pair of ions. 
Though there is no definite critical distance, the 

ir will recombine if ao/r>1 and will not 
recombine if ao/r<1. 

To these two limiting cases correspond two 
entirely different mechanisms of motion. To show 
this we shall calculate the time which the one 
jon takes to travel the distance 7 relative to the 
other, upon the assumption that either migration 
in their mutual field or diffusion alone is effective. 

If k, and k_ are the mobilities of the ions and 
if tfm<?, the time required under the mutual 
attraction alone will be" 


th= K/[3e(ks +k_)]. (2.3) 


If, on the other hand, diffusion alone were 
effective the corresponding time would be™ 


tp=?/[2(D,+D-_)], (2.4) 


D, and D_ being the coefficients of diffusion. 
If tp is very much shorter than ¢, diffusion 
will be prevalent and vice versa. Hence the ratio 


tp/th=3a0/27 (by 2.17) (2.5) 


is decisive, as stated above. 
The probability function (1.2) has a strong 
maximum for 


ro=[1/(2xn) }}. (2.6) 


Hertz has pointed out that this “‘most probable 
value” is numerically very little different from 
the ‘‘mean value” 


7F=T(4/3)-(3/(4en) }}, (2.7) 


and that both are well approximated by the 
value 5/9n! which ought to replace the value 
(1/n)* usually adopted. For reasons which will 
become evident we shall use ro instead of ?. 

Now we are prepared to treat the subsidiary 
limiting cases. If, besides \<apo, the ratio ao/ro 
is large, migration will be prevalent. Then we 
have 


ux(r) =e(ky +k_)/(Kr*), (2.8) 


and for all significant values of r, namely those 
in the neighborhood of the maximum fo, (2.2) 


"In deducing (2.3) we have tacitly assumed that the 
mobilities are not affected by the very intense fields in the 
immediate neighborhood of the ions, but it must be re- 
membered that we are only making an estimate. 

"For an explanation of the numerical factor see para- 
graph containing Eq. (2.14). e 
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yields 
p3(r) =1. (2.9) 


Inserting these expressions into (1.1) and 
making use of (1.2), we find 


a, =4re(ki+k_)/K, A<do, 1o<ao, (2.10) 


which is Langevin’s formula. 
If, on the other hand, ao/ro is small diffusion 
prevails and we find from (2.2) 


bs(r) =ao/r. 


We still have to calculate u(r) and it is this 
calculation which introduces a certain ambiguity. 
Because, according to the usual treatment of 
Brownian movement, a velocity in the proper 
sense does not exist, we must define u(r) by 
some mean value. In doing so we start from the 
point-source solution for relative Brownian 
movement, i.e.,” 


1 
[4r(D,+D_)t}} 


Xexp (-_—__). (2.12) 
4(D,+D_)t 


Arguing in the usual way we might deduce from 
(2.12) the relation 


(7?)w =6(D,+D_)t, (2.13) 


but it must be remembered that this is the mean 
square of the distance traveled in time ¢. What 
is required in our case is the mean time required 


for a given displacement. This is found from 
(2.12) to be 


(2.11) 





P(r, j= 


t=r?/[2(D4+D_)]. 


By differentiating this expression we obtain 
an approximative value for u(r). As, however, 
there remains a certain ambiguity regarding 
the numerical coefficient, we shall write 


up(r) = fi(D,+D_)/r, 


(2.14) 


(2.15) 


% We have normalized (2.12) in such a way that 
JSo"4arpdr=1, i.e., 4xr*p represents the probability of 
finding two ions at time ¢ at a distance r if, at time ¢=0, 
they are both at r=0. The solution (2.12) can be verified 
either directly as the point-source solution of the differential 
equation of relative Brownian movement (see for instance 
L. Onsager, J. Chem. Phys. 2, 599 (1934)), or it can be 
derived by multiplying the point-source solutions for a 
positive and a negative ion (relative to a fixed system of 
coordinates) and subsequent suitable transformation and 
integration. 
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where f; represents a numerical factor which is 
likely to be equal to 1 and certainly is of the 
order 1. 

Introducing (2.11) and (2.15) into (1.1) we 
obtain 


ap=fi4r(D.+D_)ao, A<Kao, roao, (2.16) 


which is, but for the numerical factor, an equa- 
tion previously deduced by Harper.” 
Making use of Townsend’s relation in the form 


(Dy +D_)/(ki+k-)=kT/e, (2.17) 
the expression (2.17) can be transformed into 
ap = fra. (2.18) 


This shows the remarkable fact, previously 
observed by Harper, that the coefficient of re- 
combination is practically (exactly, for f:=1) 
the same, whether migration in the field or dif- 
fusion is alone effective, though, of course, the 
mechanisms are far different. 

We shall see in Section IV that the case where 
ro and do are comparable is of no practical sig- 
nificance. It is, however, easy to obtain an 
approximate solution which bridges the two 
limiting cases treated so far. As a first, rough, 
approximation it will be legitimate to assume 
that u(r) will then be a simple superposition of 
the expressions (2.8) and (2.15). Introducing 
this superposition and the general expression 
(2.2) into (1.1) we obtain 


ante f (1 —¢~o/r) 
0 


XK (e(ky+h-)/K+fi(D4+D_)r)po(r)dr, (2.19) 


where 2(r) is given by (1.2). As the function 
(1.2) has a strong maximum for r=ro we obtain 
a good approximation for (2.19) by using the 
saddle-value method. Remembering (2.17) we 
can write the result in the form: 


a=4n(D,+D_)ao(1—e~*)(1+fi/z), 
A<ao, (2.20) 


where we have put 
2=4o/To. (2.21) 


It will be noticed that in the limiting cases 
z>1 and z<1 formula (2.20) leads back to (2.10) 
and (2.16), but that in all other cases a is 


18, W. R. Harper, reference 4(5). 


dependent on 2, i.e., on the concentration of the 
ions. 


III. SpectaL CASE: \ LARGE IN 
COMPARISON WITH do 


With this restriction conditions are analogous 
to those prevailing in “molecular flow.” The 
solution of problems of this type can be devel- 
oped® into a series which progresses with powers 
of A“. 

We shall split the probability p3(r) into two 
factors, the probability p3'(r) that two neigh. 
boring opposite ions which are at a distance 
will be caught into related closed orbits, and the 
probability p;’’ that, subsequently, actual Tre. 
combination takes place. The probability 9,'(r) 
could be calculated rigorously by successive 
approximations and even the curvature of the 
paths could be taken into account (reference 8, 
§4). The problem offers a comparatively simple 
application of the method developed by the 
author since no boundaries have to be considered. 
Still, some law of interaction between the ions 
and the neutral molecules would have to be 
accepted and the evaluations would be very 
involved. 

For these reasons we shall adopt an approxi- 
mative procedure which is based on the analogy 
with the theory of heat radiation. In the zero 
approximation collisions are altogether dis- 
regarded and the ions are considered to be 
“emitted” from each element of volume iso- 
tropically with their mean velocity. 

We start again by considering a positive ion 
A which is supposed to be at rest and has its 
nearest negative neighbor B at a distance between 
r and r+dr. The two ions will, on the average, 
form closed orbits if B approaches A within a 
distance less than‘ 


d=2e?/(3KkT) =2a)/3. (3.1) 


We are only interested in the case d<r. Then 
we shall have in our present approximation 
ps’ (r) = (1/40) (@x/r’)=a?/4r? (3.2) 
and 
u(r) =C\2, (3.3) 


where ¢i2 is the average relative velocity. 
Assuming provisionally p3'(r)=p3(r) and in- 
serting (3.2) and (3.3) into (1.1) we find as zero 
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approximation 
a = wdc. (3.4) 


This is, of course, the usual formula of the 
kinetic theory (for a cross section of radius d), 
and it becomes identical with the limiting value 
of Thomson's theory if we take 


C12 = (642+ ¢_*)}. (3.5) 


If, now, collisions occur, each of them is con- 
sidered an act of “‘absorption”’ and ‘‘re-emission,”’ 
the re-emission being equally probable in all 
directions and the velocity remaining, on the 
average, unaltered. Since the probability of a 
collision along an element of path ds is given by 
ds/d, the reciprocal of the mean free path must 
be considered the “‘coefficient of absorption.” 

Let us again assume that a negative ion B is 
emitted at a distance r from its nearest positive 
neighbor A at an angle 3; with the line joining 
them. We obtain a first ‘approximation by as- 
suming that one, and only one, collision occurs 
outside the sphere of radius d around A. If the 
collision occurs at a distance 7; from the starting 
point of B, at a point which is at a distance pi 
from A, B will have to travel a distance r;+ p,—d 
~ri+pi before being caught, and its average 
radial velocity (as required for (1.1)) will be 
u(r) =Ciar/(r1+ pi). 

It simplifies matters to calculate in one step 
(u(r) ps'(r))w, ie., the average of the product for 
all possible loci of the collision. We obtain" 


(u(r) ps'(r))mw 


C12? 4 en dr, Pil dp, 
‘ f connie f min, ie 
8r wo n plo pi(fitpi) 


4 It must be pointed out that, in establishing (3.6), we 
have made use of (3.2) (with p: replacing r) for the proba- 
bility of capture after the collision. However, this is not 
quite correct because it is not true that d<p, in the whole 

ce of integration. The correct value of the solid angle 
iooniaden capture would be 2=2x[1—(1 —d*/p,*)*] 
rather than wd*/p,*. We have carried out the calculations 
with this value and verified that the result is not noticeably 
different from (3.6). 

In determining the region of integration we have omitted 
the whole solid angle 2’ subtended at B by the sphere of 
radius d. This means omitting the negligible probability 
of cases in which the ion remains inside 2’ in spite of a 
collision. The space inside 2’ behind the sphere has to be 
omitted. 

Finally it must be remarked that (3.6) is not a sym- 
metrical function of the ions A and B unless } is taken as 
the average value for positive and negative ions. In a 
somewhat more rigorous deduction (3.6) would have to be 
replaced by its arithmetic mean for A, and d_. 
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where 


Pi0= [(r’— r1)?+d? }i, 
r’ =(r?—d*), (3.7) 


pu=rt+ni, 


The first integration in (3.6) is immediate, the 
second can be carried out to a good approxima- 
tion since the function r;? exp (—r./A) has a pro- 
nounced maximum for r;=\. We thus find 








(u(r) ps'(r))m 
Ci2d? r+h [(r’—d)?+d?}! 
= {tog —log 
8x21 r+2r [(r’—d)2+-d27}+2) (3.8) 





C120? dd 
= o(-. -). 
8? Ar 


In order to calculate a, we must remember that 
it consists now of two parts, the one being due to 
ions which are caught without having undergone 
a collision, and the second to ions which have 
made one collision. The first part we obtain by 
multiplying (3.4) by exp (—7ro/A), the second by 
inserting (3.8) and (1.2) into (1.1) and evaluating 
by the saddle-value method. Thus we finally 
obtain 


asardenlen4 te (S, 2) 
2? d’ To 


dd 
=nd'ash(-, “), (3.9) 
aN To 


The function ¥, approaches unity for small 
and large values of A/ro. Neglecting terms of the 
order (d/ro)? we find 


1 To ’ 
a, =1d%nfe-~-+-(—) ae eee ; 
2\A 


A/ro>i, (3.10) 


and 


r 
ani = rd*eul a A/roK1. (3.11) 
To 


For the latter case, however, our present ap- 
proximation is not adequate. The author, there- 
fore, has carried through the next approximation 
taking account of those ions B which reach A 
after two collisions. The result is too lengthy to 
be quoted in extenso. In the limiting case \/ro<1 
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it is found 


2r A\? 
an’? = rd°eua{ 1-—+6(—) +- nd , 
To To 


d/roKt. (3.12) 


The comparison with (3.11) shows that the 
convergence of our procedure for this limiting 
case may be slow, though the essential feature, 
convergence towards 2d’cy. for d/r—-0, is 
preserved. 

So far we have identified ;'(r) with p3(r), i.e., 
we have assumed that each pair of ions which is 
caught into closed orbits will ultimately recom- 
bine. It is, however, evident from J. J. Thomson’s 
well-known argument that this need not be the 
case. Therefore the a so far calculated has still 
to be multiplied by the probability p;”’ as defined 
above. 

A rigorous calculation of this p3’’ would have 
to be based on a quantum-mechanical study of 
the interaction between ions and neutral mole- 
cules on the one hand and ions of different signs 
on the other. Such a study does not seem feasible 
at the present time. We may, however, accept 
Thomson’s treatment as a first approximation, 
though it cannot hold down to zero pressure, as 
was pointed out by Loeb." 

The probability ;’’ as calculated by Thomson 
is given with sufficient accuracy by *° 


bs!’ = g2(d/d) = 2w—w?, (3.13) 


where 


d? 2d 
o-(1 -— | 1 ~en(+1) | (3.14) 
2d? d 


Hence our final expression for a, becomes 
oO = wd?cy291(d/X, d/ro) g2(d/d), Ad. (3.15) 


This result, though intimately connected with 
Thomson’s, is distinguished from it by two essen- 
tial features: It depends upon d/ro, i.e., upon the 
concentration of the ions, and it gives a marked 
maximum for a at 7>=A, as can be shown from 
the form of ¢:(d/d, d/ro). Sufficiently above 
pressures where 7o~X, and for ionizations suffi- 


15 See reference 1, p. 130. 


1% See reference 1, p. 119. The function g is found 
tabulated there. 


ciently weak to make d/ro a small number, y, 
can be approximated by 1 (see (3.11) and (3.12)) 
and our formula reduces to Thomson’s 


ay = 9d7C12¢2(d/X), A>d, d<Xro. (3. 16) 


IV. DiscussION AND COMPARISON WITH 
EXPERIMENT 


We have seen that our general expression (1.1) 
reduces, for well-defined special conditions, to 
the formulae of Langevin, Harper and Thomson. 
However, it must be emphasized that only under 
the restrictions specified by (2.10), (2.16), and 
(3.16) does the coefficient of recombination be- 
come independent of the ionic density, m. The 
transitional formulae, (2.20) and (3.15), involve 
the parameter 


g=4/ro=2ao/3ro=d(2rn)!. (4.1) 


It seems to be generally accepted that the 
Langevin formula holds with adequate accuracy 
for sufficiently high pressures. As pointed out by 
Harper, this does not prove the correctness of 
the mechanism by which Langevin deduced his 
formula, as shown by the formal similarity be- 
tween the expressions for a, and ap (see (2.18)). 

It is pretty certain that in no direct experi- 
ments on volume recombination is the migration 
of the ions under their mutual attraction the 
predominant feature. The highest observed 
values of m are of the order 10’ ions per cc. Even 
in liquid dielectrics (which behave in this regard 
like very dense gases) m never rises above 2X 10°. 
The corresponding values of ro are 2.5 10-8 and 
0.9210-*. As d=4.05 X10-°273/T, q is never 
larger than about 4X10-* and in most cases is 
considerably smaller. Under these circumstances 
diffusion prevails and the agreement of experi- 
mental data with Langevin’s formula really 
proves the diffusion formula (2.16). 

There is, however, indirect evidence for the 
correctness of Langevin’s formula, including its 
mechanism. The theory of columnar ionization, 
as developed by the author, seems to be generally 
accepted for the ionization with heavy particles. 
The fundamental concept of this theory is that, 
inside the columns, recombination (and diffusion) 
take place according to the usual laws and with 
the same coefficients as in cases of more homo- 
geneous ionization. 
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‘ 
d/A 
Fic. 1. The diagram represents a/az as a function of 
d/s. The curves for a),/az are calculated from (3.16) and 
(3.12) for various values of g, and the curve for ap/az is 
calculated from (4.3) with fi=f,=1. The observational 
points are reduced values. 


The initial concentration in the axis of a 
column” is of the order of 510° in air of atmos- 
pheric pressure and rises in liquid hexane to 
2.50108, giving values of g from 1.210 to 
9.5. This brings columnar ionization in dense 
gases within the scope of formula (2.20)'* and for 
liquid dielectrics even of (2.10). Thus the agree- 
ment of the columnar theory with experimental 
facts indirectly confirms Langevin’s theory. 

At low pressures, say up to 1 atmos., Thomson’s 
theory seems to give an adequate representation 
of experimental facts. That is to be expected 
from our theory as long as \ is somewhat larger 
than d and, at the same time the pressure is 
sufficiently high to make \ smaller than ro. 

In the region where \} becomes comparable 
with d neither the formula for ap nor that for a 
is a sufficiently good approximation for our 
general formula (1.1), as they are only first steps 
in the evaluation for \/d<1 and A/d>1, re- 
spectively. Theoretically, it would be possible 
to carry the two approximations on until they 
meet, but practically, that is out of the question. 
The failure is, however, not due to a breakdown 
of the fundamental concepts (as in the case of 
Langevin’s or Thomson’s theory beyond their 
domains of applicability) but to our inability to 
evaluate the general formula (1.1). Conditions 


7G, Jaffé, Ann. d. Physik 32, 152 (1910); H. J. van der 
Byl, Ann. d. Physik 39, 170 (1912). 

% According to (2.20) the coeff. of rec. does not vary 
much with g. Assuming f:=1, a has a maximum, for 
4o/ro~2, which is about 30 percent higher than the values 
at the extremes. 


are thus exactly the same as in the transition 
from Poiseuille’s ‘‘frictional flow’’ to Knudsen’'s 
“molecular flow.” 

As is well known, the critical region is in the 
neighborhood of atmospheric pressure, since 
A\=1X10~- (760/p)(A/Aa) where \/Xz is the ratio 
of the mean free path of the ion to that of a 
molecule in air at NTP.'® For this region of 
transition it must be tested whether the two 
approximations approach each other sufficiently 
to make a graphical interpolation possible. 

According to (3.16) a, is given as a function of 
d/X, whereas ap is represented by (2.16) as a 
function of (D,+D_). To reduce the two repre- 
sentations to common terms we must represent 
D, and D_ by atomic data and this reduction 
introduces a further ambiguity. For our purpose 
it will suffice to introduce a mean coefficient of 
diffusion, defined by D=(D,+D_)/2, a mean 
velocity ¢=(@,+2)/2, and write 


D=f.tr/3, (4.2) 


where A is also the average value for positive 
and negative ions and f: a numerical factor. 
According to the usual elementary theory fs is 
unity, but we insert the factor f. of order 1 to 
leave room for the improvements of the more 
elaborate theories. 

Making use of (4.2) and (3.1) and remembering 
that c:2.=2'2, we find 


ap = md*cy9°fifo2 /d. (4.3) 


If we introduce the limiting value a, =1d’ci_. we 
see that both ap/az and a/az, the latter given 
by (3.16), are universal functions of d/A. We have 
represented these two functions by the heavy 
lines of Fig. 1, choosing f;=f.=1. It will be seen 
that a fair representation of a/az for all values 
of d/d can be gained by the two approximations 
from either side, though, of course, we have no 
quantitative criterion to show at which particular 
values of d/d the deviations become noticeable. 

As to the agreement with empirical data we are 
of the opinion that neither the present state of 
experimental research nor our evaluations of 
Sections II and III would justify a comparison 
of absolute values. Therefore we have adopted a 
procedure of testing relative values. As a/az, 


19 L. B. Loeb, reference 1, p. 118. 
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should be a universal function of d/A, we have 
tried to reduce the a-scale, as well as the p scale, 
in such a way that the experimental values fit the 
theoretical curves. This, of course, amounts to 
taking ¢i2 (or the mass of the ion) and \ as two 
constants which are to be determined experi- 
mentally. 

All measurements of a in air in the interval 
between 1 and 8 atmos. have been included. 
These are the measurements of Langevin, 
Machler, Kraus and Sayers.?° The values adopted 
for az-10°* are, respectively, 1.83, 1.71, 1.70 and 
2.35 and the values for \/A._ 0.40, 0.41, 0.33, 0.26, 
respectively.”4 

It will be seen that the points thus determined 
seem to confirm the theoretical curve, but we do 
not want to stress this kind of agreement. Meas- 
urements in the indicated interval which include 
the determination of a and the mobilities in the 
same sample of gas are an urgent desideratum. 
The differences among the values of a, and of 
d/Xa, as assumed above, are probably due to the 
fact that different observers have been dealing 
with ions of different age and nature. 

So far we have limited our discussion to the 
case where the parameter g is very small, which 
reduces the more general expressions (2.20) and 
(3.15) to the simple forms (2.16) and (3.16). To 
illustrate the influence of g we have calculated 
a,/a, from (3.12) for the values g=0.01, 0.025 
and 0.05 and inserted the curves into the 
diagram. 

It will be seen that a/az decreases with in- 
creasing g and that means with increasing n.” A 
dependence of a on the concentration of the ions 
signifies that the simple law dn/dt = — an? ceases 
to hold, but the change of a predicted by theory 
is so slow that it seems adequate to retain the 
formal law and to speak of a variation of a with 
concentration. This procedure has always been 
adopted by experimenters. 


20 P. Langevin, Ann. Chim. Phys. 28, 289, 433 (1903); 
W. Machler, Zeits. f. Physik 104, 1 (1936); P. Kraus, 
Ann. d. Physik 29, 449 (1937); J. Sayers, Proc. Roy. Soc. 
A169, 83 (1938). 

#1 All measurements refer to ‘“‘room-temperature,” but 
the exact values are not given in some cases. As d depends 
on the temperature this ought to be known for a more 
exact evaluation. We have chosen ¢=18°. Furthermore, 
Sayers’ values had to be taken from his diagram. _ 

# A change for ap is hardly to be expected for realizable 
values of m, as long as volume ionization only is considered. 


It is well known that the coefficient of recom- 
bination, as determined experimentally, is not 
constant. Even after the stage of initial recom- 
bination is over, a goes on changing for a con- 
siderable time. This lack of constancy has been 
explained** by a gradual change in the nature of 
the ions and there can hardly be a doubt that this 
is the correct interpretation. The aging of ions 
will generally tend to decrease a with time and, 
therefore, with decreasing ionic density, as is 
also confirmed by the corresponding effect on the 
mobilities. Apart from this aging effect, however, 
all recent observers** have noticed that a is larger 
for lower values of the ionic concentration. The 
author is of the opinion that this dependence, 
which coincides in direction with that predicted 
by the present theory, is a genuine effect not due 
to secondary causes. The few quantitative data 
which are available indicate that the observed 
effect is considerably larger than the effect fore- 
seen by our theory. This, however, is only true 
as long as the ionization is supposed to be homo- 
geneous in all cases. The effect would be much 
stronger if, for the lower ionic densities, the initial 
inhomogeneity became more. marked, an effect 
for which there are distinct indications.** The 
dependence of a on ionic density is the second 
point where an extension of experimental evi- 
dence is desirable. 

We do not want to dwell here on the question 
of recombination at very low pressures. We have 
observed (p. 974), however, that, under conditions 
where \~7, something like selective recombina- 
tion should occur according to Eq. (3.15). This 
fact shows how careful one ought to be in extend- 
ing Thomson’s formula to very low pressures, 
whether in the atmosphere or the discharge tube. 
This point requires further study, but the study 
will, no doubt, be complicated by the fact that 
in the same interval of pressures the carriers of 
electricity change character. 

The author wishes to express his gratitude to 
the Department of Physics of the University of 
Colorado for courtesies extended. 


%1L.C. Marshall, Phys. Rev. 34, 618 (1929); O. Luhr, 
Phys. Rev. 35, 1394 (1930); M. E. Gardner, Phys. Rev. 
53, 75 (1938); J. Sayers, reference 20. 

* A. D. Power, J. Frank. Inst. 196, 327 (1923). 

% The maximum indicated by 3.15 is mot erased by the 
sp ee, but a second minor maximum is added 
or fo= 2d. 
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Low Temperatures 
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The magnetic susceptibility of the anhydrous chlorides of the iron group have been studied as 
a function of field strength up to 32,000 gauss and as a function of temperature down to 13.9°K. 
Those compounds whose room temperature susceptibility characteristic indicates a positive 
Curie temperature, also show a field strength dependence of the susceptibility at 13.9°K and 
20.4°K. Chromic chloride has a very steep magnetization curve at these low temperatures, with 
a large initial susceptibility increasing with decreasing temperature. Ferrous, cobaltous, and 
nickelous chlorides have an S-type magnetization curve characterized by a small initial sus- 
ceptibility which decreases with decreasing temperature, followed by a very large susceptibility 
and eventual approach to saturation, both increasing with decreasing temperature. The 
susceptibility of many of the compounds reaches a maximum at temperatures which do not 
appear to be correlated with the presence of field dependence, but these temperatures do 
correspond to the temperatures of known specific heat anomalies. 


INTRODUCTION 


HE behavior of the magnetic susceptibility 

of most paramagnetic substances at ordi- 
nary temperatures is adequately described by the 
Curie-Weiss law, x(7—©)=C, where x is the 
susceptibility, T is the absolute temperature, 0 
is the Curie temperature, and C is a constant. 
Theoretical considerations indicate that the value 
of 0 is determined by the nature and magnitude 
of the interaction within the magnetic atoms 
(spin-orbit coupling) and between the magnetic 
atoms and their surrounding atoms. This latter 
interaction may be either between the magnetic 
atoms themselves, both dipole-dipole and ex- 
change effects, or between the magnetic atoms 
and the crystalline (electric) field produced by 
the neighboring ions in the substance. When all 
interactions become negligible, as in a gas of free 
carriers of magnetic moment, © reduces to zero, 
and the susceptibility follows Curie’s law, xT =C. 
A magnetic carrier may be considered as free 
when the energy required to rotate it against 
interaction forces is small compared to kT. Thus 
the deviations from Curie’s law are sufficiently 
enhanced at low temperatures to be a valuable 
source of information concerning the interaction 
effects. 


* Contribution No. 457 from the Research Laboratory of 
Physical Chemistry. 


The anhydrous chlorides of the elements of the 
first transition series, the iron group, are of 
particular interest because the magnetic atoms 
are sufficiently concentrated to produce large 
interaction effects. In addition, the crystal struc- 
tures of this group of compounds appears to be 
isomorphic, although crystal structure studies 
have not been made for all the compounds. This 
facilitates comparison between the atoms of 
varying magnetic moment. Finally, there seems 
to be no indication of any metallic or semicon- 
ducting behavior in these compounds, so the 
complication that electrical conductivity intro- 
duces into magnetic properties does not» have to 
be considered. 

The cryomagnetic anomalies associated with 
some of these compounds have been studied at 
the Kamerlingh Onnes Laboratory of the Uni- 
versity of Leiden, with many interesting results. 
The present investigation is an elaboration of 
this work. 


METHOD OF MEASUREMENT 


The measurements were made with a Bitter 
type high field solenoid (magnet No. 1 of the 
M. I. T. magnet laboratory)! which permitted 
investigation up to 32,000 gauss. A nonhomo- 


1F, Bitter, Rev. Sci. Inst. 10, 373_(1939). 
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geneous field method (Curie) described by Bitter 
and Kaufmann? was used for the determination 
of the susceptibilities. The specimens consisted of 
the powdered compound sealed off in a thin- 
walled Pyrex ampoule. The lower part of the 
ampoule was shaped as a cylinder, 8 mm O.D. 
and 1 cm high, and this joined to a long tubular 
neck, 3 mm O.D. and 5 cm long. The cylindrical 
part of the ampoule was filled with the powder 
under inevstigation and a short glass-wool plug 
was placed in the neck to keep the powder in 
place. The ampoule was then filled with helium 
gas and sealed off at the top. 

The experiments were carried on in a long 
Dewar flask (2 cm I.D.) which fitted in the 
magnet. A cylindrical copper working chamber, 
connected to and supported by a long, 9 mm 
supernickel tube, was placed within the Dewar at 
a suitable region of the magnetic field. The 
specimen tube hung within the copper chamber 
and was supported by a gold chain connected to 
one arm of an analytical balance. This arrange- 
ment permitted independent control of the 
atmosphere surrounding the specimen tube and 
of the refrigerant contained in the Dewar flask. 
The following fixed temperatures were used in the 
measurements: solid CO: (193°K), boiling Ne 
(77.3°K), freezing Ne (63.1°K), boiling He 
(20.4°K), and freezing He (13.9°K). For the 
liquid N2 temperatures and above, He gas was 
slowly passed out of the neck of thé supernickel 
tube to prevent condensation of water and air on 
the specimen tube. For the temperatures made 
available by liquid He, helium gas was used in the 
working chamber. 


PREPARATION OF SPECIMENS?’ 


As a result of the extremely hygroscopic nature 
of the anhydrous chlorides of the iron group, all 
preparations and transfers were made in glass 
vessels, in vacuum or dry inert gases, and the 
final specimen tubes were sealed off with dry 
helium. The helium was used to insure rapid 
temperature equilibrium throughout the speci- 
men at all the temperatures used. 


(1939) Bitter and A. R. Kaufmann, Phys. Rev. 56, 1044 

’ The compounds were prepared by Professor R. C. 
Young of the M.I.T. Chemistry Department in collabora- 
tion with the first author. 
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NiCl,, CoCle, FeCl, MnCle—These com- 
pounds wére all made from the C.P. hydrates, 
The hydrates were dehydrated in vacuum at 100 
to 200°C. The resultant product, containing an 
appreciable quantity of oxychloride, was then 
heated in vacuum to about 650°C. The pure 
chlorides sublime and condense in characteristic 
crystals on the cold portions of the tube, the 
oxide impurities remaining behind. The sublimed 
crystals have a very high purity, and in the case 
of ferrous chloride, less than 0.1 percent of ferric 
ion was present. 

FeCl;, CrCl3.—These were made from com- 
mercially available C.P. anhydrous sublimed 
material. They were resublimed in vacuum, the 
FeCl; at 250°C, the CrCl; at 650°C. 

VCl;.—Made by decomposition of VCl, at 
160°C in dry COs. The VCl, was produced by 
reacting Cl. with V metal at 500°C. ; 

VCly.—Made by reduction of VCl; at 500°C 
with Ho. 

TiCls.—Made by reduction of TiCl,, com- 
mercially available, at 650°C with He and con- 
densing the TiCl; on a cold surface. 

TiCle—Made by vacuum decomposition of 
TiCl; at 475°C. 

CuCle.—Made by reacting Cl, with Cu at 
400°C. 


RESULTS 


The results of the measurements are shown in 
Figs. 1 and 2 and in Tables I and II, where x is 
the susceptibility per gram mole and a is the 
magnetization per gram mole in c.g.s. units. // is 
the field strength in gauss, and T is the absolute 
temperature in °K, unless otherwise specified. 
Figure 1 shows the reciprocal susceptibility as a 
function of temperature. Figure 2 shows the 
magnetization as a function of field strength for 
those compounds which exhibited a dependence 
of the susceptibility upon the field strength. For 
these latter compounds, the initial differential 
susceptibility, x =do/dH, is used in the reciprocal 
susceptibility plots. 

The susceptibility constants given in Table I 
have all been corrected for the diamagnetism of 
the cation and anion. The susceptibility values 
used in the graphs are the total measured values, 
except in the case of the titanium compounds 
where the corrected values are used. The values 
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Fic. 1. The recip- 
rocal susceptibility 
per gram mole as a 
function of the ab- 
solute temperature 
for FeCl;, FeCl, 
MnCh, CoChs, 





CrCls, CrCl, NiCle, 








VCl, 


VCls 














VCl;, VCle, CuCl, 
TiCl; and TiCle. 
The correction for 
+ thediamagnetism of 
the cation and 
anions in TiCl; and 
“| TiCl, has been ap- 
plied. 








given in Table I for the diamagnetic susceptibility 
per gram mole, xa, have been calculated for the 
complete compound (cation plus anion) by the 
method of Angus‘ for modifying Slater’s® formula 
for ionic diamagnetic susceptibilities. This modi- 
fication is used since the empirical agreement is 
improved, although the theoretical basis of 
Slater’s method is better. 

The value of the Curie temperature* was 


*W. R. Angus, Proc. Roy. Soc. A136, 569 (1932). 

5]. C. Slater, Phys. Rev. 36, 57 (1930). 

*In this paper the term “‘Curie temperature” refers to 
the value of © found from the extrapolation of measure- 
ments made at comparatively high temperatures. The 
term “Curie point” refers to the actual temperature at 
which the magnetic behavior sharply changes in character, 








determined from the slope of the reciprocal 
susceptibility vs. temperature curve at room 
temperature, corrected for the diamagnetic sus- 
ceptibility. The effective number of Bohr mag- 
netons, p, was calculated by means of the formula 
p=2.84[x(T—0) ]!. The temperature of the 
maximum susceptibility given in Table I corre- 
sponds to the temperature of the reciprocal 
susceptibility minima (where such exist) in Fig. 1. 

Specific heat studies of the nickel, cobalt, 
ferrous, manganous and chromic chlorides have 


as determined by actual measurements in this tempera- 
ture range. This distinction is similar to the usual distinc- 
tion between the paramagnetic and ferromagnetic Curie 
point. 
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Fic. 2. The mag- 
netization per gram 
mole as a function of 
the field strength at 
20.4°K and 13.9°K 
of CrCls, FeCl, 
NiCl,, and CoCh, 
In the case of NiCl, 
the dotted line re 








resents the straight 





line continuation of 
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been made by Trapeznikowa and Schubnikow,* 
and their results, tabulated in Table I, show that 
the anomalies in the specific heats occur at 
temperatures which correspond closely to the 
temperatures of maximum susceptibility. 

In Table II, the magnetization data are listed 
for the compounds, which showed a field de- 
pendence of the susceptibility at Hz tempera- 
tures. The differential susceptibility (slope of the 
magnetization curve) is given for the temperature 
of boiling Hz, 20.4°, and freezing He, 13.9°, and 

60. Trapeznikowa and L. Schubnikow, Phys. Zeits. 
Sowjetunion 7, 66 and 255 (1935); O. Trapeznikowa, L. 
Schubnikow and G. Miljutin, Phys. Zeits. Sowjetunion 


9, 237 (1936); O. Trapeznikowa and G. Miljutin, Phys. 
Zeits. Sowjetunion 11, 55 (1937). 


both the initial and maximum values are pre- 
sented. Ascan be seen in Fig. 2, the magnetization 
curves for 20.4° and 13.9° cross each other in the 
case of FeCl, and CoCle. The magnetization and 
field strength of this point of intersection is shown 
in Table II. 

The precision of the results was not the same 
for all the specimens, and depended chiefly upon 
the magnitude of the forces involved. The 
analytical balance could detect changes of 0.1 
milligram force, but this sensitivity could only be 
used with great difficulty because of the dis- 
turbance created by the gas escaping from the 
specimen chamber. This reduced the working 
sensitivity of the balance to about 1 milligram. 
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The forces were usually between 1 gram and 100 
grams for the highest field strengths. Generally, 
the accidental error in any measurement was 
about 0.1 percent. The sources of systematic 
errors were as follows: (1) calibration of the field 
and field gradient, (2) change of position of the 
specimen in the field caused by the stretching 
of the gold chain or thermal contractions, (3) 
errors in the weight of the compound in the 
specimen tube, and (4) error in the diamagnetic 
susceptibility correction for the Pyrex specimen 
tube and gold chain. These errors were probably 
of the order of magnitude of 1 percent. For 
specimens of very small susceptibility, such as the 
titanium compounds, the forces are very small 
and the correction (4) is relatively large, so that 
the systematic error may be 10 percent or more. 

Another source of error was found in com- 
pounds with very large susceptibilities at liquid 
H, temperatures, such as FeCl». In these cases 
the forces on the individual crystals are sufficiently 
large to pack the crystals together. This causes 
the density of the specimen to be increased 
slightly, when the magnetic field is applied, in a 
manner not completely reversible. This packing 
effect results in an apparently increased mag- 
netization of the specimen when the field is being 
decreased from its maximum value, as compared 
to the magnetizations at the same fields when the 
field is being increased to its maximum value. 
After repeated magnetizations at the highest 
fields, the measurements follow a reproducible 
hysteresis loop in which the difference between 
increasing and decreasing field points is about 1 
percent. This difference could be greatly in- 


TABLE I. Susceptibility data for the anhydrous chlorides of 
the elements of the first transition series. 








Temp. Temp. 





OF Spec. 
—Xxd Max. Hr. 
Compound X10* x(T—®) ? e x ANOMALY 
TiCl; 80.4 0.25 1.42 
TiCl, 66.7 0.521 2.05 —702 
VCl; 82.3 1.005 2.85 — 30.1 30 
VCl, 66.9 2.13 4.15 —565 
CrCl; 82.9 1.89 3.90 27.0 ~20 16.8 
CrCl, 66.1 3.26 5.13 —149 ~40 
MnCl, 65.0 4.07 5.73 — 3.3 None 
FeCl; 82.5 4.07 5.73 — 11.5 
FeCl, 63.9 3.59 5.38 48.0 24 23.5 
CoClh, 62.8 3.46 5.29 38.1 25 24.9 
NiCl, 61.8 1.36 3.32 68.2 ~50 49.6 
57.0 
8 0.536 2.08 -—109 ~70 


CuCl, 60. 








creased by shaking the specimen tube and 
loosening the cystals. Repeated magnetizations 
re-established the small hysteresis loop. For this 
reason, the curves in Fig. 2 have been plotted 
from measurements taken in the direction of 
decreasing field strength after stabilization by 
repeated magnetizations. 

It is of interest that most of these compounds 
form layer lattice crystals, so that the individual 
crystal particles take the form of flat flakes. When 


TABLE II. Magnetization data for the compounds whose 
susceptibility is a function of magnetic field strength. 








DIFFERENTIAL SUSCEPTIBILITY INTERSECTION 
Com- INITIAL MAXIMUM PoINT 
POUND 20.4° 13.9° 20.4° 13.9° o H 





CrCl; | 0.947 6.75 
FeCl, | 0.248 0.102 1.88 4.27 11,700 13,500 
CoCl, | 0.294 0.272 0.422 0.456 5900 14,800 
NiCl, | 0.0700 0.0700 0.0961 0.0981 














packed into the specimen ampoule, the flakes 
tend to lie parallel to the base of the ampoule 
(i.e., at right angles to the direction of the field) 
Since the individual crystals are undoubtedly 
anisotropic, this results in a small anisotropy of 
the specimen. Serres’ has found that this speci- 
men anisotropy is about 1 percent at room 
temperature for CrCl; crystals. This type of 
specimen anisotropy was probably present in all 
the compounds studied. 

The correction for the demagnetizing field has 
not been made, since it is not certain whether the 
specimen is behaving as a whole or each indi- 
vidual flat crystal is demagnetizing itself. The 
true magnetic field inside the crystals is H;=H, 
—DI, where H, is the applied field, D is the 
demagnetizing coefficient, and J is the mag- 
netization per unit volume. This leads to the 
formula 1/x:=1/xa—Dp/M, where x; and x. are 
the gram mole differential susceptibility (do /dH) 
true and apparent, respectively, p is the density, 
and M is the molecular weight of the compound. 
The maximum correction would apply when the 
flat crystals are demagnetizing themselves indi- 
vidually ; D=4n and p is the crystal density, and 
Dp/M=0.3 to 0.2 approximately. The minimum 
correction would apply when the specimen 
ampoule is demagnetizing as a whole; D=4 and p 
is the specimen density (approximately } the 


7 A. Serres, Ann. de physique 17, 5 (1932). 
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crystal density) and Dp/M=0.05 to 0.03, which 
is about $ as much as in the maximum case. The 
demagnetizing correction is largest for the initial 
susceptibility of CrCls at 13.9°, x,=6.75. The 
maximum assumption would be more than suffi- 
cient to make this initial susceptibility infinite. 
The minimum assumption would make the true 
initial susceptibility about 9.0. At 13.9°, FeCl. has 
a maximum differential susceptibility, x.=4.27. 
The maximum demagnetization assumption 
would make this differential susceptibility infi- 
nite. The minimum assumption would raise it to 
x:=5.1. These demagnetization corrections are 
so very large as compared with the usual cor- 
rection for paramagnetics because the above 
molar susceptibility for CrCl; is about 50 times as 
large as the value for dilute chromic compounds 
(spin-only value), and for FeCl, is about 20 times 
as large as the value for dilute ferrous compounds. 

No remanences were observed in any of the 
specimens, within the sensitivity of the appara- 
tus. The magnetizations were measured with the 
field first in one direction and then reversed. Any 
remanence would appear as a difference between 
the two sets of measurements. No such difference 
was observed. The recent Leiden measurements,® 
however, indicate that some of these compounds 
exhibit remanence. Inasmuch as the apparatus 
used for the Leiden work had a very high sensi- 
tivity their results may not be at variance with 
the present work. 


DISCUSSION AND COMPARISON OF RESULTS 


CuCle.—The specimen used contained suff- 
cient impurity in the form of CuCl to invalidate 
the absolute values of the measured susceptibility. 
The relative values agreed very well with the 
Leiden measurements,? so the absolute values of 
the results were adjusted to agree with the 
Leiden susceptibility at room temperatures. The 
CuCl contributes only a diamagnetic correction. 
No field dependence of the susceptibility was 
found. The adjusted results are shown in Fig. 1. 

NiCl2.—The room temperature value of the 
susceptibility agrees very closely with that found 


’W. J. de Haas and B. H. Schultz, Physica 6, 481 
(1939); W. J. de Haas, B. H. Schultz and J. Koolhaas, 
Physica 7, 57 (1940). 

(1934) J. de Haas and C. J. Gorter, Leiden Comm. 215a 





by Serres'® who reports a Curie temperature of 
95°K. The field dependence at hydrogen temper- 
atures agrees qualitatively with that found at 
Leiden."' The Leiden measurements show no 
dependence upon the temperature from H, 
temperatures to He temperatures, but their 
investigations went only to 20,000 gauss. 

CoCl..—The room temperature value of the 
susceptibility agrees reasonably well with that of 
Fehrenbach” who found a Curie temperature at 
26.6°. The recent measurements of deHaas and 
Schultz® do not agree with the older Leiden 
measurements," but are in good qualitative 
agreement with the present work. They find a 
Curie temperature of 19°. They also report that 
the magnetization curves at boiling and freezing 
He cross at a field of 6500 gauss, as compared 
with the 14,800 gauss in the present work. Their 
ratio of initial susceptibilities is about the same, 
although their values are lower. 

FeCl;.—The Leiden measurements™ on this 
compound are very unsatisfactory as they used 
an incomplete Gouy method, and the results do 
not permit simple interpretation in terms of 
magnetization. They found indications of a 
maximum in the susceptibility at Hz tempera- 
tures, the maximum decreasing and shifting to 
higher field strengths with decreasing tempera- 
ture. The present results are at variance with this. 

FeCl;.—The room temperature results of 
Lallemand™ agree very well with the present 
work. No field dependence was found at any 
temperature. The susceptibility deviates slightly 
at low temperatures from the Curie-Weiss law. 

MnCl,.—The present work agrees qualitatively 
with the Leiden measurements” on this com- 
pound, and is in good quantitative agreement 
with the room temperature results of Lallemand™ 
and Fehrenbach.” At the temperature of freezing 
Hz (13.9°) the expected beginning of paramag- 
netic saturation could be observed. No other field 
dependence of susceptibility was found. 


10 A. Serres, Ann. de physique 20, 441 (1933). 

1H. R. Woltjer, Leiden Comm. 173b (1925); H. R. 
aa and H. Kamerlingh Onnes, Leiden Comm. 173c 
(1925). 

12 C. Fehrenbach, J. de phys. et rad. (7) 8, 11 (1937). 

13H. R. Woltjer and E. C. Wiersma, Leiden Comm. 
201a (1929). 

14 A, Lallemand, Ann. de physique 3, 97 (1935). 

%H. Kamerlingh Onnes and E. Oosterhuis, Leiden 
Comm. 129b (1912). ' 

















CrCle.—The work of Serres’ at room tempera- 
ture agrees very well with the present work. 
Serres reports a Curie temperature of —115.5°. 
Lips'* studied this compound down to 55°K, and 
his values, which indicate a Curie temperature of 
—128°, are in excellent agreement with the 
present ones. No field dependence of the sus- 
ceptibility was found. 

CrCl3.—The Leiden results" are in good agree- 
ment with the present work on this compound. It 
is of interest that this compound was the only one 
that exhibited simple, steep magnetization curves. 
In this connection it should be noted that the 
Curie point indicated by the specific heat 
measurements is 16.8°, although the field 
dependence of the susceptibility occurs above this 
temperature as well as below. 

VClz, VCl3.—The previous work of Klemm 
and Hoschek"’ down to 90°K is in good quanti- 
tative agreement with the present results. No 
field dependence of the susceptibility was found 
at any temperature. 

TiCle, TiCl;.—These compounds have not 
been investigated before, probably because of the 
chemical difficulties involved in their production. 
The dependence of their susceptibility upon 
temperature is completely different from that of 
the other compounds. No field dependence of the 
susceptibility was found. The values shown in 
Table I for TiCl; were calculated from the room 
temperature value of the susceptibility, with 0 
assumed equal to zero. 


SUMMARY 


The salient features of the experimental results 
may be summarized as follows: 
% E. Lips, Helv. Phys. Acta 7, 537 (1934). 


@W. Klemm and E. Hoschek, Zeits. f. anorg. allgem. 
Chemie 226, 359 (1936). 
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(1) None of the compounds showed any de- 
pendence of the susceptibility upon field strength 
at temperatures above 63°K. 

(2) At He temperatures (13.9° and 20.4°) a 
dependence of the susceptibility upon field 
strength was found for NiCl., CoCls, FeCle, and 
CrCl;. These compounds are the only ones which 
exhibit positive Curie temperatures (©) at room 
temperature. 

(3) In the case of CrCl; the specific heat 
anomaly indicates a Curie point of 16.8°K. A 
field dependence of the susceptibility was found 
above this temperature, as well as below. The 
initial susceptibility increases with decreasing 
temperature and may be infinite, or very large, at 
13.9° if the demagnetization correction is applied. 

(4) FeCl, exhibits an S-type magnetization 
curve at H, temperatures. This magnetization 
curve is characterized by a very small initial 
susceptibility which decreases with decreasing 
temperature, followed by a very large differential 
susceptibility and eventual approach to magnetic 
saturation, both increasing with decreasing 
temperature. CoCl, and NiCl, exhibit this S-type 
magnetization curve in lesser degree. 

(5) MnCl, exhibits no deviations from the 
Curie-Weiss law, although its crystal structure is 
very similar to that of the other anhydrous 
dichlorides. 

(6) Many of the compounds exhibit a mini- 
mum in the reciprocal susceptibility vs. tempera- 
ture curve, but apparently this is not necessarily 
connected with any field dependence of the 
susceptibility. 

The authors express their appreciation for the 
encouragement and cooperation given by Pro- 
fessor F. G. Keyes, and their thanks for the aid 
rendered by Professor R. C. Young in the 
chemical preparation of the compounds, 
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A theoretical study is made of the experimental results reported in Part I. The predictions of 
* paramagnetic theory are compared with the room-temperature results, and an attempt is made 
to explain some of the departures from theory in terms of the crystalline field. The low tempera 
ture experimental anomalies require the development of a theory of metamagnetism. The as- 
sumption is made that at temperatures below the Curie point neighboring atoms behave as 
groups, which are spontaneously formed at the Curie point. The atomic magnetic moments 
within each group are assumed to be either parallel or antiparallel to each other, and the 
antiparallel arrangement is assumed to have the lowest energy. The theory developed on these 
assumptions satisfactorily explains the general behavior of the metamagnetic chlorides (CrCls, 
FeCl, CoCl2, and NiCl:), and also gives a qualitative explanation of the behavior of the other 
anhydrous chlorides which exhibit a minimum in the reciprocal susceptibility vs. temperature 
curve but do not show any field dependence of the susceptibility. 


INTRODUCTION 


HE results of the experimental investigation 
of the magnetic properties of all the 
anhydrous chlorides of the first transition series, 
the iron group, have been reported in the pre- 
ceding paper, Part I.! It is the purpose of this 
paper to present a qualitative theoretical analysis 
of these results, especially the anomalous mag- 
netization characteristics which are associated 
with the lowest temperatures used in the investi- 
gation. In order to facilitate the discussion of 
the results a brief review of present paramagnetic 
theory will be given. For a detailed treatment of 
this subject reference should be made to Van 
Vleck? and Stoner.’ 

The paramagnetic susceptibility of atoms in a 
given j state results from the statistical distribu- 
tion of the atoms among the various energy 
levels which are associated with the removal of 
the degeneracy of the j state by the applied 
magnetic field. The determination of the various 
energy levels and the removal of their degeneracy 


by applied fields or interaction effects is the basic 


* Contribution No. 458 from the Research Laboratory of 
Physical Chemistry. 

1 See preceding paper ‘‘Part I, Experimental.” 

2J. H. Van Vleck, The Theory of Electric and Magnetic 
Susceptibilities (Oxford University Press, New York, 1932). 

*E. C. Stoner, Magnetism and Matter (Methuen and 
Company, London, 1934). 





problem in determining the susceptibility of para- 
magnetic compounds. The normal or “ground” 
state of an atom is determined by the Hund 
stability rules. These rules specify that in the 
lowest term the spin momenta combine to give 
the maximum resultant value of s consistent with 
the Pauli exclusion principle, and then the 
orbital momenta combine to give the maximum 
resultant value of /; for an incomplete shell of 
electrons the lowest term 7=/—s when the shell 
is less than half full, and 7=/+s when the shell is 
more than half full. As a result of the various 
possible combinations of a given resultant s and I, 
j can assume the values /+s, /+s—1, ---|/—s|. 
This results in a natural multiplet of energy 
levels, in which the difference in energy between 
different j7 values is determined by the strength 
of the interaction between / and s. In the case of 
wide multiplets, the energy of transition between 
successive j states is very large compared to kT, 
the /s coupling is strong, and nearly all the atoms 
are in their ground state. In the case of narrow 
multiplets the energy of transition between 
successive j states is very small compared to kT, 
the Js coupling is very weak so that the orbital 
and spin momenta are quantized separately. It 
has been shown that the presence of an asym- 
metrical electrostatic field surrounding an atom 
may render the / moment partially or completely 
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ineffective by producing a preferred orientation 
electrostatically, and thus hinder the response to 
an applied magnetic field. This is similar to the 
Stark effect. The s moment is unaffected directly 
by the electrostatic field. In the case of para- 
magnetic solids the electrostatic ‘‘crystalline” 
field, due to the atoms surrounding the one under 
consideration, may thus be responsible for the 
“quenching” of the orbital moment, so that only 
the spin moment contributes to the magnetic 
susceptibility. 

The theoretical limiting cases for the suscepti- 
bility of paramagnetic atoms may be given by 
the following formulae. The susceptibility per 
gram mole Xm=(Nus*p?/3kT), where N is Avo- 
gadro’s number, ug is the magnetic moment of 
one Bohr magneton, k is the Boltzmann constant, 
T is the absolute temperature, and p is the 
effective number of magnetons per atom. 

Wide multiplets: hv(jj’)>kT 


' p=gliG+1))}, 


where g is the Landé splitting factor. 
Narrow multiplets: hv(jj’)<kT 


p=(1(1+1)+4s(s+1) }}. 


Spin only: p=[4s(s+1) ]}. 

In the case of wide multiplets, a temperature- 
independent susceptibility term should be added 
to the above susceptibility. If the multiplet 
width is very wide this term is negligible. When 
the multiplet width is comparable with kT, the 
statistical distribution of the atoms over all the j 
states must be considered and more detailed 
formulae must be used. These are presented in 
the detailed references. 

Table I contains the tabulation of the above 
limiting p values for the iron group ions. These 
may be compared with the experimental values 
found in Part I. The experimental p values are 
determined from the characteristics of the re- 
ciprocal susceptibility vs. temperature curve at 
room temperature. Excluding the uncertain p 
values for the titanium ions, several generaliza- 
tions may be drawn from Table I. The Mn? and 
Fe’ ions have a ®°S ground state and, theoretically, 
they should exhibit the spin-only value. The 
experimental p value is actually less than the 
spin-only value. This result cannot be attributed 
to experimental error. The excellent agreement 
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of the experimental p values of the Mn? and Fe? 
ions is a strong indication of the validity of the 
results. There is also no question concerning the 
chemical purity of the compounds (see Part I). 
Table I also shows that in the trichlorides, Cr* 
and V* have experimental » values that agree 
very well with the spin-only values, while in 
their corresponding dichlorides the experimental 
p values are considerably larger than the spin- 
only values. It would thus appear that the 
orbital magnetic moment is almost completely 
quenched in the trichlorides and only partially 
quenched in the dichlorides. When considering 
the effect of the number of electrons in the 3d 
orbit, it is of interest that although the theo- 
retical spin-only values are symmetrical about 
the half-filled shell (5 electrons), the experi- 
mental p values are larger for the ions in the 
second half of the shell (5-9 electrons) than in 
the corresponding ions (same s values) of the 
first half (1-5 electrons). 


CRYSTALLINE FIELD EFFECTS 


The theoretical treatment of the effects of the 
electrostatic crystalline field upon the orbital 
magnetic moment has been developed by Van 
Vleck.? The action of the crystalline field may be 
described as a removal of the degeneracy of the 
resultant orbital moment of the ion. As in the 
case of narrow multiplets separate quantization 
of the orbital and spin moments is assumed, 
since the crystalline field acts directly only on 
the orbital moment. If the decomposition of the 


TABLE I. Calculated and observed magneton numbers 
effective in the ions of the first transition series. N is the num- 
ber of electrons in the 3d orbit. 

















EFFECTIVE NUMBER OF BouR 
MAGNETONS 
Wipe Narrow 
GROUND Mutti- Mutti- Spin 
N State l s j PLETS PLETS ONLY Exp. | Ion 
1 Dae 2 % 3%} 1.55 3.00 1.73 os Ti® 
ss FF £ Fee 4.47 2.83 (2.05)| Ti 
2.85 | V* 
3 ‘Fin 3 3 3%} 0.77 5.20 3.87 4.15 | V? 
3.90 | Cr® 
4 Dy 2 2 O | 0.00 5.48 4.90 5.13 | Cr? 
5 %Ss2 0 5 5] 5.92 5.92 5.92 5.73 | Mn? 
5.73 | Fe? 
6 3D, 2 2 4 | 6.70 5.48 4.90 5.38/| Fe? 
7 ‘Fen 3 36 % | 6.64 5.20 3.87 5.29 | Co? 
8 *F, 3 1 445.59 447 2.83 3.32) Ni 
9 2Desre 2 1g 56 3.55 3.00 1.73 2.08 Cu? 
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orbital energy levels is very large, the ions will 
mainly reside in the lowest orbital level, and an 
applied magnetic field will not appreciably 
change the distribution among the levels. Since 
the susceptibility is determined by the change in 
the average energy per ion due to the applied 
magnetic field, the contribution of the orbital 
moment is thus effectively quenched. The spin 
moments remain effective although the freedom 
of the spins is hindered by the spin-orbit coupling. 
The magnitude of the deviation of the suscepti- 
bility from the spin-only value depends, to a 
first approximation, upon the ratio of the spin- 
orbit interaction energy (as determined from the 
free ion multiplet width) to the crystalline field 
splitting. 

The detailed calculations have been made by 
Schlapp and Penney‘ for the ions Cr’, Co*, and 
Ni?, which are all in F states (/=3). The crystal- 
line field was assumed to consist of a strong 
field of cubic symmetry and a superimposed 
weak field of rhombic symmetry. As will be 
shown later, this type of field would be expected 
in the anhydrous chlorides. The effect of the cubic 
field upon the F state (7-fold orbital degeneracy) 
is to split it into one singlet level and two triply 
degenerate levels, the two latter being resolved 
completely by the rhombic field. In the case of 
Ni? and Cr* the singlet level is lowest, and since 
the splitting due to the strong cubic field is 
large, the susceptibility approaches the spin-only 
value. In the case of Co? the triplet level is low- 
est, and since the splitting due to the weak 
rhombic field is small, the effects are more 
complicated and large deviations from the spin- 
only value would be expected. The Hamiltonian 
function in a magnetic field HT was assumed to 
be of the form 


D(x+y'+2') +AX?+By*— (A +B)2* 
+A(I-s)+usH(/+25), 


where A, B, D are constants specifying the 
crystalline field, the cubic term D predominating, 
and X is the constant of spin-orbit interaction. 
The theoretical analysis resulted in suscepti- 
bilities for Ni? and Cr* equal to the spin-only 
susceptibility multiplied by a factor depending 
on \/D. If p is the actual effective magneton 


*R. Schlapp and W. G. Penney, Phys. Rev. 42, 666 
(1932). 
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number, p, the spin-only value, and q is a positive 
constant, then this multiplying factor is 


(p/Ps) = (1—4A/5Dgq) for Ni?, and 
(p/ps) =(1—2A/5Dq)? for Cr’. 


The value of \, as determined from the free ion 
multiplet width, is equal to —335 cm™ for Ni?, 
and +87 cm for Cr*. From the experimental p 
values in Table I, Dg is equal to 710 cm~ for Ni?, 
and —4700 cm“ for Cr*. These values of Dg are 
widely different from those found by Schlapp 
and Penney‘ and Janes’ on more magnetically 
dilute salts. The value for Cr*® also has the 
wrong sign, since Gorter® has shown that if a 
metal ion is surrounded by an octahedron of 
negative charges (as is the case here), D should 
be positive. 

Theoretical studies of the crystalline field 
effect in other ions of the iron group have been 
carried out by Jordahl’ and Siegert.* The in- 
version of the splitting of the orbital levels due 
to a cubic field for pairs of ions with the same / 
value, as in Ni® and Co’, has been shown by 
Van Vleck® to be characteristic of any pair of 
ions for which the number of 3d electrons adds 
up to 5, 10, or 15, such as (Co?, Ni?), (Fe*, Cu’), 
(Co?, Cr*), (V*, Ni*); and the splitting is such 
that the quenching of the orbital moment should 
be largest for the second member of each pair 
provided the spin-orbit interaction is the same. 
However, for a given / and s, the multiplet 
width is proportional to the j value, and is, 
therefore, largest for the ions in the second half 
of the iron group. The larger spin-orbit interac- 
tion may account for the fact that the experi- 
mental p values for the ions in the second half 
of the group deviate very much more from the 
spin-only values than those of the first half. 

In the above treatment of the crystalline field 
effects, the possible influence upon the suscepti- 
bility of the dipole-dipole and exchange effects 
has been neglected. Since the crystalline field can- 
not account for the low temperature anomalies 
observed in the anhydrous chlorides (see Part I), 
and as the concentration of magnetic ions in these 


5R. B. Janes, Phys. Rev. 48, 78 (1935). 


6 C. J. Gorter, Phys. Rev. 42, 437 (1932). 
70. M. Jordahl, Phys. Rev. 45, 87 (1934). 
8 A. Siegert, Physica 3, 85 (1936); 4, 138 (1937). 
*J. H. Van Vleck, Phys. Rev. 41, 208 (1932). 
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compounds is very large, it is probable that the 
neglected interactions must be considered even at 
room temperatures. 

Van Vleck and Penney’ have discussed the 
fact that Mn? and Fe’ show less than spin-only 
values. Although they attribute this discrepancy 
to experimental error they suggest a number of 
possible causes of such an effect. As has been 
mentioned above, this is not a case of experi- 
mental error. Dipole-dipole and exchange inter- 
actions are considered too small to be a satis- 
factory explanation. The most probable suggested 
explanation assumes the presence of incipient jj 
coupling (arising from the coupling between the 
spin moment of each electron and its own orbital 
moment). Under such circumstances it is not 
sufficient to consider only the Russell-Saunders 
coupling (between the resultant s and / of the 
ion), and the Hund rule is not valid. Van Vleck" 
has pointed out that the energy of the ion in the 
crystalline field is a minimum when the resultant 
spin has its lowest possible value. The internal 
energy of the ion is a minimum when the re- 
sultant spin has its maximum value, i.e., when 
the Hund rule is valid. Therefore, the combined 
action of the crystalline field and the incipient jj 
coupling would result in decreased effective spin 
for the ion, and the experimental p value would 
be less than the spin-only value. The observed 
low temperature curvature of the reciprocal 
susceptibility vs. temperature curve of FeCl; 
(Fig. 1, Part I) may be a consequence of this 
crystalline field effect. As has been suggested 
above the crystalline field in the trichlorides is 
apparently very much more effective than in the 
dichlorides. This may account for the fact that 
MnCl, does not show the same curvature as 
FeCl;. The difference between the dichlorides and 
trichlorides may also be due in part to the 
difference in exchange effects between the mag- 
netic ions. In the trichlorides each metal ion has 
three metal-ion near neighbors, and in the di- 
chlorides there are six near neighbors. 


CRYSTAL STRUCTURES 


The crystal structure of about half of the iron 
group anhydrous chlorides have been determined 
10 J. H. Van Vieck and W. G. Penney, Phil. Mag. 17, 961 


(1934). 
uJ. H. Van Vleck, J. Chem. Phys. 3, 807 (1935). 
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from x-ray studies." These crystal structures are 
all of the layer lattice type; each layer of metal 
atoms is separated by two layers of chlorine 
atoms from the next metal atom layer. The 
dichlorides have the typical CdCl, structure, 
each metal atom layer forming a two-dimensional 
hexagonal network in which every metal atom 
has six near neighbors. The trichlorides are 
similar, except that a metal atom is missing from 
the center of each hexagon so that every metal 
atom has three near neighbors. The dimensions 
of these structures are given in Table II. CrCl, 
and FeC); differ slightly, in that the layers of Cl 
atoms are so stacked as to produce cubic close 
packing of the Cl atoms in CrCl;, and hexagonal 
close packing in FeCl. 

Each metal atom in these structures is sur- 
rounded by six Cl atoms at the corners of an 
octahedron. Three of these Cl atoms are in the 
layer above the metal atom layer, and three are 
below, so that the axes of the octahedron do not 
coincide with the hexagonal axes of the crystal 
as a whole. The octahedra are distorted slightly 
by the crystal forces, and the distortion is 
probably greater in the trichlorides than in the 
dichlorides. It can thus be seen that the crystal- 
line potential surrounding each metal atom is 
adequately described by a large cubic (octa- 
hedral) term and a small rhombic term. 


THEORY OF METAMAGNETISM 


The iron group anhydrous chlorides which 
have a positive Curie temperature (0) at room 
temperature exhibit a typically anomalous low 
temperature magnetic behavior (see summary, 
Part I). This behavior is characterized by (1) a 
large dependence of the susceptibility upon the 
field strength, (2) no apparent spontaneous 
magnetization, and (3) no large hysteresis effects. 
The first characteristic is not typical of para- 
magnetism, and the second and third are not 
typical of ferromagnetism. For this reason it has 
been suggested™ that the magnetic behavior 
of these compounds should be classified as 
“‘metamagnetism.’’ No adequate explanation of 


2R.W.G. Wyckoff, The Structure of Crystals (1931) and 
Strukturbericht, [II] (Akademische Verlagsgesellschaft, 
Leipzig, 1937). 

13 J. Becquerel and J. van den Handel, J. de phys. et rad. 
10, 10 (1939). 
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metamagnetism has been presented. Landau™ 
suggested that some of the properties of meta- 
magnetics are due to the presence of spon- 
taneously magnetized metal atom layers, oppo- 
sitely oriented. This is similar to the theories of 
antiferromagnetism suggested by Néel® and 
others.!® These suggestions, in a modified form, 
are developed in this paper into a theory which 
is applicable to the metamagnetic chlorides 
(CrCl3, FeCl2, CoCle, and NiCl.).!” 

The specific heat characteristics of the meta- 
magnetic chlorides (Part I, reference 7) indicate 
a Curie point anomaly of the cooperative type. 
Such an anomaly would result from a spon- 
taneous arrangement of magnetic moments either 
parallel, as in ferromagnetism, or antiparallel. 
This ordering process may occur either through- 
out the crystal as a whole (long range order) or 
in small groups, or clusters, of neighboring atoms 
(short range order). The fact that no spontaneous 
magnetization occurs below the Curie point 
temperature may be explained by assuming a 
group behavior, with either (1) a haphazard 
orientation of atom groups, each group having a 
permanent moment due to parallel alignment of 
the atomic magnetic moments within it, or 
(2) a zero resultant moment for each atom group 
due to some antiparallel alignment of the atomic 
magnetic moments. Because of the complexity of 
the problem it is not practically feasible to calcu- 
late the energies of all the possible magnetic 
moment configurations in a crystal, and thus to 
determine theoretically which of the above 
explanations is most likely. The answer to this 
question is suggested by the fact that the initial 
susceptibility of FeCl, decreases as the tempera- 
ture decreases. This indicates that, as the tem- 
perature is decreased, there is an increase in the 
number of atomic groups whose energy is un- 
affected by the application of a magnetic field. 
Since a decrease of temperature increases the 
population of the lowest energy levels, the 
assumption (2), that the atomic groups with zero 


14. Landau, Phys. Zeits. Sowjetunion 4, 675 (1933). 

% L. Néel, Ann. de physique 18, 5 (1932). 

16F, Bitter, Phys. Rev. 54, 79 (1938). 

17 After this paper was written, Professor Van Vleck 
brought to the author’s attention an unpublished paper by 
J. Becquerel, presented in May, 1939, at the Strasbourg 
conference on magnetism of the Institut International de 
Cooperation Intellectuelle, in which similar considerations 
are applied to the problem of paramagnetic rotatory power. 
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TABLE II. Layer lattice dimensions of the iron 











group anhydrous chlorides 
METAL Atom DISTANC 
(ANGSTROMS) “ aie or 
BETWEEN METAL NEIGHBOR 
BETWEEN NEAREST Ion DISTANCE To 
CoMPOUND LAYERS NEIGHBORS} RAapIUs Ion Rapivus 
CrCl; 5.76 3.46 0.70 4.95 
FeCl; 5.75 3.42 0.67 5.10 
MnCl, 5.86 3.682 0.83 4.43 
FeCl, 5.84 3.580 0.80 4.48 
CoCl, 5.79 3.546 0.78 4.55 
NiCl, 5.78 3.544 0.74 4.79 














resultant moment have the minimum energy, is 
substantiated. CrCl; does not behave the same 
as FeCl,, but an explanation of this effect, not 
in contradiction with the above, will be presented 
later. 

The theoretical development will be based 
upon the simple model suggested by the above 
analysis of the general character of the experi- 
mental results. Below the Curie point the atoms 
are assumed to behave as groups. Within each 
group the atomic magnetic moments may be 
parallel or antiparallel. The antiparallel con- 
figuration produces a zero resultant moment for 
the group and has the minimum energy. The 
resultant moment produced by the parallel con- 
figuration is assumed to be in the same or 
opposite direction as an applied magnetic field. 
Thus, in the presence of an applied field there 
are three group energy levels, the lowest one of 
zero resultant moment, a higher one with the 
effective resultant moment parallel to the field 
direction, and the highest one with the effective 
resultant moment antiparallel to the field direc- 
tion. Above the Curie point the spontaneously 
formed groups no longer exist, although the 
interaction between near neighbors is such that 
the antiparallel and parallel energy levels are 
applicable to the atoms individually. The further 
assumption is made that the average field acting 
on any one group due to the neighboring groups 
is proportional to the average magnetization of 
the crystal. Finally, only one zero moment 
configuration is assumed to exist. This model is 
obviously oversimplified by the restrictions upon 
the direction of the resultant moments and upon 
the number of zero moment configurations. 
However, from the point of view of the analysis 
and general consideration of the experimental 
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results, such an idealized model has many ad- 
vantages as compared with a more rigorous and 
elegant theoretical treatment. 

Consider the above model at a temperature 
below that of the Curie point. Designate N the 
number of atoms per group, u the atomic mag- 
netic moment, o the average magnetization per 
atom (averaged over the whole crystal), ¢ the 
energy difference between groups with a re- 
sultant moment (parallel atomic moments) and 
groups with zero resultant moment (antiparallel 
atomic moments), and H the magnetic field 
applied to the groups. The average magnetization 
per group, 


No=kTa(log Z)/0H, 
where the partition function 


ZH=1 pe NuMD kT 4 p(t NBH)IET, 
Let a=e/kT and B=NyuH/kT. Then 


o 2 sinh B 
~= — —_, (1) 
uw e*+2 cosh B 


This is plotted in Fig. 1 for various values of a. 
When a=-—~© this equation reduces to the 
usual hyperbolic tangent relationship. The mag- 
netic field acting on each group is assumed to be 
H=H,+De where H, is the applied field and D 
is the internal field coefficient. The differential 
atomic susceptibility is 


do Nuo(ctnh B—o/u) 
x=—= 5 
dH, kT—DNuyuo(ctnh B—o/u) 





(2) 


At a temperature above that of the Curie point 
Eq. (1) is still valid, but the atoms now act 
individually, N=1, and e is the difference in 
energy between the parallel and antiparallel 
alignment of an atom with its near neighbor. 
When £ is large (>2), Eq. (1) reduces to 


o/u=1/(er*+1). (3) 


This equation is plotted in Fig. 2, and repre- 
sents the general S type form of the magnetiza- 
tion curve at low temperatures. Since a—£ 
=(e—NuH)/kT, the effect of decreasing tem- 
perature is to increase the step-like character of 
the curve of o vs. H, so that at absolute zero the 


magnetization changes abruptly from zero to 
saturation when NuH=e. For this value of the 
field the magnetization is independent of tem- 
perature, and the curves of o vs. H for different 
temperatures should all cross at the point where 
¢ is one-half saturation. This is approximately 
what was found for FeCl, (Fig. 2, Part I). 

When the applied field approaches zero, 8-0, 
and the initial susceptibility is 


xi=9/[D(T—9)], (4) 


where 0=2Ny?D/k(2+e*). The usual form of 
expressing the atomic susceptibility is x(T—@) 
=(pus)*/3k where p is the effective number of 
Bohr magnetons and uz is the magnetic moment 
of the Bohr magneton. In comparison with the 
usual form, © is as given and 


bus=ul6N/(2+e*) }}. 


Because of the presence of e*, @ and p are func- 
tions of temperature, decreasing as the tempera- 
ture decreases. Eq. (4) may also beYput into 


the form mand 
1 Nu? 2+e* Nu?D 





(4a) 
Xi € 2a € 


Equation (4a), without the constant term 
NwD/e, is plotted in Fig. 3 where RT /e=1/a. 

For high temperatures or small a, Eq. (4) 
reduces to 


xi(T+6€/3k—2Ny2D/3k)=2Nyu?/3k. (5S) 


If NywD/e is greater than 3, then the © as 
measured at high temperatures will be ‘positive, 
as is found in the metamagnetic chlorides. The 
dotted line of Fig. 3 illustrates Eq. (5) when 
D=0. 

If the difference in energy per atom between 
a parallel and an antiparallel arrangement is 4, 
then e= N6. Figure 3 is then the plot of y?/xé vs. 
kT/Né. Assume that 4 is independent of the size 
of the group. Above the Curie point the atoms 
act individually, N=1, and k7/é is large. The 
applicable part of the reciprocal susceptibility 
curve has a positive slope (right-hand part of 
Fig. 3). At the Curie point, group formation 
spontaneously occurs, N is large, and k7/e is 
now a fraction of its former value. This places 
the reciprocal susceptibility in the region of 








Cc. 








Fic. 1. The 
relative mag- 
netization as 
a function of 
8B, for vari- 
ous values of 
Qa. 














negative slope (left-hand part of Fig. 3). Such a 
process would explain the reciprocal suscepti- 
bility characteristic of a metamagnetic chloride 
such as FeCl, (Fig. 1, Part I). The whole of 
Fig. 3 may be applicable to those anhydrous 
chlorides whose experimental characteristics are 
similar to those of chromium, nickel, copper and 
vanadium chlorides shown in Part I, Fig. 1. 
These will be discussed later. 


COMPARISON WITH EXPERIMENT 


The constants used in the theoretical develop- 
ment may be approximately determined from the 
experimental results by means of the appropriate 
theoretical equations. Below the Curie point yu, 
N, and ¢ are assumed independent of temperature 
and field strength. The experimental values for 
FeCl, (Table II, Part I) result in «/k=59°. 
With this value of «/k the magnetization curves 
determined from Eq. (1) are plotted in Fig. 4 
for the two temperatures, 20.4° and 13.9°. As was 
noted above, the two curves cross for 


_ NuH=e=59k and oa/p=}. 


Figure 4 should be compared with the experi- 
mental curves for FeCl: of Fig. 2, Part I. 

The outstanding differences between the theo- 
retical and experimental curves are the accentu- 
ated step-like character of the 13.9° experi- 
mental magnetization curve, as compared with 
the theoretical, and the fact that the 20.4° 
experimental curve has an apparent saturation 
value appreciably less than that of the 13.9° 
curve. 

If only the spin magnetic moment is effective 
in these compounds h=4yg for the Fe? ion. 
The saturation moment per gram mole is then 
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22,200 and the theoretical intersection point of 
the two curves is 11,100. The experimental] 
intersection point is 11,700, in good agreement 
with the theory. The experimental initial sys. 
ceptibilities give N=9 atoms per group, approxi- 
mately. The value of H at the intersection point 
is then 24,500 gauss. The applied field at the 
intersection point H,=13,500 gauss. The differ- 
ence may be due to the internal field or to the 
fact that N and «¢ are not independent of tem- 
perature and field strength as has been assumed. 
It is also possible that the magnetization of the 
groups with an effective resultant moment in- 
creases with decreasing temperature, as in the 
theory of ferromagnetism. The possible ferro- 
magnetic behavior of the small groups postulated 
above is not known. It may be of interest in this 
connection that Peierls'* has shown that in the 
Ising model, which is similar to the above, all 
lattices except the linear chain may be ferro- 
magnetic. 

Information as to the nature of the interaction 
effects which give rise to the internal field 
coefficient D above the Curie point may be found 
from its value. The high temperature © = 2y°D/3k 
approximately, as can be seen from Eq. (5). For 
FeCl,, @=48° at room temperature (Part I, 
Table I). The experimental result is D=7X10*. 
Actual calculation of D assuming magnetic 
dipole-dipole interaction in a two-dimensional 
hexagonal network gives a value of 10” approxi- 
mately. The Lorentz field (spherical cavity) 
value is even less. It is apparent, therefore, that 
the actual interaction is more than 10 times as 
strong as magnetic dipole-dipole interaction. In 
view of the small distances between near neigh- 
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18 R. Peierls, Proc. Camb. Phil. Soc. 32, 477 (1936). 
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bors in these lattices (Table II), it is probable 
that the exchange interaction is the principal 
cause of these effects, as in ferromagnetics. 
Since the ratio of nearest neighbor distance to 
the ion radius in these compounds is only slightly 
greater than in the ferromagnetic metals (ratio 
is about four in Ni), it appears reasonable that 
the exchange interaction should be appreciable. 

The magnitude of the specific heat jump at 
the Curie point (Part I, reference 7) is of interest. 
The experimental values in cal./mole were 4.54 
for FeClse, 1.50 for CoCls, 0.66 for CrCl;, and 
0.13 and 0.18 for NiCl». If only the spin moments 
are effective, the theoretically predicted values 
for a simple ferromagnetic model* are 4.57 for 
FeCle, 4.38 for CoCls, and CrCl;, and 3.97 for 
NiCle. The agreement in the case of FeCle is 
very good but for the others the experimental 
value is too low. This latter discrepancy may be 
due to the incipient formation of groups above 
the Curie point, so that the total change is 
spread out. The specific heat anomaly is approxi- 
mately the same for an antiparallel group forma- 
tion as it is for a parallel group, since the entropy 
change is of the same magnitude. 

The behavior of CoCl, and NiCl, is deities to 
that of FeCl., although the effects are appreci- 
ably smaller. Calculations, as above, for CoCl: 
result in an ¢/k of 28°. The predicted intersection 
point magnetization per mole is 8300 and the 
experimental is 5900. The initial susceptibilities 
give N=5. The resulting intersection value of 
H= 28,000 gauss, while the applied field at this 
point is /7,=14,800 gauss. Similar calculations 
for NiCl, result in an €/k equal to 24°. No inter- 
section point can be chosen from the results. 
The initial susceptibility gives N= 3. 
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The behavior of CrCl; (Figs. 1 and 2, Part 1) 
differs from that of the other metamagnetic 
chlorides. The magnetization characteristics of 
CrCl; require that @ be either very small or 
negative, as can be seen from Fig. 1. The fact 
that the slope of the experimental reciprocal 
susceptibility vs. temperature curve is constant 
down to low temperatures indicates that 0 is 
independent of temperature above the low tem- 
perature region. This requires that e* be constant 
or negligible if Eq. (4) is to be applicable. The 
behavior of CrCl; may, therefore, be made 
consistent with the theory by assuming that 
e/k=0 or — ©. Experimentally a field of 30,000 
gauss at 13.9° results in a magnetization approxi- 
mately 0.8 of saturation. This gives N=4 for 
e/k=0, and N=2 for e/k= — ~. Since the Curie 
point from the specific heat anomaly is 16.8° 
and a field dependence occurs at 20.4°, it is 
probable that incipient group formation occurs 
above the actual Curie point temperature. 

The anhydrous chlorides CrCl.,, CuCle, VCl» 
and VCl; (Fig. 1, Part I) are apparently not 
metamagnetic. However, they all do exhibit the 
behavior indicated by Fig. 3. Group formation 
apparently does not occur in these compounds 
and, therefore, they have no Curie points. The 
applied magnetic fields are only sufficiently large 
to indicate the beginning of paramagnetic satura- 
tion. These compounds all show the negative 0 
predicted by Eq. (5). The theory and experiment 
do not agree quantitatively, but that may be 
due to the simplicity of the theory and the fact 
that other effects, such as the crystalline field, 
while comparatively negligible in the meta- 
magnetic chlorides, are not negligible in these 
compounds. 

The theory does not explain the absence of 











large interaction effects in MnCl, and FeCl. 
The crystal structures of these compounds are 
the same as that of the corresponding metamag- 
netic chlorides. The only apparent point of 
difference is that the ground state of these ions 
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is a 8S state. This would indicate that the ex. 
change interaction is negligible for a °S state, 
even though the electron shell is incomplete. 
The author thanks Professor F. Bitter and 
Professor J. H. Van Vleck for their suggestions. 
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The conditions in an investigation of the silent electric discharge from sharp points and fine 
wires in hydrogen at atmospheric pressure are such that the discharge itself purifies the gas from 
electron-attaching impurities, so that the conditions of the discharge are productive of free 
electron conduction. Measurements are made of the generalized coefficient of mobility for the 
drift of free electrons in hydrogen, and of the mobility of hydrogen positive ions in hydrogen. 
Past observers have usually had doubts concerning the purity of the gas prepared by chemical 
means, but the results here reported involve a demonstrably sustained purification during the 


measurement itself. 


HEN the point is negative in a point-to- 
plane arrangement of electrodes, it is sup- 

posed that the silent discharge is caused by 
electrons arising at the tip of the point and 
producing ionizations by collision in avalanches, 
until the avalanche fronts have moved away 
from the tip of the point to a place where the 
electric field is no longer sufficient to support 
ionizations by collisions. The number of electrons 
freed by ionizations in an avalanche may be con- 
sidered to be equal to 2" where N is the number 
of ionizations an original electron makes. By con- 
sidering directions of propagation of an avalanche 
increasingly removed from the axis of symmetry 
such as, for example, along the direction of the 
line, B, in Fig. 1, the electric field is less intense at 
any given distance from the metal than along a 
direction such as A, and so the number of 
ionizations that the original electron can produce 
is less. For example, if the number of ionizations 
along a direction A, is 30, while the number along 
a direction like B, is 20, an avalanche along the 
direction B will produce less than 1/1000 as much 
charge as an avalanche along the direction A 
(i.e., in the ratio 2!°). Thus, ionizations which are 
responsible for most of the current occur within a 





region approximating a hemispherical shell, such 
as CC. The region between the point and the shell 
and including the shell may be called the 
“ionizing sheath.”” In hydrogen, this sheath ap- 
pears to be about 0.2 cm thick. Such a hemi- 
spherical cap is generally observed in this form of 
discharge. 


STABILITY OF SHEATH 


A rather detailed discussion of the negative 
point-to-plane discharge in air has been published 
by Trichel, and by Loeb and Kip.! 

These observers report on a pulsating character 
of the negative corona discharge, which they 
attribute to the attachment of electrons to oxy- 
gen molecules to form negative molecular ions, 
and, due to the much smaller mobility of ions 
than electrons, the suppression of further ioniza- 
tions in the sheath until ions have had time to 
move away. Their supposition is that each step 
in the sequence: avalanching, electron attaching, 
movement away of ions, et seg. must advance well 
toward completion before the next step in the 
sequence can begin. 


1G. W. Trichel, Phys. Rev. 54, 1078 (1938); L. B. Loeb 
and A. F. Kip, J. App. Phys. 10, 142 (1939). 
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MOBILITIES IN HYDROGEN 


The writer also made a study of the corona 
discharges in air several years ago, and had 
observed a pulsating character of the discharge 
both on the negative and also on the positive 
polarities. 

In our measurements of the negative discharge 
in hydrogen, observations were made to detect 
any pulsations which might occur. A circuit 
similar to Trichel’s was used, except that in most 
of our work both on hydrogen and air, as well as 
other gases, a higher amplification was available 
and the signal was taken from a high frequency 
inductance in series with the discharge circuit. In 
this way, the bias of the first grid of the amplifier 
circuit was not thrown off by fluctuations of cur- 
rent during unsteady parts of a characteristic. It 
was found that a hydrogen discharge initially 
has more rapid and smaller pulsations than the 
pulsations in an air discharge, and occasionally 
these come in regular sequences of pulses. These 
pulsations gradually disappear with operation, 
and by the time they have disappeared, the cur- 
rent for any given voltage has risen to much 
higher values. The steady pulsation-free dis- 
charge current in purified hydrogen must be a 
space charge limited current with the greater 
density of charge near the ionization sheath at its 
outer boundary. 

If the tube was allowed to stand idle for a few 
hours, the reapplication of potential showed that 
the tube had become partially deconditioned. It 
was necessary to pass current for a short time 
before the tube returned to the high current con- 
dition which had been attained by sustained 
operation previously. This reconditioning, how- 
ever, did not require nearly so long as it did 
initially. Allowing the tube to stand idle again, 
the tube again might become deconditioned and 
if so, it could be reconditioned by operation for a 

still shorter time. The tube eventually reached a 
condition such that it remained in the high cur- 
rent condition even when allowed to stand idle. 
This showed that the deconditioning of the tube 
was not due to any leak. It was discovered that 
the initial conditioning of the tube was involved 
in some way with a sharpening of the tungsten 
wire ends. Examination with a microscope showed 
that the end of the tungsten wire was etched into 
a regular conical shape, and if, for example, a 
wire 0.0075 cm diameter was being used, the tip 
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of the conical point had a radius of curvature 
which was less than 0.0010 cm and was usually 
of the order 0.0005 cm. 

It had been observed in our earlier work on the 
corona discharges in air, that for points having a 
radius of curvature of 0.0010 cm or less, and of 
suitable metal to remain with that small a 
diameter, the pulsating character of the discharge 
reported by Trichel, Loeb, and Kip! did not 
occur. There were in air, however, some much 
smaller pulsations of an irregular character which 
did occur with these extremely sharp points. The 
fact that the discharge in hydrogen loses its 
pulsating character with the production of the 
extremely sharp conical tip, is probably due to 
the same thing, namely, the very sharp tip of the 
point. However, this sharpening of the tip of the 
point does not completely explain the decon- 
ditioning with standing, because the tip of the 
point did not change its shape any simply by 
turning off the current. It has been well estab- 
lished in high vacuum thermionic investigations? 
that tungsten in the presence of both oxygen and 
hydrogen acquires a mono-molecular layer of 
oxygen over which there is a mobile partial 
second layer of oxygen. The continued operation 
of the tungsten at an elevated temperature does 
not remove this oxygen except by the slow 
formation of tungsten trioxide and all of the 
oxygen is eliminated from the gas phase by this 
kind of reaction before the oxygen mono-molec- 
ular layer begins to disappear. In other words, 
this is an extremely quantitative means for 
eliminating oxygen from the gas phase. 











2 Langmuir, J. Chem. Soc., 511-543 (1940). 
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WIRE DISTANCE POINTS VOLTAGE 
CURVE DIAMETER TO PLANE SUBTRACTED 

1 0.0010 cm 0.27 cm 520 volts 
2 1.18 590 
3 3.10 660 
4 10.0 810 
5 0.0025 0.29 570 
6 0.98 645 
7 3.00 720 
& 10.0 880 
9 0.0075 0.38 790 
10 1.04 880 
11 3.06 970 
12 10.0 1170 
13 0.020 0.32 1080 
14 1.67 1200 
15 3.24 1320 
16 10.0 1620 


Fic. 2. Log-log plots for point negative. 


In the present corona investigations in hydro- 
gen, the faintest trace of oxygen would go to the 
tungsten wire end, and adhere so as to change the 
work function of the tip of the point. This change 
of work function at the tip of the point where 
field conditions are so extremely critical for the 
initiation of the avalanches responsible for the 
conductivity is observed as the conditioning of 
the tube up to the high current condition (or as 
explained later, a decrease in potential drop in 
the ionizing sheath) and so, the high current 
operation of the tube is really a very sensitive 
and reliable index of whether the oxygen has been 
removed from the gas phase. It seems reasonable 
to suppose that the other electron-attaching im- 
purities should behave more or less like oxygen 
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: WIRE DISTANCE POINTS VOLTAGE 
CURVE DIAMETER TO PLANE SUBTRACTED 
17 0.0010 cm 0.27 cm 590 volts 
18 1.18 690 
19 3.10 810 
20 10.0 950 
21 0.0025 0.29 640 
22 0.98 740 
23 3.00 860 
24 10.0 1100 
25 0.0075 0.38 900 
26 1.04 1000 
27 3.06 1120 
28 0.020 0.32 1350 
29 1.67 1500 
30 3.24 1650 


Fic. 3. Log-log plots for point positive. 


with respect to adsorption on tungsten, and so 
the disappearance of pulsations accompanied by 
the rise of the current to the high current con- 
dition is a reliable index of whether the con- 
duction is really free electron conduction. The 
effects of sharpening of the points by etching are 
being separated from the effects of work function 
at the tip in some investigations now being made, 
and it is expected that a detailed report of these 
two different effects can be made at a later date. 

Brass targets were used in most of this as well 
as all of the later work with cylinders because it 
was found that a harder surface sometimes 
acquired specks of material which gave a back 
ionization when the wire was negative and this 
back ionization might produce serious errors in 
measurements of mobilities. 
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MOBILITIES 


THE CURRENT VOLTAGE RELATION 


Although the shape of electrodes used in the 
first half of the work was that of a fine wire 
discharging toward a plane electrode which does 
not permit an exact solution for the space charge 
limited current, a dimensional argument can be 
used as follows. It will be supposed that the drift 
velocity of electrons is proportional to the mth 
power of the electric field intensity. 


u=KE", 


where u is the drift velocity, E is the electrical 
field intensity and K is a generalized coefficient 
of mobility. There is supposed to be a potential 
drop between the wire and the hemispherical shell 
V, so that the potential difference applied be- 
tween the shell and the collecting plane equals 
the total applied voltage Vi, minus the sheath 
drop Vo. 

Between the shell and the plane, the charge 
satisfies Poisson’s equation. 


p=V°V/4r. 


It is supposed that for some value of current and 
voltage, there exists a solution, V(x, y, z). Multi- 
plying the potential at every point by a constant 
multiplier \ gives a function \- V which satisfies 
Poisson’s equation and which satisfies boundary 
conditions where the potential applied between 
the sheath and the plane is greater by the same 
constant multiplier, A(Vi—Vo). The charge 
density is likewise proportional to V and so is 
increased by the ratio \ to dA-p. 

The total current flowing can be found by 
integrating the current density over any equi- 
potential surface. If 7 is the current density so 
that j=p-u, the total current is: 


I= fi-as= ff fo-was= ff [o-K-r-as, 


which gives the value of the current for the con- 
ditions for which it was assumed that there 
exists a solution. For the circumstances resulting 
in the potential distribution \V the current is 


= f froxazyas=r f fip-K-Bras, 


The value of the current when the potential 
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applied across the gap was increased by the 
factor, A, is 
h = AetY 
or 
IT=a( Vi- Vo)"*! 


gives the current as a function of the applied 
voltage, where a is a proportionality constant. 


DISCHARGE FROM SHARP POINTS 


The data plotted in Figs. 2 and 3 were taken 
with tungsten wire in hydrogen after continued 
operation at 0.1 ma rising to about 1 ma for each 
set-up, or until pulsations had completely disap- 
peared and the current had reached a steady 
value. Wires with diameters of 0.001, 0.0025, 
0.0075, and 0.02 cm were used at various dis- 
tances from point to plate, 0.3, 1, 3 and 10 cm. 
The diameter of the plane target was 15 cm so 
that the deviation from a point-to-plane arrange- 
ment, as well as the erratic effect of charges on 
the walls, distorts somewhat the observations for 
the 10 cm distance. It was further observed that 
the sheath did not present the smooth regular 
appearance with the 0.02 cm wire that it did with 
the finer wires. There appeared to be a tendency 
toward channeling through the sheath at the 
larger current for the largest wire. The assump- 
tion. of a substantially constant sheath drop is 
probably in error for this largest wire. 

It will be seen from Fig. 2 that the data 
generally can be adjusted in a log-log plot to 











Fic. 4. Diagram of tubes. 











straight line relations only by subtracting 
amounts assumed for the sheath drop which gives 
slopes of 3. The current voltage relations have the 
form: 


I_=a(Vi- Vo)?. 


Observations were also made when the point 
was positive with respect to the plane, and these 
data are plotted in Fig. 3. Here the only good fits 
to straight lines could be obtained by assuming 
sheath drops yielding slopes of 2. The current 
voltage relations are 


I, =B(Vi— Vo)’. 


These observations show that the conduction 
in hydrogen which is purified by sustained dis- 
charge from negative tungsten points until the 
pulsations have disappeared, is due to electrical 
charge whose drift velocity is given by an 
expression of the form. 


u=KE'* 


and that the discharge when the point is positive 
is due to charge whose drift velocity is given by an 
expression of the form: 


v= CE. 


It should be noticed that the value of sheath drop 
changes very little, if at all, for relatively very 
considerable changes in wire size, and target 
distance except for the largest size and largest 
target distance where irregularities of obvious 
nature can cause very considerable distortions of 
apparent sheath drop. 


DISCHARGE FROM A STRAIGHT WIRE 


Measurements were made of the negative and 
positive discharges from 0.0075 cm _ straight 


TABLE I. Values of K calculated from the data. 











I Ve Vi Ve-—Vi K 
(ma) (VOLTS) (VOLTS) (VOLTs) (E.s.U.) 

10 3940 2395 1545 6.8 x 10° 

8 3660 2330 1330 6.8 

7 3540 2290 1250 6.6 

6 3360 2265 1095 6.9 

5 3160 2200 960 7.0 

4 2950 2160 790 7.5 

3 2735 2100 635 6.7 

2 2600 2010 590 5.8 

1 2425 1895 530 3.4 

0.6 2300 1830 470 4.1 
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tungsten wires held axially in cylinders as shown 
in Fig. 4. Guard cylinders were used in order to 
obtain cylindrical symmetry. Cylinders with two 
different radii were used in order to separate out 
the conditions in the ionization sheath from those 
in the rest of the gap. The two arms of the tube 
were connected to insure the same composition 
of hydrogen in the tube for corresponding 
measurements. 

If the ionization sheath, with wire negative, 
has a radius 7 the field at the outer surface of the 
sheath is Eo, the potential drop through the 
sheath is Vo, and the electrons have a drift 
velocity in the rest of the gap proportional to the 
square root of the field: 


u=KE}, 


the current J to a guarded cylinder with radius r 
and length L, is related to the voltage between 
the wires and the cylinder by: 


21 \3 1 KLE,} 
V- v.-(—) (rr {14+-(1- ) 
KL 3 27 


KLE,!\? ro(ro+r) 
x—"—-_(1-——*) ++} 
r+(rro)' 18 ZI 





To 1 





r: 

For sufficiently large values of current, that is 
I>KLE,}, 

the series is convergent and can be written as 


V—Vo=(2I/KL)!r 
X {1—0.722ro/r—0.333(ro/r)i +--+}. 


The radius of the larger cylinder used was 2.70 
cm, and the smaller cylinder was 1.15 cm in 
radius. Each guarded section was 15.15 cm long. 
If the voltage applied to the larger cylinder is V2 
and that to the smaller cylinder Vi, when the 
same current is passing to the central cylinder of 
each, J, 


TABLE II. Calculated values of C. 











I V2 Vw Ve-Vi 

(ma) (VOLTS) (VOLTS) (VOLTS) (E.s.v.) 
0.5 7750 4115 3635 3.2 108 
0.4 7100 3900 3200 3.0 

0.3 6500 3650 2850 3.2 

0.2 5600 3340 2260 3.4 

0.1 4725 2990 1735 2.8 

0.01 3140 2500 640 2.1 
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Ve- Vo= (2I1/KL) ir, 
x { 1 —0.722r0/r2—0.333(ro/r2)§+ eee i. 


Vi— Vo=(2I/KL)in 
x {1 —0.722ro/r1 —0.333(ro/r1)'+ now a 


so the difference in these voltages is 


2I\3 0.33370! 
v.-ti=(—) (er) 1+ 
KL rery'—riret 
The mobility coefficient can be obtained from the 
data by 





I eT 0.5r3 


K= ; 
(V2—V,)! L retiit+rire! 


which is relatively insensitive to the assumed 
value of 7o. Putting an outside limit on ro of 


0.6 cm, 
0.579! 
1-+—————_ ¢ 1.05. 
rotyit+rire! 


It was observed that the discharge with the 
0.0075-cm wire negative was localized in spots, 
and that increasing the current produced more 
discharging spots and crowded them more 
closely together. 

Table I shows the data and calculated value 
of K. 

The values for K do not show a monotonous 
change with decreasing current so it appears 
that the assumption of J large enough was 
justified. However, the fluctuations in values of K 
for the smaller currents is easily understood from 
the randomness of likelihood that the discharging 
spots should fall inside the collecting cylinders if 
near an end. From the data for the five largest 
current values, 


K=6.8 X10 e.s.u. 


The current with wire positive, limited by 
space charge where drift velocity is proportional 
to field (v= CE) is related to voltage by 
V—Vo=(2I/CL)*{r—1.542ro 

+0.500r0(r0/r) +0.042ro(ro/r)?§— +++}, 


if I>CLE,*. Again in terms of voltages for the 
two cylinders for the same current as before, 


2I\3 0.579? 
v.-Vi=(—) (nr) 1— te], 
CL rife 
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Fic. 5. Comparison of the present value of 6.8 X 10° e.s.u. 
for the mobility of electrons in hydrogen with the results of 
Bradbury and Nielsen. The experimental points are from 
Bradbury and Nielsen. The dotted curve is calculated from 
the above value of the mobility. 


Table II shows the data and calculated values 
of C. 


DISCUSSION 


The results for the mobility of electrons in 
hydrogen may be compared with those of 
Bradbury and Neilsen* as shown in Fig. 5 where 
the dotted line represents the relation 


u=6.8X105E! (for p=760 mm, E in e.s.u.) 
or 
u=10.8X105(X/p)! (in volts/cm and mm) 


and agrees well with these data for field in- 
tensities up to X/p=5. 

Larger field strengths than about 4000 volts 
per centimeter produced arc-over so that a 
mobility of electrons at greater than about 
X/p=5 could not be measured. 

The result for the mobility of positive ions in 
hydrogen 


C=3.2X10' e.s.u. 
= 10.7 cm/sec. per volt/cm 


agrees much better with Bradbury’s‘ value of 
8.2 rather than with those of Mitchell, which 
Tyndall’ notes. 


( 936) E. Bradbury and R. A. Nielsen, Phys. Rev. 49, 391 
1 le 

‘N. E. Bradbury, Phys. Rev. 40, 519 (1932). 

& A. M. Tyndall, The Mobility of Positive Ions in Gases 
(Cambridge University Press, 1938), p. 35. 








DECEMBER 1, 1940 


PHYSICAL REVIEW 


VOLUME 5g 


Interpretation of Torsional Frequencies of Crystal Specimens 
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(Received September 12, 1940) 


A crystalline cylinder when subjected to a harmonically 
varying twisting moment in general not only twists but 
bends. A correct interpretation of the response frequencies 
is important in dynamical determinations of the elastic 
constants of crystals, and it is desirable that the measure- 
ments be made under such conditions that a simple 
formula, similar to that for isotropic materials, may be 
used without fear of serious error. In static measurements 
two such formulas are required, one for twist with bending 
prevented and the other for twist with bending unre- 
strained. It is shown here that each of these is a good 
approximation in dynamic measurements over certain 


ranges of frequencies and of specimen lengths, and formulas 
are derived by which the errors may be estimated. A 
simple straight-line plot based on approximate values of 
the elastic constants gives all the information that is 
needed for the selection of good working conditions in the 
final measurements. For the sodium and copper-gold 
crystals on which measurements have been reported, the 
pure torsion formula is a good approximation when no 
parasitic response frequencies are observable; but it would 
cease to be a good approximation for specimen lengths 
below about 3 cm. 





N a crystalline cylinder whose axis is oriented 
arbitrarily with respect to the crystal axes, 
application of either a bending or a twisting 
moment alone produces both bending and twist- 
ing. In static experiments the effective Young’s 
modulus E’ and rigidity G’ measured under these 
conditions differ from the values E and G 
obtained when only pure bending or pure twist is 
allowed to occur.! If the specimen axis is Oz and 
if the plane of the torsion-coupled bending is 
taken as the zx plane, then in terms of the elastic 
moduli s;;’ referred to the xyz axes 


F’'=E(1 —e) = 1/533’, 
G’ =G(1— 6) =2/(s44’ +555’), (1) 
€ = $34'?/S33' (Saa’ +555’). 


In dynamic measurements of the elastic con- 
stants of crystals, it has usually been possible to 
choose an orientation for which e=0; the inter- 
pretation of the observed resonance frequencies 
under torsional excitation then offers no diffi- 
culty.? In recent measurements on metal crystals 
this has not been possible.*4 The experimental 
results have been consistent with the assumption 
that the effective rigidity was G, but a theoretical 
analysis of the problem seems desirable. The 


1W. Voigt, Lehrbuch der Kristallphysik (Teubner, 1928), 
pp. 638, 734, 738. 

2 F. Rose, Phys. Rev. 49, 50 (1936); M. Durand, Phys. 
Rev. 50, 449 (1936). 

3S. L. Quimby and S. Siegel, Phys. Rev. 54, 293 (1938). 

4S. Siegel, Phys. Rev. 57, 537 (1940). 


formulas derived by Goens® are unsatisfactory 
because in their rigorous form they are incon- 
venient for numerical calculation, and the first- 
order correction derived from them is not always 
reliable. 


THEORY 


The elastic properties of the specimen are 
completely described by three quantities E or E’, 
G or G’, and e. Let E* and G* represent whichever 
values of Young’s modulus and the rigidity are 
chosen for this purpose. Then the actual prop- 
erties may be regarded as produced by a gradual 
increase of ¢ from zero to its actual value in an 
originally isotropic specimen with the same E* 
and G*. The actual frequencies emerge during 
this process from original unperturbed frequencies 
which may be classified unambiguously as 
flexural or torsional but whose values will depend 
on the meanings of E* and G*. Since the “free” 
frequencies (E* = E’, Gt =G’) may be as much as 
20 percent lower than the “pure’’ frequencies 
(E*=E, G*=G), it is clear that a first-order 
correction for the perturbation, if it is adequate 
at all, can be so only when the “unperturbed” 
frequencies are properly chosen. 

In terms of E and G, the differential equations 
of motion of a cylinder of density p and radius a 


5 E. Goens, Ann. d. Physik 15, 455 and 902 (1932). See 
ta" Kimura, Proc. Phys.-Math. Soc. Japan 21, 686 
1 . 
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may be written® 
po = — 1a*Ed'd/ds'+2Ca*y/dz?, 
pv = —}a°Cd'/dz'+Gary/dz?, 
C=(eEG/2)!. (3) 


where 


(Here rotatory inertia about the transverse axis 
Oy has been neglected.) The variables ¢ and y 
are the bending and twisting moments exerted 
across a cross section, divided by the moments 
of inertia }wa* and }2a‘ of the cross section 
about the corresponding axes. These are more 
convenient variables than the transverse dis- 
placement and angular twist because they make 
it possible to write the boundary conditions in a 
form not containing ¢. For a rod with free ends 
at z=+//2 the boundary conditions are 


o=~=0, d¢/dz=0; (4) 


only this case will be considered in detail. 

The differential Eqs. (2) may be written in 
three other forms, corresponding to other choices 
of E* and G*; these may be obtained from (2) by 
means of (1). If E* and G* are held constant but e 
is set equal to zero, the equations reduce to the 
form for isotropic specimens, and for vibrations 
of frequency p/27 


pp?o={a°E*d'g/dz*, pp>y= —G*d*y/dz*. (5) 


The corresponding flexural frequencies p;*/2r 
and torsional frequencies w;*/2m are given by’ 


raf E*\} n (G*} 
pt=ni—-(—) ’ w*=7(—) , (6) 
I22\ p IX p 


jis an integer, and n, is the ith nonvanishing root 


of 
tanh? (mxz/2) = tan? (mx/2) (7) 


(m,=1.506, nm;=i+} to four figures for 7>1). 
If the frequencies are plotted logarithmically as 
functions of specimen length /, the torsional fre- 
quencies lie on straight lines of slope —1 and 
the flexural on lines of slope —2; these are the 
fine lines (both solid and dashed) in Fig. 1. The 
free frequencies (E*= E’, p;* =p,’ etc.) are lower 
than the pure (E*=E, p;* =); etc.) in the ratio 
(1—e)}. 
® Reference 1, p. 673, or reference 5. 


7H. Lamb, Dynamical Theory of Sound (Arnold, 1910), 
p. 124. 
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The normalized eigenfunctions ¢;, ~; corre- 
sponding to p;* and w;* are 


¢,=l-} [cosh kz, ‘cosh (Ril, 2) 
—cos k,z/cos (kil/2)], todd, 

















, — 8 
=]-' [sinh k,z/sinh (k,l/2) (8) 
—sin k,z/sin (kil/2)], ieven; 
¥;= (2/1)! cos jxz/l, 7 odd, (9) 
=(2/l)' sin jxz/l, j even; 
k,=n<w/l. (10) 
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LENGTH OF SPECIMEN IN CENTIMETERS 


Fic. 1. Unperturbed and perturbed torsional and flexural 
frequencies for a hypothetical sodium crystal cylinder of 
the diameter used by Quimby and Siegel, with specimen 
axis in an orientation for which the coupling between 
torsion and bending is large. E=5.47 X 10", G=2.804 x 10" 
dyne/cm?*, «=0.3628. The fine solid and dashed lines repre- 
sent the unperturbed frequencies as functions of specimen 
length: the ’s are flexural, the w's torsional; the primed 
are “‘free,”” the unprimed ‘‘pure.” Of these, the ones repre- 
sented by solid lines are a better approximation to the 
actual frequencies, which are given by the heavy curves and 
the circles. The circles have been calculated by a rigorous 
but laborious method, the curves by the approximate 
methods of this article. The subscripts are the number of 
moment loops in the specimen; the strongest interaction 
is between torsional and flexural modes of vibration for 
which this number is the same. 


The solution for «#0 may be expressed in the 
form 


= Adi, Y=L By,j. (11) 


Equations determining the A,’s and B,’s may 
be found by substituting (11) in (2) (with 
0?/dt?= —p*), using the differential equations 
satisfied by the ¢,’s and y,’s, multiplying the first 
of the resulting equations by ¢, and the second 
by y¥., and integrating over the length of the 
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specimen. This gives 


(p?— px*)Az—(2C/G) DY cxjw;?B;=0, 
i (12) 
—(C/E) X pircivA + (p?—w,*)Bi=0, 


where 


1/2 
— f odds. (13) 
1/2 


For i and j both even or both odd, 
| Ciz| = (4V2/m)n;27/(ni*#— 7"); (14) 


the other c;;’s vanish. Thus Eqs. (12) break up 
into two sets. One set contains only odd 7’s and 
j’s and refers to modes of vibration in which ¢ 
and y are symmetrical about z=0 (see Fig. 2); 
the other contains only even 7’s and j’s and refers 
to antisymmetrical modes. There is no inter- 
action between a flexural mode of one type and a 
torsional mode of the other type. 

From (14), ¢:?=0.972, c2x?=0.952, css*=0.919, 

*+, ¢;70.810 as j--~; and since Siici? 

= > ci =1, the other c;;’s are small. Therefore a 
torsional mode interacts chiefly with the flexural 
mode for which the number of moment loops in 
the specimen is the same (1=/, Fig. 2), and only 
weakly with the others. As a first approximation, 
account may be taken of this strongest interac- 
tion alone by setting |c;|=1 for i=j, =0 
otherwise. Then Eqs. (12) separate into pairs ; the 
kth pair contains only A; and B,, and for non- 
vanishing values of these 
(p?— pr*) (P?— wi?) —epiror?=0, — (15) 
whence p*?= p;? or w,;” for p2>w,’?, and p?= p,” 
or w,” for p2<w,?. Thus if the strongly inter- 
acting frequencies are not too close together, the 
effect of their interaction is that the lower of 
them is approximately free and the upper ap- 
proximately pure. These frequencies correspond 
to the fine solid lines in Fig. 1. 

To obtain a better approximation, the best 
choice of unperturbed frequencies is p; and w;’ for 
lengths below that at which the frequency under 
consideration and its strongest perturbing fre- 
quency become equal ; p,’ and w; for lengths above 
this. The fact that the lengths actually used in 
the experiments lie in the latter range explains 
the success of the pure torsion formula. The 
following corrections to the unperturbed fre- 
quencies may be derived by standard methods 
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Fic. 2. Unperturbed flexural and torsional modes o 

vibration and relative degrees of interaction between them 

@ is proportional to the bending moment and y to the 


twisting moment. The double arrows indicate strong 
interaction, the single arrows weak interaction. 


from the equations that now replace (12); the 
range in which each formula is to be used will be 
evident from the frequencies occurring in it. 
Very near a perturbing frequency, provided 
the effect of other perturbing frequencies is 
negligible, 
p?=43[b+ (b’—4c)#], 


b= pi? + uP +nei7pi" 
or (16) 
b= pitw;?+nci7u;”, 


where 


and 
C= pi?w?=piw;”, 


n=e/(1—e). 
Far from all perturbing frequencies, 
(p?—w;?)/w? = 2X cij*pi'*/(wj?— pi’), 
(p?— pi’*)/pi? =n LD cij*pi'?/(i? —w/?), 
(p?—w;'*)/w;*=4 y Cij7w;'*/(wj'*— pi*), 
(p?— p;*)/p2=n r €i;*w;'?/(pi?—w;'*). 


2 


(17) 


Occasionally it is necessary to take account 
simultaneously of two perturbations, one due toa 
weakly interacting frequency too close for use of 
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(17), the other to a strongly interacting frequency 
still not distant enough to neglect. This is true for 
instance in Fig. 1 of frequencies produced by 
interaction of w; and p;’ near their crossing point, 
where ~:’ and ws still produce a considerable 
effect. Formulas for such cases need not be given 
here; in the application to experiments the 
unfavorableness of such a situation is evident 
without detailed calculation. 


APPLICATION 


In the experiments mentioned, the specimen is 
joined end-to-end to a quartz crystal equipped 
with electrodes, and torsional vibrations are 
excited piezoelectrically. The frequencies ob- 
served are the resonance frequencies of this 
composite system, but the lengths are so chosen 
that the junction of the two cylinders is approxi- 
mately a moment node for torsion. The fre- 
quencies are therefore almost the same as the 
natural frequencies of vibration of the specimen, 
as long as there is no torsion-flexion coupling. 
The effect of such coupling on the observed 
frequencies cannot be computed accurately by 
using the formulas derived above for the effect on 
the natural frequencies of the specimen by itself. 
If, however, that calculation shows that the 
natural frequency of the specimen alone differs 
negligibly from the “‘pure”’ frequency, it is safe 
to assume that the observed frequency also 
differs negligibly from the one that would be 
observed if bending could be prevented. 

The simplest procedure is therefore the follow- 
ing. From the observed frequencies, the elastic 
constants may be computed by means of the pure 
torsion formula; the values obtained will be at 
least approximately correct. From the elastic 
constants, the unperturbed frequencies may be 
calculated and plotted as in Fig. 1, for each 
specimen used. An inspection of the plot will 
reveal whether any of the experimental points 
fall in a range where the error is likely to be 
excessive, and in such ranges the error for a 
single free rod can be calculated by the formulas 
given above. Points for which the error is 
objectionably large may be rejected and the 
elastic constants recomputed without them; or 
the whole preliminary measurement and calcu- 
lation may be used merely as a means of selecting 
specimen lengths for the final set of measurements. 
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Figure 1 has been calculated for a hypothetical 
sodium specimen in a particularly unfavorable 
orientation, from the values of the elastic con- 
stants given by Quimby and Siegel.** The heavy 
curves are corrected frequencies of free vibration. 
These have been calculated by the method of this 
article. The circles are check points computed by 
means of a rigorous formula similar to that of 
Goens ; this calculation consists in solving a cubic 
equation for some value of p, calculating and 
plotting a function f of the three roots and of p 
and / for a number of values of /, and thus finding 
a combination of / and p for which f=0. The 
calculation for given 1 would be still more 
laborious. 

It can be seen that the fundamental torsional 
frequency is practically identical with the 
unperturbed “‘pure’’ frequency for lengths greater 
than 3 cm, except when the perturbing frequency 
is very close ; in this case two response frequencies 
would probably be observed and the data would 
automatically be rejected.* The only possibility of 
serious error, therefore, lies in the use of too short 
a length. For the lengths used by Quimby and 
Siegel one would expect the pure torsion formula 
to hold quite well, as they found it did. The 
length range for which the free torsion formula 
would be valid is of no interest, at least for 
specimens of this diameter, for at such short 
lengths the specimen can no longer be regarded as 
thin (a?</P). 

Siegel in his torsional measurements on copper- 
gold crystals used two specimens.*?° One had an 
orientation so close to [100] that the frequency 
correction at room temperature is only 0.004 
percent. In the other w2, which was measured, 
was rather close to p,’ and was still appreciably 
perturbed by 2’; the effects are 0.7 percent and 
0.3 percent, respectively, and additive. The 
analogous situation in Fig. 1 is a measurement of 
w, at about /=6. Here a correct choice of E* and 
G* is essential ; for if ps3 instead of p;’ is chosen as 
unperturbed frequency, the perturbation is 
shifted to the right and the whole calculation is 
erroneous. 


*For cubic crystals, in the notation of reference 3, 
1/E’=sy—2sT, 1/G’ =su+4sT, ¢=2s8E’G’ (I —4I*+-3,x). 

* Reference 3, footnote 12. 

1®The data required for this calculation were kindly 
provided by Dr. Siegel. 
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not hold itself responsible for the cpinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


A Study of the Frequencies v. and v, in the Infra-Red 
Spectra of SiH, and GeH, 


The infra-red spectra of the molecules SiH, and GeH, 
have been investigated by Steward and Nielsen! who 
report intense regions of absorption near 4.54 and 11.0u 
in both cases. These bands were found, under higher 
dispersion, to resemble those occurring in the spectrum of 
methane in a general way and were accordingly identified 
as the optically active frequencies vs and »4, respectively, 
in each case. These measurements revealed, moreover, 
that the rotational absorption pattern was much more 
complex than could be accounted for on the basis of the 
prevalent theory for such molecules. In the case of SiH, a 
secondary peak was found in the spectrum near 10.34 to 
which no real interpretation was given. 

In a series of papers by Childs and Jahn’ and by Jahn? 
this complex structure of v4 was accounted for on the basis 
that there exists a resonance Coriolis interaction between 
the frequency », and the optically inactive frequency v2 
both of which have nearly the same values here. These 
results have been verified more recently by Shaffer, 
Nielsen and Thomas‘ who have presented a theory of the 
rotation-vibration problem of the tetrahedrally symmetric 
xy4 type molecule accurate and complete to second order 
of approximation. The Coriolis interaction manifests itself 
in somewhat the following manner. If », is larger than v2 
the effective moments of inertia in these states will be, 
respectively, larger and smaller than the moments of 
inertia in the normal vibration state. This should give rise 
to convergence of the rotation lines in », toward higher 
frequencies and in the frequency v2, if it were active, toward 
lesser frequencies. If, on the other hand, v4 is smaller than 
v2 a situation just the reverse of the above should prevail. 

It was even suggested by Childs and Jahn with regard 
to the SiH, molecule where », and v2 are separated only 
by some 60 cm that a mixing of the wave functions of 
the two states might occur because of this Coriolis per- 
turbation, sufficient to induce an optical activity into the 


TABLE I. Values of the frequencies v2 and v4 in the spectra of SiH« 
and GeHa. 








SiH« GeHs 
v2 975 cm™ 932 cm™ 
va 910 cm™ 818 cm™ 
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otherwise inactive frequency ». In this manner, they 
suggest an explanation of the weaker maximum found by 
Steward and Nielsen near 10.3.4. 

In the light of what has gone before it has been con- 
sidered important to re-examine the spectra of the above 
molecules under the higher resolving power which is now 
practicable. 

In the case of the silane spectrum the absorption pattern 
verifies in a general way the pattern reported by S. and N. 
except that here almost complete resolution has been 
accomplished, separations of the order of 0.25 cm™ being 
frequently observed. In the present data the secondary 
peak near 10.34 is clearly discernible as a Q branch of a 
much weaker band extending to 9.0u. It is apparent that 
although the Q branch near 11.04 is unsymmetrically 
broadened toward smaller frequencies this Q branch is 
broadened unsymmetrically toward higher frequencies. 
This makes it possible with certainty to identify the bands 
at 11.04 and 10.34 with the frequencies », and w, re- 
spectively. 

With regard to the spectrum of GeH,, the present meas- 
urements again verify those reported by S. and N. near 
10u. They show, however, that there exists immediately 
beyond the region studied by S. and N. at longer wave- 
lengths a region of far more intense absorption. The con- 
clusion drawn is evidently that the region originally studied 
was not v4, but in reality vz. The complete region from 94 
to 13u has been measured and as in the case of SiH, the 
bands are almost completely resolved. They reveal, as 
before, a convergence toward higher frequencies in the 
band », and toward smaller frequencies in the band ». 

The results of this preliminary investigation are sum- 
marized in Table I. 

Complete details of these measurements with wave- 
lengths will be published in a later communication. 

J. W. STRALEY 
C. H. TrnpaL 


H. H. NIeELsen 


Mendenhall Laboratory of Physics, 
Ohio State University, 
Columbus, Ohio, 
November 7, 1940. 


+5 ‘enone and H. H. Nielsen, Phys. Rev. 47, 828 (1935); 48, 
861 (1935). 
2 W. H. J. Childs and H. A. Jahn, Nature 141, 916 (1938); Proc. Roy. 
Soc. A169, 451 (1939). 

3H. A. Jahn, Proc. Roy. Soc. A168, 469 (1938); A168, 495 (1938). 

4 W. H. Shaffer, H. H. Nielsen and L. H. Thomas, Phys. Rev. 56, 895 
(1939) ; 56, 1051 (1939). 
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Rotatory Power of Nickel Sulphate at Low Temperatures* 


An investigation of the natural and magnetic rotatory 
powers of crystalline a—NiSO,-6H:20 at liquid-air tem- 
perature has been undertaken in this laboratory. Pre- 
liminary results for several visible wave-lengths are 
reported below. 

The crystals were cooled by direct immersion in liquid 
“air” which was assumed to be largely oxygen. Observa- 
tions were made through the double cylindrical Pyrex glass 
walls of the container (outside diameter 1.3 cm). The 
rotations actually measured ranged from 1° to 7° and 
satisfactory agreement was obtained between visual ob- 
servations with a Lippich half-shade polarimeter, and 
photographic measurements by a method previously 
described.1 Control measurements indicated that no 
serious errors arose from double refraction in the glass 
walls, or from lack of parallelism of the light rays in the 
crystal. 

The results so far obtained, with corresponding room- 
temperature data for comparison, are given in Table I. 


TABLE I. Rotatory power of nickel sul phate. 








WAVE-LENGTH IN ANGSTROMS 4358 5461 5780 


Natural rotatory power, in degrees/mm: 
At 90°K 





+1.2 —2.0 —2.3 
At room temperature! +1.7 —0.75 —1.25 
Difference —0.5 —1.25 —1.05 
Magnetic rotatory power, in min./cm-oerst.: 
°K +0.072 | +0.044 | +0.034 
At room temperature? +0.044 | +0.026 | +0.022 
Ratio 1.64 1.69 1.54 














1F. G. Slack and P. Rudnick, Phil. Mag. (7) 28, 241 (1939); N. 
Underwood, F. G. Slack, and E. B. Nelson, Phys. Rev. 54, 355 (1938). 
( We Slack, R. T. Lageman, and N. Underwood, Phys. Rev. 54, 358 

1 ’ 











The natural rotations are thus seen to be displaced in 
the negative (i.e. ‘“abnormal”) direction at the low tem- 
perature. The magnetic rotations, although in the positive 
(diamagnetic) direction, are increased at the low tem- 
perature, as in the presumably parallel case of NiSiF.-6H:0 
which has recently been observed and discussed by 
Becquerel and collaborators.? The two values just given 
for the Verdet constant at 5461A can be represented by the 
relation: : 

V=0.018+2.3/T, 


which agrees qualitatively with the earlier observations 
(see footnote 2 to Table I) made over a much more re- 
stricted temperature range. 
We acknowledge gratefully the assistance of the Amer- 
ican Philosophical Society in support of this investigation. 
Puitip RuDNICK 
F. G. SLack 
Joun O'Connor 
Vanderbilt University, 


Nashville, Tennessee, 
November 7, 1940. 


* Work conducted with the aid of a grant from the Penrose Fund of 
the American Philosophical Society. 

1L. R. Ingersoll, P. Rudnick, F. G. Slack, and N. Underwood, Phys. 
Rev. 57, 1145 (1940). 

2J. Becquerel and J. Van Den Handel, Physica 6, 1034 (1939); J. 
Becquerel and W. Opechowski, Physica 6, 1039 (1939). 
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Scattering of 20° Neutrons in Ortho- and Parahydrogen 


Since our last report on this subject,’ we have completely 
rebuilt our scattering chamber and have incorporated 
several new features in the neutron monochromator.? A 
rotating cadmium shutter was placed immediately in front 
of the BF; chamber to keep neutrons with energies above 
50°K from entering the detector. The scattering hydrogen 
was in the gas phase to eliminate the complicating effects of 
liquid forces. A jacket of liquid hydrogen boiling at 
atmospheric pressure maintained the gas at 20.4°K. The 
scattering chamber was 10 cm in diameter and 46 cm long; 
it was placed midway between source and detector, which 
were seven meters apart. The intensity of the beam was 
monitored by passing it alternately through the scattering 
chamber and through a dummy chamber. The ratio of 
neutron transmission through these two chambers was 
measured for 20°K neutrons to be 1.00+0.01. 

Fifty thousand counts were recorded during the ex- 
periment and the probable errors computed from the 
fluctuations agreed well with those expected from the total 
number of counts. The linearity of counting rate with 
cyclotron beam intensity was checked before and after each 
run, and all parts of the experimental set-up worked ex- 
ceedingly well throughout the course of the experiment. 
This fact is worth noting in view of the complex nature of 
the apparatus. The neutron temperature calculated from 
boron absorption measurements, assuming a 1/v depend- 
ence, agreed exactly with that obtained from the velocity 
measurements, i.e., 20.2°K. (This is in sharp disagreement 
with the work of Fertel, Gibbs, Moon, Thomson and Wynn- 
Williams,* who were unable to observe a change in the 
boron absorption cross section with neutron velocity.) 

The parahydrogen was prepared by allowing liquid 
hydrogen to remain in contact with active charcoal for 36 
hours. The pressure difference between the charcoal cham- 
ber and a surrounding liquid hydrogen bath was observed 
at frequent intervals. This indicated that the equilibrium 
state was approached very closely. The parahydrogen was 
then vaporized and the gas passed immediately into the 
scattering chamber where it was maintained at 20.4°K for 
the duration of the experiment. We believe that the 
composition of this sample was certainly within 0.2 percent 
of the equilibrium value, 99.8 percent parahydrogen. 

With one atmosphere of n-hydrogen in the scattering 
chamber J/J,>=0.238+0.005. With an equal pressure of 
parahydrogen, I/J>=0.868+0.01. The corresponding mo- 
lecular cross sections are oo= 103 X 10- cm? and ¢,=7.40 
x10 cm*. To obtain the true scattering cross sections, 
the capture cross section of two protons per molecule must 
be subtracted. If one assumes a 1/v law, this is 2.24 10-™* 
cm* at 20°K. The scattering cross sections for 20°K 
neutrons in ortho- and parahydrogen are then 


oo= (100+3) X10-* cm? 
op = (5.2+0.6) X 10-* cm? 
oo/o,=19. 


The earlier work in this field*~* gave values of o/o, from 
two to five. It may be concluded, as was suggested by 
those who performed these experiments that their ratios 











1004 


were reduced by excess scattering of faster neutrons in 
parahydrogen. The theoretical implications of these data 
will be discussed in a companion note by Dr. Schwinger. 
We had planned to repeat the work, to improve the 
statistical accuracy and to search for possible systematic 
errors, but pressure of other work now makes that im- 
possible for some time. 

Our thanks are due the Research Corporation for finan- 
cial assistance in this work. 
Luis W. ALVAREZ 


KENNETH S. PITZER 
Radiation Loy 5 


a be Physics an and Chemistry, 


Berkeley, California, © 
November 8, 1940. 
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Neutron Scattering in Ortho- and Parahydrogen and the 
Range of Nuclear Forces 


Experiments on slow neutron scattering in ortho- and 
parahydrogen' have proved useful in distinguishing be- 
tween the alternative possibilities of a real or a “virtual” 
singlet state of the deuteron. It is the purpose of this note 
to stress that experiments of the nature described by 
Alvarez and Pitzer? provide invaluable information con- 
cerning the neutron-proton interaction in the triplet state, 
in particular the effective range of the forces. The informa- 
tion ultimately available from such investigations is 
embodied in values for the singlet and triplet scattering 
amplitudes ao and a, which are related to the corresponding 
cross section, oo and 1, by oo,1=4aj,,. The sign of da» 
relative to that of a; is the experimental criterion for the 
nature of the singlet state. That the relative sign is negative 
is the essential result of previous investigation. The 
magnitude of a» coupled with a detailed assumption con- 
cerning the range and shape of the singlet interaction 
potential serves to determine the magnitude of the interac- 
tion energy. The magnitude of a; and the triplet binding 
energy will determine both magnitude and range of the 
triplet interaction energy with a definite assumption 
concerning the shape of the interaction potential. Thus the 
essence of the investigation is an accurate determination of 
the triplet scattering amplitude. 

The transition cross sections among the various levels of 
the hydrogen molecule involve different combinations of 
@; and dp. The cross sections for transitions among the para- 
levels are proportional to (3a;+4»)*, transitions among the 
ortho-levels are proportional to the combination (3a;+ 49)? 
+2(a,—a»)*, while transitions between the ortho- and 
para-level systems involve (2:—4»)?. Only the quantity 
(3a;+a0)? is sensitive to the precise numerical value 
of a;, for a» is large and of opposite sign compared to 
@;(a@9/a,;~—4). Parahydrogen at the low temperatures 
required for the nonexcitation of rotational levels, irradi- 
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ated by neutrons of energy insufficient for excitation js 
therefore the experimental medium appropriate to the 
determination of this quantity; for this enables one to 
measure the elastic scattering cross section of parahydrogen 
in its ground rotational level, free from contamination by 
other, large, cross sections, save for the small capture cross 
section. To avoid the disturbing influence of molecular 
binding forces, it is necessary to use molecular hydrogen 
in the gaseous phase. This, however, introduces a compli- 
cation occasioned by the large thermal motion of the 
molecules. If the target molecules are not stationary, but 
possess a distribution in the molecular velocity u, the 
effective cross section for a neutron of velocity v, cor¢(v), 
is not the usual cross section o(v), but rather 


oett(¥) =*S|v—ule(v—u) M(u)(du). 


Here N(u), the molecular velocity distribution function, 
represents the Maxwellian distribution appropriate to the 
temperature of the gas, 7. Thus, were the actual cross 
section independent of velocity,’ the effective cross section 
would exhibit a velocity dependence represented by‘ 


catE)=e{ (455 Jo@tetter|, g=(2E/kT)) 


2? 
considered as a function of the neutron energy E. Under 
the conditions employed by Alvarez and Pitzer (E/k=T 
=20°K), this yields o.f;=1.247¢. The actual correction 
factor for parahydrogen will be only somewhat smaller, 
since the cross section decreases but slightly in the small 
energy range under consideration (E =0.00173 ev). 

The cross section curves necessary for the evaluation of 
cetr have been computed by Schwinger and Teller,’ with 
extensions and improvements by Hamermesh.* The 
effective cross section for the 0—0 rotational transition, the 
elastic scattering of para Hz, is simply 

Cpara = F040 = 6.473 (3a, +p). 


The cross section for scattering by orthohydrogen will 
consist additively of the elastic scattering cross section 
(1-1) and the cross section for the inelastic process in 
which the molecule is converted to parahydrogen in its 
ground state (1—>0). For these, one obtains 


O11 = 6.291 [ (3a; +4d)?+2(a,— ao)*] 
T1»9 = 1 .447(a; — do)’. 


Combined with the experimental results of Alvarez and 
Pitzer, 
Tpara= (5.2 +0.6) x 10-* cm?, 
Fortho = (100+3) X10 cm, 

these formulae imply that 

(3a:+ao)* =0.803 X10-* cm*, (a;—a9)?=6.77 X10™ cm’. 
These equations have four sets of solutions but may, by the 
theoretical prediction that a, be of negative sign, be 
restricted to the two sets: 


(I) a=—4.26X10-"% cm, ao=21.8X10-" cm, 
(II) a:=—8.74X10-% cm, a@o=17.3X10-" cm. 


Had we used the lower limit of the experimental para cross 
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section (4.6 10-* cm?), the triplet amplitude of the first 
set would have been —4.49X10-" cm. The theoretical 
triplet amplitude increases in magnitude with increasing 
range and has the value 


—1/a( = (h?/M| Eo|)4) = —4.37X 10-" cm 


for zero range. Thus to within the statistical uncertainty 
of the observations the two admissible values of a; are 
(I) a:=—1/a, (Il) a:=—2/a. These correspond, respec- 
tively, to (I): zero range, (II): a rectangular-well range 
~8X10-" cm. Both values are theoretically untenable, the 
expected range being ~2.8X10-" cm. A further difficulty 
is revealed by computing the cross section for slow neutron 
scattering by free protons, viz.: 


301+ 400, 
which may be written 
da[ (3a; +40)?+3 (a: —ao)*] = 16.6 X 10- cm’. 


This result is considerably less than the directly measured 
value’ of (20+1) x 10 cm?. 

The consequences of these ortho-para measurements are 
at such complete variance with present theoretical concepts 
that it would be highly desirable to repeat these measure- 
ments and search for possible systematic errors. 

JULIAN SCHWINGER* 


Department of Physics, 
University of California, 
Berkeley, California, 
November 14, 1940. 


1F. G. Brickwedde, J. R. Dunning, H. J. Hoge, and J. H. Manley, 
Phys. Rev. 54, 266 (1938). 

?L. W. Alvarez and K. S. Pitzer, Cf. accompanying note. 

+ This situation is realized in the scattering of slow neutrons by 
helium. The cross section (1.5 X10-* cm*) obtained by H. Carroll and 
J. R. Dunning, Phys. Rev. 54, 541 (1938), is not the true cross section 
o but oer averaged with respect to the thermal energy distribution of the 
neutrons. The corrected value is ¢ =1.25 X10™™ cm?, 

4 (g) =2e4 fite-#*dx. 

‘ yo and E. Teller, Phys. Rev. 52, 286 (1937). 

*M. Hamermesh and J. Schwinger, Hg Rev. 55, 679 (1939). 

7V. W. Cohen, H. H. Goidsmith and J. Schwinger, Phys. Rev. 55, 106 
(1939); H. B. Hanstein, Phys. Rev. 57, 1045 (1940). 

* This work was done while the author held a National Research 
Fellowship, 1939-1940. 





Spectroscopically Pure Mercury (198) 


For some time, spectroscopists in the national labora- 
tories have been searching without much success for a line 
more nearly monochromatic than the red line of cadmium, 
to use as a standard of length. Professor W. E. Williams 
pointed out to us that if it were ever possible by some 
means to separate the isotopes of mercury, the green line 
45461 produced by one of the even isotopes would be 
admirably suited for the purpose. There would be no 
hyperfine structure, no isotope shift, and little Doppler 
broadening because of the high mass. 

We have bombarded gold with slow neutrons from the 
60” cyclotron, and have collected enough of the trans- 
mutation product, mercury, to observe its spectrum. Since 
gold has only one isotope, 197, slow neutron capture gives 
rise to a single radioactive isotope, Au'®*. This artificially 
radioactive product emits negative beta-rays with a half- 
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Fic. 1. Microphotometer traces for \4047 from ordinary Hg and 
rom Hg!%, 


life of 2.7 days and therefore turns into Hg'**, one of the 
stable isotopes of Hg. The experimental procedure is as 
follows: 

A cylinder of gold 15 cm long, 2.5 cm in diameter, and 
with a wall thickness of 0.2 mm is placed in a quartz tube 
of slightly larger diameter. To one end of this tube is fused 
a quartz capillary with an inside diameter of 2 mm. The 
whole system is evacuated and heated for 36 hours in a 
furnace almost to the melting point of gold. The gold is 
thus freed of any ordinary mercury contamination. 
Spectroscopically pure argon is then admitted to a pressure 
of 6 mm of Hg, and the quartz system is sealed off from the 
pumps. The gold cylinder in its quartz container is now 
placed in a paraffin-lined box near the target of the 
cyclotron, where it is bombarded with “stray” neutrons for 
about a month. At the end of this time the gold is again 
heated, while the end of the capillary tube is cooled in 
liquid air. After an hour of this treatment, a 3-cm length of 
the cooled capillary is sealed off. When the spectrum of the 
gas in this tube is excited by a 3-meter oscillator, the 
mercury lines are quite brilliant, but the argon spectrum is 
quenched. The mercury lines are visible after a neutron 
bombardment of a few hours, but they last for only a few 
seconds; under these conditions the Hg vapor is driven into 
the walls by the discharge. With a bombardment of a 
month, however, equilibrium between gas space and walls is 
apparently attained, so the spectrum is visible for some 
time. A microphotometer trace of a Fabry-Perot etalon 
spectrogram of the line \4047 is shown in Fig. 1. The 
absence of the hyperfine components shows that the 
mercury is actually a transmutation product. 

Since the Hg"* is a by-product of bombardments for 
biological purposes, no expenditure of ‘‘cyclotron time” is 
involved in its preparation. We will therefore be able to 
satisfy a reasonable demand for tubes filled with pure 
Hg", and we invite requests for such tubes. We gratefully 
acknowledge the support given to this work by the 
Research Corporation. 

Jacos WIENs 
Luts W. ALVAREZ 


Radiation Laboratory, 
nt of Physics, 
University of California, 
Berkeley, ifornia, 
November 9, 1940. 
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On the Nature of the Nuclear Forces 


It was shown recently by one of us,' that the properties 
of the wave equation for a particle with spin 1 in a Coulomb 
field are essentially different from those of the wave equa- 
tion for particles with spin 0 or }. It turns out that the 
solutions of this equation corresponding to the total 
momenta j21 and orbital momenta /=j+1 have an 
essential singularity in the coordinate origin. The wave 
equations are in the vicinity of the origin of the form 


v" +f(r)y=0 (1) 


with f~(e/yc*)(1/r*) (u is the mass of the particle). 

However, it is necessary to hold in mind that the 
applicability of the Proca equations is limited by distances 
of the order ro=e*/uc?, as was shown by one of us.’ If 
we note that (1) is of the form analogous to that of the 
Schrédinger equation the qualitative picture of the motion 
can be formally described as the motion of a fictitious non- 
relativistic particle with mass 1 in a potential hole with 
the radius ro and of the depth U~h?/r,?. It is known that 
the levels corresponding to the motion of the particle 
inside the hole exist if U > *h?/8r,?. As in the present case 

U is itself of the order h?/ro?, it is a priori impossible to 
decide if there are states which correspond to the motion 
limited by the distances of the order ro. In particular such 
states can prove to be possible only for certain values of 
the momentum j. 

If we identify, as is commonly done, mesotrons with 
particles with spin 1, and assume that the neutron consists 
of a proton and a negative mesotron, then we can conclude 
from the foregoing considerations that there is no specific 
interaction between the proton and the mesotron and the 
formation of the neutron is the result of the ordinary elec- 
trical interaction. If such a point of view is adopted we 
must assume that the system proton + negative mesotron 
possesses for the total momentum of the system equal to 
4 and the total momentum of the mesotron j=1 a level 
which corresponds to the relative motion of both particles 
at the distances of the order ro. 

So we are led to the view that the nuclear forces are to 
be regarded as originating from the electrical interaction 
of the nuclear particles by means of the mesotrons, and 
are connected with the peculiar properties of the equations 
of motion of the mesotrons. 

It could appear that such a point of view is inconsistent 
with the symmetry in the properties of the neutrons and 
the protons, which is substantially supported by experi- 
ment. However, it follows from the mere fact that the 
system proton + negative mesotron possesses a level with 
the binding energy ca. uc*, that the system neutron + posi- 
tive electron forms a proton with the liberation of nearly 
the same energy uc? (the possibility of the mesotron pair 
creation must be held in mind). 

As the radius e*/uc* of the proton and of the neutron is, 
in this theory, small as compared with that of the nucleus, 
the possibility arises of calculating various effects (the 
scattering of the mesotrons on the nuclei, the mesotron 
photo-effect, etc.) without going into the detailed con- 
sideration of the unknown interaction inside the protons 
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and the neutrons themselves (as is done in the theory of 
the deuteron).’ 

A detailed discussion of these problems will be published 
in the Journal of Physics (edited by the Academy of 
Sciences of the U. S. S. R.). 

L. Lanpau 

The Institute for Physical Problems 

of the Academy of Sciences of the U. S. S. R., 
I. TAMM 


The Physical Institute of the 
Academy of Sciences of the U. S. S. R., 
Moscow, U.S. S. R., 
September 25, 1940. 


11, Tamm, Phys. Rev. 58, 952 (1940). 
2? L. Landau, Phys. Rev. 58, 1006 (1940). 
3H. Bethe and R. Peierls, Proc. Roy. Soc. A148, 146 (1935). 





On the “Radius” of the Elementary Particles 


The limit e*/me* of the applicability of classical electro- 
dynamics can be deduced most convincingly from the con- 
dition that the reaction of the field of the electron must be 
small.! But the real limit of the classical electrodynamics 
is, as is well known, \>>h/me which lies much higher 
than A\~e?/mc*. As for the limits of applicability of the 
electrodynamical conceptions in quantum mechanics 
there are no reasons to assume that they are given by 
e?/mc?, and some other way must be found to determine 
the real “radius” of elementary particles. 

The previous discussions of this problem have always 
included some additional physical assumptions,* which 
makes the results quite unconvincing. It is our purpose 
here to show that in the quantum theory a simple method 
can be developed to find the limit beyond which the 
theory comes into contradiction with itself because of the 
neglect of the reaction of the field. 

The neglect of the reaction of the field means that we 
assume the particle to move under the influence of the 
external field without taking into account the variation of 
this field by the electron itself. In the mathematical treat- 
ment this is expressed by calculating the matrix elements 
of the interaction of the charged particle and a photon 
with the help of free particle wave functions. Hence the 
smallness of the reaction force corresponds in the quantum 
theory to the applicability of the perturbation theory to 
the collision of a charged particle and a photon. 

The effective cross section for the scattering of particles 
with an angular momentum / in the frame of reference 
moving with the center of inertia is o;=4(2/+1)xX? sin? 8), 
X being the wave-length (divided by 27), 5:- the phase 
of the wave function at infinity. The applicability of the 
perturbation theory requires 5<1. As the scattering 
decreases with increasing /, it suffices to require this con- 
dition to be fulfilled for /~1; the corresponding o; we denote 
as oo. Thus we come finally to the condition 


os? (1) 


(it is to be noted that the condition o; <4(2/+1)#X? must 
be fulfilled quite independently of the applicability of the 
perturbation theory; the opposite would mean that the 
number of scattered particles exceeds the number of cor- 
responding incident particles, which is impossible). 
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the correct results. As in the classical theory all the scat- 
tering is dipole scattering, we must simply insert for oo 
the full Thomson’s scattering cross section, this leading to 
e/me<x. 

We now apply (1) to the quantum theory of the electron. 
We must first of all transform the Klein-Nishina formula 
to the proper frame of reference and then take part of the 
scattering corresponding to the angular momentum of the 
order 1. This gives oo~(e?/mc*)*mc?/E (E being the energy 
of the photon in the frame of reference where the electron 
is at rest). oo differs from the total cross section by the 
absence of the factor log E/mc* which arises from the 
scattering with large angular momenta. The energy e of 
the photon in the frame of reference moving with the center 
of inertia is given by 2@=mcE (when E>mc) so that 
X~hc/e~hc/(meE)'. The energy E disappears now from 
(1) and we get e*<he which is always fulfilled. This means 
that in the quantum theory of the electron there are no 
limits of its applicability, arising from itself. The “radius” 
of the electron in quantum electrodynamics is in some 
sense equal to zero. 

One gets the same result for particles with zero spin, if 
the corresponding expression for Compton effect‘ is used. 

The Compton effect for particles with spin 1 (mass ») 
was recently calculated by F. Booth and A. H. Wilson® 
and by J. Smorodinski.* The result is co~(e?/uc*?)*E/yc*. 
Hence the condition (1) gives E<ycthc/e*. The “radius” 
of the particle is equal to the corresponding wave-length 
in the frame of reference in which the particle is at rest, i.e. 


r~he/E~e*/ue*. (2) 


Hence the “radius” of the particle with spin 1 is equal to 
its classical value. 

A more detailed discussion of the problem with the 
application to the scattering of particles with spin 1 by a 
charged particle, is to be published in Journal of Physics, 
volume 3, 1940. 


L. LANDAU 


The Institute for Physical Problems 
of the Academy of Sciences of U. S. S. R., 
Moscow, U.S. S. 
September 19, 1940. 


1 Cf. W. Heitler, The, Quantum Theory of Radiation (Oxford, 1936), 
33, 232. 
PP; 2W. Heisenberg, Zeits. f. Physik vO = (1938). 
3 V. Weisskopf, Phys. Rev. 56, 73 (1939). 
4W. Gordon, Zeits. f. Physik 40, 117 (1927). 
5‘ F. Booth, and A. H. W ilson, Proc. Roy. Soc. A175, 483 (1940). 
6 J. Smorodinski, J. Phys. 3 (1940). 





Relative Intensities of 'S—*P and 'S—'!P 
Transitions in Mg 


The relative intensity of singlet-triplet and singlet- 


singlet transitions in cases where the spin-orbit interaction 
is weak is given by 


f= |. S3|2(vs/mi)fi, (1) 


where the subscripts 3 and 1 denote singlet-triplet and 
singlet singlet values, respectively. The factor S; is the 
coefficient which represents the contamination of the 
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It is easy to see that in the classical theory (1) leads to 
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singlet state by the triplet state; that is, approximately 
¥(P) = fo('P)+Ssvo(*P), (2) 


where the y's are Russell-Saunders wave functions. One 
has from perturbation theory 


| Ss! = H(ij)/hvji, (3) 


where H(ij) is the matrix element of the Hamiltonian. 
Furthermore, the deviation of the central component of the 
triplet from the interval rule is 


AW =(A(aj) P/hvji. (4) 
According to Houston,! we may take for P states 
HA(ij) =(1(1+1)}\A = v2A, (S) 


if A is the constant of the spin-orbit interaction AL-S. In 
practice A is determined from the over-all width of the 
triplet, 3A. The theoretical expression for A is 


h? 1aV 
~igeaad Rie(* ala "dr. (6) 


The radial wave function, R3p, of the triplet state is here 
assumed to be identical with that of the singlet state. 
Since the relative intensities found from (1), (3), and (5) 
are in only fair agreement with observation, King and Van 
Vleck? have suggested that an undetermined parameter \ 
be inserted in (5) and (6) to account for the difference 
between R3p and Rip which was previously neglected. 


Then A refers to the integral (6) with Rip in the integrand 
while \A denotes the integral with R3pRip. The value of 
is determined so as to make the deviation (4) from the 
interval rule correct. In this way King and Van Vleck* 
find for the atoms Hg, Cd, Zn, Ba, Sr, and Ca that \~0.8. 
If, however, we try to apply their procedure to the transi- 
tion 3s?!'\S—3s3p*P of Mg, we find \=2.4, a value out of 
line with those previously found. It seems, then, that in 
this case the deviations from the interval rule are not due 
simply to perturbations by 'P but that other effects enter. 
Mg is so light that spin-spin interaction may contribute 
appreciably to the observed deviation, 0.3 cm=', from the 
interval rule. One should therefore not place much trust in 
the value of \ for Mg derived above. 

It is thus uncertain whether to use the value \=0.8 by 
analogy with previous calculations, or merely to use the 
Houston formula, where \=1. The relative intensities 
found in the two cases are 


fi/fs=340,000 for A=1, (7 
fi/fs=530,000 for =0.8. ) 


We may consider these values as correct within a factor 2 
or 3. 

I am indebted to Professor D. H. Menzel for suggesting 
this problem and Professor J. H. Van Vleck for his guidance 
and advice. 


PEARL J. RUBENSTEIN 
Harvard College Observatory, 
Cambridge, Massachusetts, 
November 8, 1940 


1W. V. Houston, Phys. Rev. 33, 297 (1929). 
2G. W. King and J. H. Van Vleck, Phys. Rev. 56, 464 (1939). 
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Gamma-Radiation from Rhenium 


We have produced by proton bombardment of tungsten 
(purity 99.95 percent) the long period rhenium isotope 
which Fajans and Sullivan! obtained by W(D,n) and 
Re(n,2m) and assigned to Re'. We have also produced 
the 90-hour and 30-minute periods by proton bombard- 
ment and checked the fact that they follow rhenium 
chemistry. 

After removal from the cyclotron, the tungsten strip 
which had been bombarded intermittently for four months 
with currents of the order of 20 microamperes, was allowed 
to age for two months before observations began, so that 
the shorter periods might die out, after which the activity 
was followed for four months. The original intensity of the 
long period was about 100 divisions per minute on the 
Lauritsen electroscope. Our value of the half-life is 54+2 
days. See Fig. 1. 

Absorption measurements of the radiation showed 
negative electrons of three kinds, of ranges 12, 50, and 350 
mg/cm? of aluminum, corresponding to energies of roughly 
0.1, 0.22, and 0.86 Mev, respectively, and gamma-radiation 
of about 1 Mev, as seen from Fig. 2. The highest energy 
group may be conversion electrons from this gamma-ray, 
which was also reported by the above authors.' The inter- 
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' 
0.096 Mev 
Fic. 3. Conversion electrons from tungsten bombarded with protons, 


mediate energy electrons were too few for careful study, 
and indeed may be from a tungsten isotope produced from 
the_neutrons around the cyclotron. 

Previous experience had shown that a beta-ray spectro- 
graph could be used to photograph a line of electrons of 
around 100 kev energy if the intensity was such that the 
equivalent of about 10° divisions of the electroscope were 
available. This corresponds roughly to 10° particles. Since 
the fraction of these striking the photographic film in our 
spectrograph is about 10-*, as an estimate 10° electrons are 
required to produce a usable image. Because various 
workers? have shown the reciprocity law to be valid for 
electrons, it seemed worth while to attempt a long exposure 
to measure the low energy group from our target. A twenty- 
one-day exposure showed a line at 0.096 Mev, reproduced 
in Fig. 3. If these are K-conversion electrons, there must 
be a gamma-ray from the 54-day rhenium isotope of 
0.096+0.071 (Re K edge) =0.17 Mev. 

E. Creutz 
W. H. Barkas* 


N. H. FurmMant 


Palmer Physical Laboratory, 
Princeton University, 
Princeton, New Jersey, 
November 13, 1940. 


* Now with Bureau of Ordnance, U. S. Navy, Thirteenth Naval 
District, Seattle, Washington. 

+ Department of Chemistry, Princeton University. 

1K. Fajans and W. H. Sullivan, Phys. Rev. 58, 276 (1940). 

2 A. Becker and E. Kipphan, Ann. d. Physik 10, 15 (1931); also M. J. 
Nacken, Physik. Zeits. 31, 296 (1930). 





The Electric Fields in Vibrating Polar Crystals 


Within a vibrating polar crystal there are electric forces 
caused by the polarization, P, associated with the dis- 
placements of the ions of the crystal during the vibration. 
These forces give rise to an electric field, E, and an electric 
displacement, D, which, according to Lyddane and 
Herzfeld,' are connected with the polarization by 


=—42P, D=0, (1) 


if the wave-length of the vibration is short compared to 
the size of the crystal but long compared to the lattice 
distance. This result is stated to be independent of whether 
the waves in the crystal are longitudinal or transverse. On 
the other hand, Fréhlich and Mott? have given formulae 
that agree with (1) for the longitudinal displacement but 


become 
E=0, D=47P (2) 


for the transverse case. 

That the latter result is the correct one becomes apparent 
from the following direct calculation based on formulae 
given by Lyddane and Herzfeld! for the electric force, F, 
within a crystal. When the wave-length is long compared 
to the dimension of a crystal slab, these formulae are: 
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(1a) F=—(8/3)P if P is perpendicular to the face of 
the slab, 
(1b) F=(42/3)P if P is parallel to the face of the slab. 


Should the wave-length be short compared to the dimen- 
sions of the slab but long compared to the lattice distance: 


(2a) F=-—(8/3)P for the longitudinal vibration 
(2b) F=(4x/3)P for the transverse vibration. 


According to the conventional definition, the electric 
field is the electric force in a slot in the crystal with its 
short dimension perpendicular to the direction of the 
electric force. We consider a slot of dimensions short com- 
pared to the wave-length. Then E is obtained by subtract- 
ing from (2a) or (2b) the effect of the material that has 
been removed to form the slot. Since the electric force is 
always parallel to P, the effect of this material is just the 
force (1b). Therefore 


E = —4nxP for the longitudinal vibration, 
E=0 for the transverse. 


The displacement, D, is given by the electric force in the 
slot when its short dimension is parallel to the polarization. 
Therefore, subtracting (la) as the effect of the material in 
the slot from (2a) and (2b) 


D=0, D=4xP 


for the longitudinal and transverse vibrations, respectively. 

These formulae are consistent with D=E+4-xP, and 
they also lead, according to (2a) and (2b), to the Lorentz- 
Lorenz force for both types of vibration; that is: 


F=E-+ (4x/3)P. 


At first sight, the fact that the electric displacement does 
not vanish for the transverse vibration may appear to be 
inconsistent with the absence of an impressed external 
field. However, the polarized crystal itself gives rise to an 
external field which is just 4%P at the surface and vanishes 
exponentially with distance from the surface. This may be 
seen by considering a charge distribution on a plane 
surface that is periodic in one direction with wave-length 
\ and constant in the perpendicular direction. Such a 
charge distribution gives rise to a field proportional to- 
exp (—22x/A), if x is perpendicular to the surface. A 
transversely displaced crystal with the polarization per- 
pendicular to the face of the slab may be considered as 
made up of a series of such planes, and it is apparent that 
there will be an exponentially decreasing external field 
over a distance approximately given by the wave-length 
of the vibration under consideration. The result obtained 
directly from the lattice sums by Kellerman is in agree- 
ment with this qualitative argument. 

R. H. LyppDANE 


University of North Carolina, 
Chapel Hill, North Carolina. 
K. F. HERZFELD 
Catholic University, 
Washington, D. C. 
R. G. Sacus 
George Washington University, 
Washington, D. C., 


November 1, 1940. 
1R. H. Lyddane and K. F. Herzfeld, Phys. Rev. 54, 846 (1938). 


? H. Frohlich and N. F. Mott, Proc. Roy. Soc. A171, 496 (1939). 
+E. W. Kellerman, Phil. Trans. Roy. Soc. A238, 513 (1940). 
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The C“H,— N" Mass Difference 


Three years ago the author, in an attempt to obtain 
more accurate mass values, undertook to construct a new 
type of mass spectrograph having a very high resolving 
power and dispersion. This instrument was completed last 
spring and reported' at the Washington meeting of the 
American Physical Society. ! 

Its constants are given in Table I together with those of 
other instruments being used at the present time to study 
the light element doublets. 

The mass scale is correction free, mass differences being 
obtained directly by taking the mean of the ratios of the 
doublet separation to the forward and backward dispersion 
line separations. In the case of the C"H:—N™*™ doublet, 
this method of calculation gives the mass difference directly 
to 0.000015 mass unit, a figure which is several times 
smaller than the error introduced in making a single 
doublet distance measurement. In addition, exact methods 
have been worked out for correcting this small deviation 
and it is estimated that the resulting error due to the 
method of computation alone is approximately three to 
four units in the sixth decimal place. The above relation 
holds true regardless of the plate position within rather 
wide limits. 

Recently a considerable number of matched as well as 
unmatched C®H,—N*" doublets were obtained in a series 
of spectra which were photographed under different con- 
ditions and used to verify the theoretically expected rela- 
tions. The values of the mass differences corresponding to 
these doublets are listed in Table II in order that a birds- 
eye view of the variations may be obtained. These data 
result from measurements taken on three different plates, 
three to five measurements of the doublet distance being 
made in each instance. Although the difference in density 
between the doublet lines is not large in any case, only 
about half the doublets corresponding to the values listed 
are matched. It so happens that if the matched doublets 
alone are considered, the same average value is obtained, 
the probable error also being approximately the same. 

For comparison purposes, this mass difference, together 
with other published values for the light element doublets 
is listed in Table III. Although no cross check doublets 


TABLE I. Mass-spectrograph constants. 








DISPERSION (MM 
FOR 1% MASS DIFF.) 


RESOLVING 


MASS SPECTROGRAPH Power M/AM 





Aston 2,000 4.5 
Mattauch 6,500 1.4 
Asada and others 17,000 4.65 - 
Jordan 30,000 14.6 











TABLE II. The C"H2—N"™ doublel mass differences. 








0.01255 0.01255 0.01257 0.01258 0.01254 
0.01258 0.01254 0.01263 0.01254 0.01255 
0.01256 0.01252 0.01258 0.01257 0.01251 
0.01253 0.01260 0.01256 0.01265 0.01252 
0.01254 0.01256 

Average 0.01256 +0.000015* 














_ * This error is three times the probable error computed from the 
internal consistency of the data. 
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TABLE III. Mass differences in terms of 1074 mass unit. 











DouBLeEt H? —D? | D? —C#++! C&Hy —O1%8 | CY Hy’ —N“ 
Aston 15.2 423.6 360.1 124.5 
+0.4 +1.8 +1.6 +0.7 
Bainbridge 15.3 421.9 364.9 127.4 
and Jordan +0.4 +0.5 +0.8 +1.1 
Mattauch 15.39 422.39 364.06 125.81 
+0.071 +0.21 +0.40 +0.23 
Asada 364.2 125.7 
and Others +0.9 +0.6 
Jordan 125.6 
+0.15 




















have been measured as yet, it seemed worth while to 
publish this value in view of the fact that it was obtained 
on an instrument which not only has extremely high 
constants but also differs radically from any of the others 
in use at the present time. Other measurements on the 
light element doublets are in progress. 

E. B. JorDAN 


Physics Department, 
University of Illinois, 
Urbana, Illinois, 
28, 1940. 


1 E. B. Jordan, Phys. Rev. 57, 1072A (1940). 





Cosmic Rays and Comets 


The conjecture that the comets are contraterrene bodies,' 
which accounts, if only in a rough qualitative way, for 
some of the more conspicuous features of comets, can be 
properly tested as follows: (a) by investigating its useful- 
ness in detailed quantitative interpretation of the now 
known properties of comets; (b) by performing experi- 
ments, involving cosmic-ray measurements, during meteor 
showers known to be of cometary origin; and (c) by 
performing appropriate cosmic-ray experiments while a 
sufficiently active comet is sufficiently near to the earth. 
Method (a) is indirect ; and since it involves detailed study 
of the intricacies of cometary behavior, as well as a 
knowledge of the distribution of meteoric material within 
the solar system, it is difficult. Method (b) implies the 
acceptance of the idea that certain meteor showers are of 
cometary origin, but is otherwise quite direct; and the 
requisite experiments can be done, for example, almost 
every August. Method (c) is in principle the most direct ; 
but the occasions on which it can be tried are very rare. 

One purpose of this note is to call attention to the fact 
that if the comet recently discovered by Cunningham 
proves sufficiently active, it might, within the next few 
weeks, provide opportunities for performing the experi- 
ments involved in method (c). Lacking the appropriate 
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data, however, the writer is at present unable to venture an 
estimate whether, if comets are contraterrene, the cosmic- 
ray effects of Cunningham’s comet would be sufficiently 
intense for experimental purposes. 

If comets are contraterrene, their activity depends not 
only on the solar energy that they receive, but also on the 
rate at which they encounter ordinary meteoric matter. In 
the absence of a theory of interaction of terrene and 
contraterrene matter, it is impossible to say with certainty 
what the primary products of their mutual annihilation 
would be. It is reasonably safe to guess, however, that 
among these products there should be photons having 
energies of about a billion electron volts. Therefore 
experiments designed to determine whether any high 
energy photons do originate in a comet should best be 
performed in the upper atmosphere. 

The second purpose of this note is to mention that among 
the products of mutual annihilation of terrene and contra- 
terrene matter there may be free mesons. This possibility 
is of interest not only in view of the present conjecture 
concerning comets, but also in connection with cosmic rays 
in general; for it suggests that: (1) The continual release of 
mesons at the top of the earth’s atmosphere may be caused 
by contraterrene matter coming from interstellar space and 
impinging upon the atmosphere, and (2) the non-ionizing 
particles which release mesons deeper in the atmosphere 
may be contraterrene neutrons, which have themselves 
been set free in the earlier stages of the annihilation of the 
impinging contraterrene atomic nuclei.? 

V. ROJANSKY 

Union College, 

Schenectady, New York, 
November 4, 1940. 


1V. Rojansky, Astrophys. J. 91, 257 (1940). We call ‘‘contraterrene"’ 
a body composed of hypothetical atoms consisting of negatively charged 
nuclei surrounded by positrons; an ordinary body we call ‘‘terrene.”’ 

2 The possibility of production of high energy photons in the upper 
atmosphere by the annihilation of negative protons arriving from 
interstellar s was discu by F. Zwicky, Phys. Rev. 48, 169 
— > also F. Zwicky, Proc. Nat. Acad. Sci. 22, 266 (1936), esp. 
Section F. 





Erratum: On the Angular Distribution of Fast Neutrons 
Scattered by Hydrogen, Deuterium and Helium 


(Phys. Rev. 58, 590 (1940)) 


The captions of Fig. 3 and Fig. 4 on page 592 should be 
interchanged. 
H. H. BARSCHALL 
M. H. KANNER 


Palmer Physical Laboratory, 
Princeton University, 
Princeton, N. J., 
November 15, 1940. 
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